Home
last modified time | relevance | path

Searched refs:dfdy (Results 1 – 25 of 231) sorted by relevance

12345678910

/dports/math/gsl/gsl-2.7/ode-initval2/
H A Dtest.c976 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 Drk2imp.c55 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 Drk1imp.c50 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 Drk4imp.c56 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 Dbsimp.c70 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 Dtest.c54 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 Drk2simp.c46 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 Dbsimp.c69 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 DDifferentiator.cpp268 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 DFuncGradDisc.h74 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 DFuncHeaviside.h37 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 DCubic_spline_calc.cpp921 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 Dodeiv.py23 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 Dodeiv2.py71 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 Dtest_dx.F18 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 Dtest_dxy.F20 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 Dsubmaxima_dx.max6 print(" subroutine submaxima_dx(x,y,dfdx,dfdy)"),
8 print(" double precision x,y,dfdx,dfdy"),
10 Fortran(dfdy = df1y(x,y)),
H A Dsubmaxima_dxy.max9 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 Dtest_dx.F18 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 Dtest_dxy.F20 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 Dsubmaxima_dx.max6 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 Dgsl_ode-initval__rk2simp.c46 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 Dgsl_ode-initval__bsimp.c69 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 Dodeiv_test.py21 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 Dtest_functions.py139 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])

12345678910