Lines Matching full:fjac
15 int fcn_chkder(const VectorXd &x, VectorXd &fvec, MatrixXd &fjac, int iflag)
49 fjac(i,0) = -1.;
50 fjac(i,1) = tmp1*tmp2/tmp4;
51 fjac(i,2) = tmp1*tmp3/tmp4;
62 MatrixXd fjac(m,n);
69 internal::chkder(x, fvec, fjac, xp, fvecp, 1, err);
70 fcn_chkder(x, fvec, fjac, 1);
71 fcn_chkder(x, fvec, fjac, 2);
72 fcn_chkder(xp, fvecp, fjac, 1);
73 internal::chkder(x, fvec, fjac, xp, fvecp, 2, err);
147 int df(const VectorXd &x, MatrixXd &fjac) const
156 fjac(i,0) = -1;
157 fjac(i,1) = tmp1*tmp2/tmp4;
158 fjac(i,2) = tmp1*tmp3/tmp4;
223 internal::covar(lm.fjac, lm.permutation.indices()); // TODO : move this as a function of lm
231 // std::cout << fjac*covfac << std::endl;
234 cov = covfac*lm.fjac.topLeftCorner<n,n>();
237 // VERIFY_IS_APPROX( covfac*fjac.topLeftCorner<n,n>() , cov_ref);
260 int df(const VectorXd &x, MatrixXd &fjac)
263 assert(fjac.rows()==n);
264 assert(fjac.cols()==n);
268 fjac(k,j) = 0.;
269 fjac(k,k) = 3.- 4.*x[k];
270 if (k) fjac(k,k-1) = -1.;
271 if (k != n-1) fjac(k,k+1) = -2.;
614 internal::covar(lm.fjac, lm.permutation.indices()); // TODO : move this as a function of lm
622 // std::cout << fjac*covfac << std::endl;
625 cov = covfac*lm.fjac.topLeftCorner<n,n>();
628 // VERIFY_IS_APPROX( covfac*fjac.topLeftCorner<n,n>() , cov_ref);
648 int df(const VectorXd &b, MatrixXd &fjac)
651 assert(fjac.rows()==54);
652 assert(fjac
657 fjac(i,0) = -x*e*factor;
658 fjac(i,1) = -e*factor*factor;
659 fjac(i,2) = -x*e*factor*factor;
732 int df(const VectorXd &b, MatrixXd &fjac)
735 assert(fjac.rows()==14);
736 assert(fjac.cols()==2);
738 fjac(i,0) = (1.-exp(-b[1]*m_x[i]));
739 fjac(i,1) = (b[0]*m_x[i]*exp(-b[1]*m_x[i]));
811 int df(const VectorXd &b, MatrixXd &fjac)
814 assert(fjac.rows()==236);
815 assert(fjac.cols()==7);
819 fjac(i,0) = 1.*fact;
820 fjac(i,1) = x*fact;
821 fjac(i,2) = xx*fact;
822 fjac(i,3) = xxx*fact;
824 fjac(i,4) = x*fact;
825 fjac(i,5) = xx*fact;
826 fjac(i,6) = xxx*fact;
903 int df(const VectorXd &b, MatrixXd &fjac)
906 assert(fjac.rows()==14);
907 assert(fjac.cols()==2);
910 fjac(i,0) = b[1]*x[i] / den;
911 fjac(i,1) = b[0]*x[i]*(den-b[1]*x[i])/den/den;
978 int df(const VectorXd &b, MatrixXd &fjac)
981 assert(fjac.rows()==24);
982 assert(fjac.cols()==6);
984 fjac(i,0) = exp(-b[1]*x[i]);
985 fjac(i,1) = -b[0]*x[i]*exp(-b[1]*x[i]);
986 fjac(i,2) = exp(-b[3]*x[i]);
987 fjac(i,3) = -b[2]*x[i]*exp(-b[3]*x[i]);
988 fjac(i,4) = exp(-b[5]*x[i]);
989 fjac(i,5) = -b[4]*x[i]*exp(-b[5]*x[i]);
1066 int df(const VectorXd &b, MatrixXd &fjac)
1069 assert(fjac.rows()==9);
1070 assert(fjac.cols()==3);
1073 fjac(i,0) = 1./(1.+e);
1074 fjac(i,1) = -b[0]*e/(1.+e)/(1.+e);
1075 fjac(i,2) = +b[0]*e*x[i]/(1.+e)/(1.+e);
1143 int df(const VectorXd &b, MatrixXd &fjac)
1146 assert(fjac.rows()==16);
1147 assert(fjac.cols()==3);
1151 fjac(i,0) = e;
1152 fjac(i,1) = b[0]*factor*e;
1153 fjac(i,2) = -b[1]*b[0]*factor*factor*e;
1222 int df(const VectorXd &b, MatrixXd &fjac)
1225 assert(fjac.rows()==6);
1226 assert(fjac.cols()==2);
1229 fjac(i,0) = 1.-e;
1230 fjac(i,1) = b[0]*x[i]*e;
1301 int df(const VectorXd &b, MatrixXd &fjac)
1304 assert(fjac.rows()==33);
1305 assert(fjac.cols()==5);
1307 fjac(i,0) = 1.;
1308 fjac(i,1) = exp(-b[3]*x[i]);
1309 fjac(i,2) = exp(-b[4]*x[i]);
1310 fjac(i,3) = -x[i]*b[1]*exp(-b[3]*x[i]);
1311 fjac(i,4) = -x[i]*b[2]*exp(-b[4]*x[i]);
1389 int df(const VectorXd &b, MatrixXd &fjac)
1392 assert(fjac.rows()==11);
1393 assert(fjac.cols()==4);
1397 fjac(i,0) = (xx+x*b[1]) * factor;
1398 fjac(i,1) = b[0]*x* factor;
1399 fjac(i,2) = - b[0]*(xx+x*b[1]) * x * factor * factor;
1400 fjac(i,3) = - b[0]*(xx+x*b[1]) * factor * factor;
1474 int df(const VectorXd &b, MatrixXd &fjac)
1477 assert(fjac.rows()==154);
1478 assert(fjac.cols()==3);
1481 fjac(i,0) = e;
1482 fjac(i,1) = - b[0]*e/b[2]/(b[1]+x[i]);
1483 fjac(i,2) = b[0]*e*log(b[1]+x[i])/b[2]/b[2];
1555 int df(const VectorXd &b, MatrixXd &fjac)
1558 assert(fjac.rows()==37);
1559 assert(fjac.cols()==7);
1563 fjac(i,0) = 1.*fact;
1564 fjac(i,1) = x*fact;
1565 fjac(i,2) = xx*fact;
1566 fjac(i,3) = xxx*fact;
1568 fjac(i,4) = x*fact;
1569 fjac(i,5) = xx*fact;
1570 fjac(i,6) = xxx*fact;
1651 int df(const VectorXd &b, MatrixXd &fjac)
1654 assert(fjac.rows()==15);
1655 assert(fjac.cols()==4);
1659 fjac(i,0) = pow(1.+e, power);
1660 fjac(i,1) = power*b[0]*e*pow(1.+e, power-1.);
1661 fjac(i,2) = -power*b[0]*e*x[i]*pow(1.+e, power-1.);
1662 fjac(i,3) = b[0]*power*power*log(1.+e)*pow(1.+e, power);
1739 int df(const VectorXd &b, MatrixXd &fjac)
1742 assert(fjac.rows()==35);
1743 assert(fjac.cols()==3);
1747 fjac(i,0) = e / b[1];
1748 fjac(i,1) = ((x[i]-b[2])*(x[i]-b[2])/b12-1.) * b[0]*e/b12;
1749 fjac(i,2) = (x[i]-b[2])*e*b[0]/b[1]/b12;
