Lines Matching refs:fjac
18 int fcn_chkder(const VectorXd &x, VectorXd &fvec, MatrixXd &fjac, int iflag) in fcn_chkder() argument
52 fjac(i,0) = -1.; in fcn_chkder()
53 fjac(i,1) = tmp1*tmp2/tmp4; in fcn_chkder()
54 fjac(i,2) = tmp1*tmp3/tmp4; in fcn_chkder()
65 MatrixXd fjac(m,n); in testChkder() local
72 internal::chkder(x, fvec, fjac, xp, fvecp, 1, err); in testChkder()
73 fcn_chkder(x, fvec, fjac, 1); in testChkder()
74 fcn_chkder(x, fvec, fjac, 2); in testChkder()
75 fcn_chkder(xp, fvecp, fjac, 1); in testChkder()
76 internal::chkder(x, fvec, fjac, xp, fvecp, 2, err); in testChkder()
150 int df(const VectorXd &x, MatrixXd &fjac) const in df()
159 fjac(i,0) = -1; in df()
160 fjac(i,1) = tmp1*tmp2/tmp4; in df()
161 fjac(i,2) = tmp1*tmp3/tmp4; in df()
226 internal::covar(lm.fjac, lm.permutation.indices()); // TODO : move this as a function of lm in testLmder()
237 cov = covfac*lm.fjac.topLeftCorner<n,n>(); in testLmder()
263 int df(const VectorXd &x, MatrixXd &fjac) in df()
266 assert(fjac.rows()==n); in df()
267 assert(fjac.cols()==n); in df()
271 fjac(k,j) = 0.; in df()
272 fjac(k,k) = 3.- 4.*x[k]; in df()
273 if (k) fjac(k,k-1) = -1.; in df()
274 if (k != n-1) fjac(k,k+1) = -2.; in df()
617 internal::covar(lm.fjac, lm.permutation.indices()); // TODO : move this as a function of lm in testLmdif()
628 cov = covfac*lm.fjac.topLeftCorner<n,n>(); in testLmdif()
651 int df(const VectorXd &b, MatrixXd &fjac) in df()
654 assert(fjac.rows()==54); in df()
655 assert(fjac.cols()==3); in df()
660 fjac(i,0) = -x*e*factor; in df()
661 fjac(i,1) = -e*factor*factor; in df()
662 fjac(i,2) = -x*e*factor*factor; in df()
735 int df(const VectorXd &b, MatrixXd &fjac) in df()
738 assert(fjac.rows()==14); in df()
739 assert(fjac.cols()==2); in df()
741 fjac(i,0) = (1.-exp(-b[1]*m_x[i])); in df()
742 fjac(i,1) = (b[0]*m_x[i]*exp(-b[1]*m_x[i])); in df()
816 int df(const VectorXd &b, MatrixXd &fjac) in df()
819 assert(fjac.rows()==236); in df()
820 assert(fjac.cols()==7); in df()
824 fjac(i,0) = 1.*fact; in df()
825 fjac(i,1) = x*fact; in df()
826 fjac(i,2) = xx*fact; in df()
827 fjac(i,3) = xxx*fact; in df()
829 fjac(i,4) = x*fact; in df()
830 fjac(i,5) = xx*fact; in df()
831 fjac(i,6) = xxx*fact; in df()
910 int df(const VectorXd &b, MatrixXd &fjac) in df()
913 assert(fjac.rows()==14); in df()
914 assert(fjac.cols()==2); in df()
917 fjac(i,0) = b[1]*x[i] / den; in df()
918 fjac(i,1) = b[0]*x[i]*(den-b[1]*x[i])/den/den; in df()
985 int df(const VectorXd &b, MatrixXd &fjac) in df()
988 assert(fjac.rows()==24); in df()
989 assert(fjac.cols()==6); in df()
991 fjac(i,0) = exp(-b[1]*x[i]); in df()
992 fjac(i,1) = -b[0]*x[i]*exp(-b[1]*x[i]); in df()
993 fjac(i,2) = exp(-b[3]*x[i]); in df()
994 fjac(i,3) = -b[2]*x[i]*exp(-b[3]*x[i]); in df()
995 fjac(i,4) = exp(-b[5]*x[i]); in df()
996 fjac(i,5) = -b[4]*x[i]*exp(-b[5]*x[i]); in df()
1075 int df(const VectorXd &b, MatrixXd &fjac) in df()
1078 assert(fjac.rows()==9); in df()
1079 assert(fjac.cols()==3); in df()
1082 fjac(i,0) = 1./(1.+e); in df()
1083 fjac(i,1) = -b[0]*e/(1.+e)/(1.+e); in df()
1084 fjac(i,2) = +b[0]*e*x[i]/(1.+e)/(1.+e); in df()
1152 int df(const VectorXd &b, MatrixXd &fjac) in df()
1155 assert(fjac.rows()==16); in df()
1156 assert(fjac.cols()==3); in df()
1160 fjac(i,0) = e; in df()
1161 fjac(i,1) = b[0]*factor*e; in df()
1162 fjac(i,2) = -b[1]*b[0]*factor*factor*e; in df()
1231 int df(const VectorXd &b, MatrixXd &fjac) in df()
1234 assert(fjac.rows()==6); in df()
1235 assert(fjac.cols()==2); in df()
1238 fjac(i,0) = 1.-e; in df()
1239 fjac(i,1) = b[0]*x[i]*e; in df()
1310 int df(const VectorXd &b, MatrixXd &fjac) in df()
1313 assert(fjac.rows()==33); in df()
1314 assert(fjac.cols()==5); in df()
1316 fjac(i,0) = 1.; in df()
1317 fjac(i,1) = exp(-b[3]*x[i]); in df()
1318 fjac(i,2) = exp(-b[4]*x[i]); in df()
1319 fjac(i,3) = -x[i]*b[1]*exp(-b[3]*x[i]); in df()
1320 fjac(i,4) = -x[i]*b[2]*exp(-b[4]*x[i]); in df()
1403 int df(const VectorXd &b, MatrixXd &fjac) in df()
1406 assert(fjac.rows()==11); in df()
1407 assert(fjac.cols()==4); in df()
1411 fjac(i,0) = (xx+x*b[1]) * factor; in df()
1412 fjac(i,1) = b[0]*x* factor; in df()
1413 fjac(i,2) = - b[0]*(xx+x*b[1]) * x * factor * factor; in df()
1414 fjac(i,3) = - b[0]*(xx+x*b[1]) * factor * factor; in df()
1488 int df(const VectorXd &b, MatrixXd &fjac) in df()
1491 assert(fjac.rows()==154); in df()
1492 assert(fjac.cols()==3); in df()
1495 fjac(i,0) = e; in df()
1496 fjac(i,1) = - b[0]*e/b[2]/(b[1]+x[i]); in df()
1497 fjac(i,2) = b[0]*e*log(b[1]+x[i])/b[2]/b[2]; in df()
1572 int df(const VectorXd &b, MatrixXd &fjac) in df()
1575 assert(fjac.rows()==37); in df()
1576 assert(fjac.cols()==7); in df()
1580 fjac(i,0) = 1.*fact; in df()
1581 fjac(i,1) = x*fact; in df()
1582 fjac(i,2) = xx*fact; in df()
1583 fjac(i,3) = xxx*fact; in df()
1585 fjac(i,4) = x*fact; in df()
1586 fjac(i,5) = xx*fact; in df()
1587 fjac(i,6) = xxx*fact; in df()
1668 int df(const VectorXd &b, MatrixXd &fjac) in df()
1671 assert(fjac.rows()==15); in df()
1672 assert(fjac.cols()==4); in df()
1676 fjac(i,0) = pow(1.+e, power); in df()
1677 fjac(i,1) = power*b[0]*e*pow(1.+e, power-1.); in df()
1678 fjac(i,2) = -power*b[0]*e*x[i]*pow(1.+e, power-1.); in df()
1679 fjac(i,3) = b[0]*power*power*log(1.+e)*pow(1.+e, power); in df()
1756 int df(const VectorXd &b, MatrixXd &fjac) in df()
1759 assert(fjac.rows()==35); in df()
1760 assert(fjac.cols()==3); in df()
1764 fjac(i,0) = e / b[1]; in df()
1765 fjac(i,1) = ((x[i]-b[2])*(x[i]-b[2])/b12-1.) * b[0]*e/b12; in df()
1766 fjac(i,2) = (x[i]-b[2])*e*b[0]/b[1]/b12; in df()