/dports/math/gsl/gsl-2.7/ode-initval2/ |
H A D | test.c | 976 dfdy[0 * NRINGMOD + 10] = -dfdy[0 * NRINGMOD + 9]; in jac_ringmod() 977 dfdy[0 * NRINGMOD + 13] = dfdy[0 * NRINGMOD + 7]; in jac_ringmod() 978 dfdy[1 * NRINGMOD + 1] = dfdy[0 * NRINGMOD + 0]; in jac_ringmod() 979 dfdy[1 * NRINGMOD + 8] = dfdy[0 * NRINGMOD + 7]; in jac_ringmod() 980 dfdy[1 * NRINGMOD + 11] = dfdy[0 * NRINGMOD + 9]; in jac_ringmod() 981 dfdy[1 * NRINGMOD + 12] = dfdy[0 * NRINGMOD + 10]; in jac_ringmod() 982 dfdy[1 * NRINGMOD + 14] = dfdy[0 * NRINGMOD + 13]; in jac_ringmod() 1009 dfdy[8 * NRINGMOD + 1] = dfdy[7 * NRINGMOD + 0]; in jac_ringmod() 1016 dfdy[11 * NRINGMOD + 1] = dfdy[9 * NRINGMOD + 0]; in jac_ringmod() 1017 dfdy[11 * NRINGMOD + 4] = dfdy[9 * NRINGMOD + 2]; in jac_ringmod() [all …]
|
H A D | rk2imp.c | 55 gsl_matrix *dfdy; /* Jacobian matrix */ member 166 state->dfdy = gsl_matrix_alloc (dim, dim); in rk2imp_alloc() 168 if (state->dfdy == 0) in rk2imp_alloc() 186 gsl_matrix_free (state->dfdy); in rk2imp_alloc() 204 gsl_matrix_free (state->dfdy); in rk2imp_alloc() 240 gsl_matrix *const dfdy = state->dfdy; in rk2imp_apply() local 304 int s = GSL_ODEIV_JA_EVAL (sys, t, y, dfdy->data, dfdt); in rk2imp_apply() 319 M = dfdy->size1; in rk2imp_apply() 320 N = dfdy->size2; in rk2imp_apply() 326 double aij = gsl_matrix_get (dfdy, i, j); in rk2imp_apply() [all …]
|
H A D | rk1imp.c | 50 gsl_matrix *dfdy; /* Jacobian matrix */ member 161 state->dfdy = gsl_matrix_alloc (dim, dim); in rk1imp_alloc() 163 if (state->dfdy == 0) in rk1imp_alloc() 181 gsl_matrix_free (state->dfdy); in rk1imp_alloc() 199 gsl_matrix_free (state->dfdy); in rk1imp_alloc() 235 gsl_matrix *const dfdy = state->dfdy; in rk1imp_apply() local 283 int s = GSL_ODEIV_JA_EVAL (sys, t, y, dfdy->data, dfdt); in rk1imp_apply() 294 int s = modnewton1_init ((void *) esol, A, h, dfdy, sys); in rk1imp_apply() 331 int s = modnewton1_init ((void *) esol, A, h / 2.0, dfdy, sys); in rk1imp_apply() 475 gsl_matrix_free (state->dfdy); in rk1imp_free()
|
H A D | rk4imp.c | 56 gsl_matrix *dfdy; /* Jacobian matrix */ member 167 state->dfdy = gsl_matrix_alloc (dim, dim); in rk4imp_alloc() 169 if (state->dfdy == 0) in rk4imp_alloc() 187 gsl_matrix_free (state->dfdy); in rk4imp_alloc() 205 gsl_matrix_free (state->dfdy); in rk4imp_alloc() 241 gsl_matrix *const dfdy = state->dfdy; in rk4imp_apply() local 294 int s = GSL_ODEIV_JA_EVAL (sys, t, y, dfdy->data, dfdt); in rk4imp_apply() 305 int s = modnewton1_init ((void *) esol, A, h, dfdy, sys); in rk4imp_apply() 350 int s = modnewton1_init ((void *) esol, A, h / 2.0, dfdy, sys); in rk4imp_apply() 508 gsl_matrix_free (state->dfdy); in rk4imp_free()
|
H A D | bsimp.c | 70 gsl_matrix *dfdy; member 204 const gsl_matrix * dfdy, in bsimp_step_local() argument 243 gsl_matrix_set (a_mat, i, j, -h * gsl_matrix_get (dfdy, i, j)); in bsimp_step_local() 372 state->dfdy = gsl_matrix_alloc (dim, dim); in bsimp_alloc() 413 gsl_matrix *dfdy = state->dfdy; in bsimp_apply() local 446 int s = GSL_ODEIV_JA_EVAL (sys, t_local, y, dfdy->data, dfdt); in bsimp_apply() 467 dfdt, dfdy, in bsimp_apply() 540 gsl_matrix_free (state->dfdy); in bsimp_free()
|
/dports/math/gsl/gsl-2.7/ode-initval/ |
H A D | test.c | 54 dfdy[0] = 0.0; in jac_linear() 81 dfdy[0] = y[0]; in jac_exp() 111 dfdy[0] = 0.0; in jac_sin() 112 dfdy[1] = -1.0; in jac_sin() 113 dfdy[2] = 1.0; in jac_sin() 114 dfdy[3] = 0.0; in jac_sin() 180 dfdy[0] = 0.0; in jac_xsin() 182 dfdy[2] = 1.0; in jac_xsin() 183 dfdy[3] = 0.0; in jac_xsin() 313 dfdy[6] = c3; in jac_oregonator() [all …]
|
H A D | rk2simp.c | 46 double *dfdy; /* Jacobian */ member 103 state->dfdy = (double *) malloc (dim * dim * sizeof (double)); in rk2simp_alloc() 105 if (state->dfdy == 0) in rk2simp_alloc() 123 free (state->dfdy); in rk2simp_alloc() 136 free (state->dfdy); in rk2simp_alloc() 150 free (state->dfdy); in rk2simp_alloc() 183 gsl_matrix_view J = gsl_matrix_view_array (state->dfdy, dim, dim); in rk2simp_step() 191 s = GSL_ODEIV_JA_EVAL (sys, t, y0, state->dfdy, state->dfdt); in rk2simp_step() 365 free (state->dfdy); in rk2simp_free()
|
H A D | bsimp.c | 69 gsl_matrix *dfdy; member 203 const gsl_matrix * dfdy, in bsimp_step_local() argument 243 gsl_matrix_set (a_mat, i, j, -h * gsl_matrix_get (dfdy, i, j)); in bsimp_step_local() 369 state->dfdy = gsl_matrix_alloc (dim, dim); in bsimp_alloc() 411 gsl_matrix *dfdy = state->dfdy; in bsimp_apply() local 444 int s = GSL_ODEIV_JA_EVAL (sys, t_local, y, dfdy->data, dfdt); in bsimp_apply() 465 dfdt, dfdy, in bsimp_apply() 537 gsl_matrix_free (state->dfdy); in bsimp_free()
|
/dports/science/simbody/simbody-Simbody-3.7/SimTKmath/src/ |
H A D | Differentiator.cpp | 268 dfdy.resize(1,1); in calcJacobian() 334 dfdy.resize(1,getNumParameters()); in calcJacobian() 472 Real dfdy; in calcDerivative() local 476 return dfdy; in calcDerivative() 542 rep->frep.calcJacobian(*rep,m,y0,0,dfdy); in calcJacobian() 545 return dfdy; in calcJacobian() 718 dfdy = (fyplus-fy0)/h; in calcDerivative() 722 dfdy = (fyplus-fyminus)/(2*h); in calcDerivative() 771 dfdy.resize(NFunctions,NParameters); in calcJacobian() 782 dfdy(i) = (fyptmp-fy0)/h; in calcJacobian() [all …]
|
/dports/cad/gmsh/gmsh-4.9.2-source/Solver/ |
H A D | FuncGradDisc.h | 74 void gradient(double x, double y, double z, double &dfdx, double &dfdy, in gradient() argument 96 dfdy = 0; in gradient() 111 dfdy = dfdy + std::abs((*_ls)(v->x(), v->y(), v->z())) * dNdy; in gradient() 112 dfdy = dfdy - (*_ls)(v->x(), v->y(), v->z()) * dNdy; in gradient() 129 dfdy = dfdy + std::abs((*_ls)(v->x(), v->y(), v->z())) * dNdy; in gradient() 130 dfdy = dfdy + (*_ls)(v->x(), v->y(), v->z()) * dNdy; in gradient()
|
H A D | FuncHeaviside.h | 37 double &dfdy, double &dfdz) const in gradient() argument 39 dfdx = dfdy = dfdz = 0.0; in gradient() 42 double &dfdy, double &dfdz, MElement *e) const in gradient() argument 44 dfdx = dfdy = dfdz = 0.0; in gradient()
|
/dports/science/qwalk/mainline-1.0.1-300-g1b7e381/src/basis_function/ |
H A D | Cubic_spline_calc.cpp | 921 dfdy= -6.*x*y*z; in calcLap() 1001 dfdy= -6*x*y*z; in calcLap() 1021 dfdy= 3*(v1 - v2)*z; in calcLap() 1060 dfdy= (v1 - 3*v2)*x; in calcLap() 1149 *(ptr++)=dfdy; in calcHessian() 1162 *(ptr++)=x*dfdy; in calcHessian() 1167 *(ptr++)=dfdy+x*dfdxy; in calcHessian() 1195 *(ptr++)=dfdy+z*dfdyz; in calcHessian() 1202 *(ptr++)=v*dfdy; in calcHessian() 1230 *(ptr++)=v*dfdy; in calcHessian() [all …]
|
/dports/math/py-pygsl/pygsl-2.3.0/examples/ |
H A D | odeiv.py | 23 dfdy = numx.ones((2,2),) * 1. 24 dfdy[0, 0] = 0.0 25 dfdy[0, 1] = 1.0 26 dfdy[1, 0] = -2.0 * mu * y[0] * y[1] - 1.0 27 dfdy[1, 1] = -mu * (y[0]**2 - 1.0) 29 return dfdy, dfdt
|
H A D | odeiv2.py | 71 dfdy = numx.zeros((l,l), numx.float_) 72 dfdy[0,0] = 0.0 73 dfdy[0,1] = 1.0 74 dfdy[1,0] = - 2.0 * mu * y0 * y1 - 1.0 75 dfdy[1,1] = - mu * (y0**2 - 1.0) 81 return dfdy, dfdt
|
/dports/science/nwchem/nwchem-7b21660b82ebd85ef659f6fba7e1e73433b0bd0a/src/nwxc/nwad/unit_tests/ |
H A D | test_dx.F | 18 double precision dfdx, dfdy 26 call submaxima_dx(tmp(1),tmp(2),dfdx,dfdy) 39 if (abs((dfdy-ady)/(dfdy+ady)).gt.tol) then 40 write(*,*)"DY : fail:",i,tmp(1),tmp(2),ady,dfdy
|
H A D | test_dxy.F | 20 double precision dfdx, dfdy, d2fdx2, d2fdxy, d2fdy2 29 call submaxima_dxy(tmp(1),tmp(2),dfdx,dfdy,d2fdx2,d2fdxy,d2fdy2) 48 if (abs((dfdy-ady)/(dfdy+ady)).gt.tol) then 49 write(*,*)"DY : fail:",i,tmp(1),tmp(2),ady,dfdy
|
H A D | submaxima_dx.max | 6 print(" subroutine submaxima_dx(x,y,dfdx,dfdy)"), 8 print(" double precision x,y,dfdx,dfdy"), 10 Fortran(dfdy = df1y(x,y)),
|
H A D | submaxima_dxy.max | 9 print(" subroutine submaxima_dxy(x,y,dfdx,dfdy,"), 12 print(" double precision x,y,dfdx,dfdy"), 15 Fortran(dfdy = df1y(x,y)),
|
/dports/science/nwchem-data/nwchem-7.0.2-release/src/nwxc/nwad/unit_tests/ |
H A D | test_dx.F | 18 double precision dfdx, dfdy 26 call submaxima_dx(tmp(1),tmp(2),dfdx,dfdy) 39 if (abs((dfdy-ady)/(dfdy+ady)).gt.tol) then 40 write(*,*)"DY : fail:",i,tmp(1),tmp(2),ady,dfdy
|
H A D | test_dxy.F | 20 double precision dfdx, dfdy, d2fdx2, d2fdxy, d2fdy2 29 call submaxima_dxy(tmp(1),tmp(2),dfdx,dfdy,d2fdx2,d2fdxy,d2fdy2) 48 if (abs((dfdy-ady)/(dfdy+ady)).gt.tol) then 49 write(*,*)"DY : fail:",i,tmp(1),tmp(2),ady,dfdy
|
H A D | submaxima_dx.max | 6 print(" subroutine submaxima_dx(x,y,dfdx,dfdy)"), 8 print(" double precision x,y,dfdx,dfdy"), 10 Fortran(dfdy = df1y(x,y)),
|
/dports/audio/praat/praat-6.2.03/external/gsl/ |
H A D | gsl_ode-initval__rk2simp.c | 46 double *dfdy; /* Jacobian */ member 103 state->dfdy = (double *) malloc (dim * dim * sizeof (double)); in rk2simp_alloc() 105 if (state->dfdy == 0) in rk2simp_alloc() 123 free (state->dfdy); in rk2simp_alloc() 136 free (state->dfdy); in rk2simp_alloc() 150 free (state->dfdy); in rk2simp_alloc() 182 gsl_matrix_view J = gsl_matrix_view_array (state->dfdy, dim, dim); in rk2simp_step() 190 s = GSL_ODEIV_JA_EVAL (sys, t, y0, state->dfdy, state->dfdt); in rk2simp_step() 362 free (state->dfdy); in rk2simp_free()
|
H A D | gsl_ode-initval__bsimp.c | 69 gsl_matrix *dfdy; member 203 const gsl_matrix * dfdy, in bsimp_step_local() argument 243 gsl_matrix_set (a_mat, i, j, -h * gsl_matrix_get (dfdy, i, j)); in bsimp_step_local() 369 state->dfdy = gsl_matrix_alloc (dim, dim); in bsimp_alloc() 411 gsl_matrix *dfdy = state->dfdy; in bsimp_apply() local 444 int s = GSL_ODEIV_JA_EVAL (sys, t_local, y, dfdy->data, dfdt); in bsimp_apply() 465 dfdt, dfdy, in bsimp_apply() 537 gsl_matrix_free (state->dfdy); in bsimp_free()
|
/dports/math/py-pygsl/pygsl-2.3.0/tests/ |
H A D | odeiv_test.py | 21 dfdy = numx.ones((2,2), Float) 22 dfdy[0, 0] = 0.0 23 dfdy[0, 1] = 1.0 24 dfdy[1, 0] = -2.0 * mu * y[0] * y[1] - 1.0 25 dfdy[1, 1] = -mu * (y[0]**2 - 1.0) 31 return dfdy, dfdt
|
/dports/science/py-scipy/scipy-1.7.1/benchmarks/benchmarks/ |
H A D | test_functions.py | 139 dfdy = 4. * (x + 2. * y - 7.) + 2. * (2. * x + y - 5.) 140 return np.array([dfdx, dfdy]) 161 dfdy = (2. * (1.5 - x + x * y) * (x) + 164 return np.array([dfdx, dfdy]) 241 dfdy = 2. * deriv1 * x[1] - term2 * pi * sin(2. * pi * x[1]) 243 return np.array([dfdx, dfdy]) 264 dfdy = ((x[0] - 1.)**2 * 2. * 3. * pi * cos(3. * pi * x[1]) * sin(3. * 269 return np.array([dfdx, dfdy])
|