lbfgs - problem
- Subject: lbfgs - problem
- From: Robert Gloeckner
- Date: Thu, 28 Feb 2008 03:08:41 +0100
Hello,
if i try to fit [2] function-parameters to my data-set [3], i get an
error [1].
I do not understand what the error message means - can somebody give me
a hint what is going wrong?
Thanks,
Bernd
/*------- [1] errors from maxima -------*/
Maxima encountered a Lisp error:
WRITE-CHAR on #<CLOSED OUTPUT BUFFERED FILE-STREAM CHARACTER
#P"/dev/fd/1"> is illegal
Automatically continuing.
To reenable the Lisp debugger set *debugger-hook* to nil.
/*------- [2] maxima-code -----------*/
kill(all);
/*__ Carreau */
eq1: eta = a_T * P_1 / (1 + a_T * P_2 - dtgamma )^P_3;
/*__ WLF */
eq2: log(a_T) = (8.86 * (P_4 - P_5) / (101.6 + P_4 - P_5)) - (8.86 * (T
- P_5) / (101.6 + T - P_5));
/*__ Arrenius */
eq3: a_T = exp( (P_4 / R) * ((1 / (T + 273.15)) - (1 / P_5)));
/*__ Carreau-WLF */
eq100: radcan(exp(lhs(eq2))) = radcan(exp(rhs(eq2)));
eq101: ev( eq1, eq100);
eq102: ev( eq101, R = 8.31415);
fortran( radcan( eq102));
eq103: radcan( eq102);
Carreau_WLF: eq103;
/*__ Carreau-Arrenius */
eq200: ev( eq1, eq3);
eq201: ev( eq200, R = 8.31415);
eq202: radcan( eq201);
Carreau_Arrenius: eq202;
/*__ Fitting */
load (descriptive);
load (numericalio);
load (lbfgs);
m : read_matrix("~/Sourcen/workdir/rheofit-cvs/rheofit-data-frank.dat");
F(dtgamma,T) := ev( rhs(Carreau_WLF));
/* F(dtgamma,T) := ev( rhs(Carreau_Arrenius)); */
/* data-file-format is m[i,1] = dtgamma, m[i, 2] = eta, m[i,3] = T */
FOM : '(sum ((F(m[i, 1],m[i ,3]) - m[i, 2])^2, i, 1, length(m)) /
length(m));
estimates : lbfgs (FOM, '[P_1,P_2,P_3,P_4,P_5], [1, 1, 1, 1, 1], 1e-2,
[1, 0]);
/*--------- [3] data-file [dtgamma, eta, T] ----------*/
882.13 151.17 200.0
436.70 267.73 200.0
199.64 508.15 200.0
84.61 1039.57 200.0
34.59 2125.71 200.0
14.08 4184.46 200.0
5.77 7733.52 200.0
2.47 12466.0 200.0
917.02 106.88 230.0
365.20 221.23 230.0
149.02 456.17 230.0
62.03 859.14 230.0
26.26 1484.19 230.0
11.33 2322.23 230.0
4.92 3403.47 230.0
2.14 4634.94 230.0
775.53 104.28 250.0
308.00 210.07 250.0
127.86 385.42 250.0
54.31 655.88 250.0
23.67 996.30 250.0
10.36 1400.29 250.0
4.61 1786.14 250.0
2.04 2139.91 250.0