Home
last modified time | relevance | path

Searched refs:kthCol (Results 1 – 24 of 24) sorted by relevance

/dports/science/cantera/cantera-2.5.1-611-gc4d6ecc15/ext/sundials/examples/cvodes/serial/
H A DcvsAdvDiff_ASAi_bnd.c343 realtype *kthCol, hordc, horac, verdc; in Jac() local
364 kthCol = SUNBandMatrix_Column(J,k); in Jac()
368 SM_COLUMN_ELEMENT_B(kthCol,k,k) = -TWO*(verdc+hordc); in Jac()
369 if (i != 1) SM_COLUMN_ELEMENT_B(kthCol,k-MY,k) = hordc + horac; in Jac()
371 if (j != 1) SM_COLUMN_ELEMENT_B(kthCol,k-1,k) = verdc; in Jac()
372 if (j != MY) SM_COLUMN_ELEMENT_B(kthCol,k+1,k) = verdc; in Jac()
437 realtype *kthCol, hordc, horac, verdc; in JacB() local
450 kthCol = SUNBandMatrix_Column(JB,k); in JacB()
454 SM_COLUMN_ELEMENT_B(kthCol,k,k) = TWO*(verdc+hordc); in JacB()
457 if (j != 1) SM_COLUMN_ELEMENT_B(kthCol,k-1,k) = - verdc; in JacB()
[all …]
H A DcvsAdvDiff_bndL.c263 realtype *kthCol, hordc, horac, verdc; in Jac() local
285 kthCol = SUNBandMatrix_Column(J,k); in Jac()
289 SM_COLUMN_ELEMENT_B(kthCol,k,k) = -TWO*(verdc+hordc); in Jac()
290 if (i != 1) SM_COLUMN_ELEMENT_B(kthCol,k-MY,k) = hordc + horac; in Jac()
291 if (i != MX) SM_COLUMN_ELEMENT_B(kthCol,k+MY,k) = hordc - horac; in Jac()
292 if (j != 1) SM_COLUMN_ELEMENT_B(kthCol,k-1,k) = verdc; in Jac()
293 if (j != MY) SM_COLUMN_ELEMENT_B(kthCol,k+1,k) = verdc; in Jac()
H A DcvsAdvDiff_bnd.c263 realtype *kthCol, hordc, horac, verdc; in Jac() local
284 kthCol = SUNBandMatrix_Column(J,k); in Jac()
288 SM_COLUMN_ELEMENT_B(kthCol,k,k) = -TWO*(verdc+hordc); in Jac()
289 if (i != 1) SM_COLUMN_ELEMENT_B(kthCol,k-MY,k) = hordc + horac; in Jac()
290 if (i != MX) SM_COLUMN_ELEMENT_B(kthCol,k+MY,k) = hordc - horac; in Jac()
291 if (j != 1) SM_COLUMN_ELEMENT_B(kthCol,k-1,k) = verdc; in Jac()
292 if (j != MY) SM_COLUMN_ELEMENT_B(kthCol,k+1,k) = verdc; in Jac()
H A DcvsDirectDemo_ls.c608 realtype *kthCol; in Jac2() local
629 kthCol = SM_COLUMN_B(J,k); in Jac2()
630 SM_COLUMN_ELEMENT_B(kthCol,k,k) = -TWO; in Jac2()
631 if (i != P2_MESHX-1) SM_COLUMN_ELEMENT_B(kthCol,k+1,k) = P2_ALPH1; in Jac2()
632 if (j != P2_MESHY-1) SM_COLUMN_ELEMENT_B(kthCol,k+P2_MESHX,k) = P2_ALPH2; in Jac2()
/dports/math/sundials/sundials-5.7.0/examples/cvodes/serial/
H A DcvsAdvDiff_ASAi_bnd.c343 realtype *kthCol, hordc, horac, verdc; in Jac() local
364 kthCol = SUNBandMatrix_Column(J,k); in Jac()
368 SM_COLUMN_ELEMENT_B(kthCol,k,k) = -TWO*(verdc+hordc); in Jac()
369 if (i != 1) SM_COLUMN_ELEMENT_B(kthCol,k-MY,k) = hordc + horac; in Jac()
371 if (j != 1) SM_COLUMN_ELEMENT_B(kthCol,k-1,k) = verdc; in Jac()
372 if (j != MY) SM_COLUMN_ELEMENT_B(kthCol,k+1,k) = verdc; in Jac()
437 realtype *kthCol, hordc, horac, verdc; in JacB() local
450 kthCol = SUNBandMatrix_Column(JB,k); in JacB()
454 SM_COLUMN_ELEMENT_B(kthCol,k,k) = TWO*(verdc+hordc); in JacB()
457 if (j != 1) SM_COLUMN_ELEMENT_B(kthCol,k-1,k) = - verdc; in JacB()
[all …]
H A DcvsAdvDiff_bndL.c263 realtype *kthCol, hordc, horac, verdc; in Jac() local
285 kthCol = SUNBandMatrix_Column(J,k); in Jac()
289 SM_COLUMN_ELEMENT_B(kthCol,k,k) = -TWO*(verdc+hordc); in Jac()
290 if (i != 1) SM_COLUMN_ELEMENT_B(kthCol,k-MY,k) = hordc + horac; in Jac()
291 if (i != MX) SM_COLUMN_ELEMENT_B(kthCol,k+MY,k) = hordc - horac; in Jac()
292 if (j != 1) SM_COLUMN_ELEMENT_B(kthCol,k-1,k) = verdc; in Jac()
293 if (j != MY) SM_COLUMN_ELEMENT_B(kthCol,k+1,k) = verdc; in Jac()
H A DcvsAdvDiff_bnd.c263 realtype *kthCol, hordc, horac, verdc; in Jac() local
284 kthCol = SUNBandMatrix_Column(J,k); in Jac()
288 SM_COLUMN_ELEMENT_B(kthCol,k,k) = -TWO*(verdc+hordc); in Jac()
289 if (i != 1) SM_COLUMN_ELEMENT_B(kthCol,k-MY,k) = hordc + horac; in Jac()
290 if (i != MX) SM_COLUMN_ELEMENT_B(kthCol,k+MY,k) = hordc - horac; in Jac()
291 if (j != 1) SM_COLUMN_ELEMENT_B(kthCol,k-1,k) = verdc; in Jac()
292 if (j != MY) SM_COLUMN_ELEMENT_B(kthCol,k+1,k) = verdc; in Jac()
H A DcvsDirectDemo_ls.c608 realtype *kthCol; in Jac2() local
629 kthCol = SM_COLUMN_B(J,k); in Jac2()
630 SM_COLUMN_ELEMENT_B(kthCol,k,k) = -TWO; in Jac2()
631 if (i != P2_MESHX-1) SM_COLUMN_ELEMENT_B(kthCol,k+1,k) = P2_ALPH1; in Jac2()
632 if (j != P2_MESHY-1) SM_COLUMN_ELEMENT_B(kthCol,k+P2_MESHX,k) = P2_ALPH2; in Jac2()
/dports/science/cantera/cantera-2.5.1-611-gc4d6ecc15/ext/sundials/examples/kinsol/serial/
H A DkinLaplace_picard_bnd.c275 realtype *kthCol; in jac() local
303 kthCol = SUNBandMatrix_Column(J, k); in jac()
304 SM_COLUMN_ELEMENT_B(kthCol,k,k) = -2.0*hdc - 2.0*vdc; in jac()
305 if ( i != (NX-1) ) SM_COLUMN_ELEMENT_B(kthCol,k+1,k) = hdc; in jac()
306 if ( i != 0 ) SM_COLUMN_ELEMENT_B(kthCol,k-1,k) = hdc; in jac()
307 if ( j != (NY-1) ) SM_COLUMN_ELEMENT_B(kthCol,k+NX,k) = vdc; in jac()
308 if ( j != 0 ) SM_COLUMN_ELEMENT_B(kthCol,k-NX,k) = vdc; in jac()
/dports/math/sundials/sundials-5.7.0/examples/kinsol/serial/
H A DkinLaplace_picard_bnd.c285 realtype *kthCol; in jac() local
313 kthCol = SUNBandMatrix_Column(J, k); in jac()
314 SM_COLUMN_ELEMENT_B(kthCol,k,k) = -2.0*hdc - 2.0*vdc; in jac()
315 if ( i != (NX-1) ) SM_COLUMN_ELEMENT_B(kthCol,k+1,k) = hdc; in jac()
316 if ( i != 0 ) SM_COLUMN_ELEMENT_B(kthCol,k-1,k) = hdc; in jac()
317 if ( j != (NY-1) ) SM_COLUMN_ELEMENT_B(kthCol,k+NX,k) = vdc; in jac()
318 if ( j != 0 ) SM_COLUMN_ELEMENT_B(kthCol,k-NX,k) = vdc; in jac()
/dports/science/simbody/simbody-Simbody-3.7/SimTKmath/Integrators/src/CPodes/sundials/tests/
H A DcpsAdvDiff_bnd.c157 realtype *kthCol, hordc, horac, verdc; in Jac() local
178 kthCol = BAND_COL(J,k); in Jac()
182 BAND_COL_ELEM(kthCol,k,k) = -TWO*(verdc+hordc); in Jac()
183 if (i != 1) BAND_COL_ELEM(kthCol,k-MY,k) = hordc + horac; in Jac()
184 if (i != MX) BAND_COL_ELEM(kthCol,k+MY,k) = hordc - horac; in Jac()
185 if (j != 1) BAND_COL_ELEM(kthCol,k-1,k) = verdc; in Jac()
186 if (j != MY) BAND_COL_ELEM(kthCol,k+1,k) = verdc; in Jac()
H A DcpsAdvDiff_bndL.c157 realtype *kthCol, hordc, horac, verdc; in Jac() local
178 kthCol = BAND_COL(J,k); in Jac()
182 BAND_COL_ELEM(kthCol,k,k) = -TWO*(verdc+hordc); in Jac()
183 if (i != 1) BAND_COL_ELEM(kthCol,k-MY,k) = hordc + horac; in Jac()
184 if (i != MX) BAND_COL_ELEM(kthCol,k+MY,k) = hordc - horac; in Jac()
185 if (j != 1) BAND_COL_ELEM(kthCol,k-1,k) = verdc; in Jac()
186 if (j != MY) BAND_COL_ELEM(kthCol,k+1,k) = verdc; in Jac()
H A Dcpsbanx_lap.c157 realtype *kthCol, hordc, horac, verdc; in Jac() local
178 kthCol = BAND_COL(J,k); in Jac()
182 BAND_COL_ELEM(kthCol,k,k) = -TWO*(verdc+hordc); in Jac()
183 if (i != 1) BAND_COL_ELEM(kthCol,k-MY,k) = hordc + horac; in Jac()
184 if (i != MX) BAND_COL_ELEM(kthCol,k+MY,k) = hordc - horac; in Jac()
185 if (j != 1) BAND_COL_ELEM(kthCol,k-1,k) = verdc; in Jac()
186 if (j != MY) BAND_COL_ELEM(kthCol,k+1,k) = verdc; in Jac()
H A Dcpsbanx.c157 realtype *kthCol, hordc, horac, verdc; in Jac() local
178 kthCol = BAND_COL(J,k); in Jac()
182 BAND_COL_ELEM(kthCol,k,k) = -TWO*(verdc+hordc); in Jac()
183 if (i != 1) BAND_COL_ELEM(kthCol,k-MY,k) = hordc + horac; in Jac()
184 if (i != MX) BAND_COL_ELEM(kthCol,k+MY,k) = hordc - horac; in Jac()
185 if (j != 1) BAND_COL_ELEM(kthCol,k-1,k) = verdc; in Jac()
186 if (j != MY) BAND_COL_ELEM(kthCol,k+1,k) = verdc; in Jac()
/dports/science/cantera/cantera-2.5.1-611-gc4d6ecc15/ext/sundials/examples/cvode/serial/
H A DcvAdvDiff_bndL.c263 realtype *kthCol, hordc, horac, verdc; in Jac() local
285 kthCol = SUNBandMatrix_Column(J,k); in Jac()
289 SM_COLUMN_ELEMENT_B(kthCol,k,k) = -TWO*(verdc+hordc); in Jac()
290 if (i != 1) SM_COLUMN_ELEMENT_B(kthCol,k-MY,k) = hordc + horac; in Jac()
291 if (i != MX) SM_COLUMN_ELEMENT_B(kthCol,k+MY,k) = hordc - horac; in Jac()
292 if (j != 1) SM_COLUMN_ELEMENT_B(kthCol,k-1,k) = verdc; in Jac()
293 if (j != MY) SM_COLUMN_ELEMENT_B(kthCol,k+1,k) = verdc; in Jac()
H A DcvAdvDiff_bnd.c264 realtype *kthCol, hordc, horac, verdc; in Jac() local
286 kthCol = SUNBandMatrix_Column(J,k); in Jac()
290 SM_COLUMN_ELEMENT_B(kthCol,k,k) = -TWO*(verdc+hordc); in Jac()
291 if (i != 1) SM_COLUMN_ELEMENT_B(kthCol,k-MY,k) = hordc + horac; in Jac()
292 if (i != MX) SM_COLUMN_ELEMENT_B(kthCol,k+MY,k) = hordc - horac; in Jac()
293 if (j != 1) SM_COLUMN_ELEMENT_B(kthCol,k-1,k) = verdc; in Jac()
294 if (j != MY) SM_COLUMN_ELEMENT_B(kthCol,k+1,k) = verdc; in Jac()
H A DcvDirectDemo_ls.c607 realtype *kthCol; in Jac2() local
628 kthCol = SM_COLUMN_B(J,k); in Jac2()
629 SM_COLUMN_ELEMENT_B(kthCol,k,k) = -TWO; in Jac2()
630 if (i != P2_MESHX-1) SM_COLUMN_ELEMENT_B(kthCol,k+1,k) = P2_ALPH1; in Jac2()
631 if (j != P2_MESHY-1) SM_COLUMN_ELEMENT_B(kthCol,k+P2_MESHX,k) = P2_ALPH2; in Jac2()
/dports/math/sundials/sundials-5.7.0/examples/cvode/serial/
H A DcvAdvDiff_bndL.c263 realtype *kthCol, hordc, horac, verdc; in Jac() local
285 kthCol = SUNBandMatrix_Column(J,k); in Jac()
289 SM_COLUMN_ELEMENT_B(kthCol,k,k) = -TWO*(verdc+hordc); in Jac()
290 if (i != 1) SM_COLUMN_ELEMENT_B(kthCol,k-MY,k) = hordc + horac; in Jac()
291 if (i != MX) SM_COLUMN_ELEMENT_B(kthCol,k+MY,k) = hordc - horac; in Jac()
292 if (j != 1) SM_COLUMN_ELEMENT_B(kthCol,k-1,k) = verdc; in Jac()
293 if (j != MY) SM_COLUMN_ELEMENT_B(kthCol,k+1,k) = verdc; in Jac()
H A DcvAdvDiff_bnd.c264 realtype *kthCol, hordc, horac, verdc; in Jac() local
286 kthCol = SUNBandMatrix_Column(J,k); in Jac()
290 SM_COLUMN_ELEMENT_B(kthCol,k,k) = -TWO*(verdc+hordc); in Jac()
291 if (i != 1) SM_COLUMN_ELEMENT_B(kthCol,k-MY,k) = hordc + horac; in Jac()
292 if (i != MX) SM_COLUMN_ELEMENT_B(kthCol,k+MY,k) = hordc - horac; in Jac()
293 if (j != 1) SM_COLUMN_ELEMENT_B(kthCol,k-1,k) = verdc; in Jac()
294 if (j != MY) SM_COLUMN_ELEMENT_B(kthCol,k+1,k) = verdc; in Jac()
H A DcvDirectDemo_ls.c607 realtype *kthCol; in Jac2() local
628 kthCol = SM_COLUMN_B(J,k); in Jac2()
629 SM_COLUMN_ELEMENT_B(kthCol,k,k) = -TWO; in Jac2()
630 if (i != P2_MESHX-1) SM_COLUMN_ELEMENT_B(kthCol,k+1,k) = P2_ALPH1; in Jac2()
631 if (j != P2_MESHY-1) SM_COLUMN_ELEMENT_B(kthCol,k+P2_MESHX,k) = P2_ALPH2; in Jac2()
/dports/science/cantera/cantera-2.5.1-611-gc4d6ecc15/ext/sundials/examples/cvodes/C_openmp/
H A DcvsAdvDiff_bnd_omp.c295 realtype *kthCol, hordc, horac, verdc; in Jac() local
315 #pragma omp parallel for collapse(2) default(shared) private(i, j, k, kthCol) num_threads(data->nth… in Jac()
319 kthCol = SUNBandMatrix_Column(J,k); in Jac()
323 SM_COLUMN_ELEMENT_B(kthCol,k,k) = -TWO*(verdc+hordc); in Jac()
324 if (i != 1) SM_COLUMN_ELEMENT_B(kthCol,k-MY,k) = hordc + horac; in Jac()
325 if (i != MX) SM_COLUMN_ELEMENT_B(kthCol,k+MY,k) = hordc - horac; in Jac()
326 if (j != 1) SM_COLUMN_ELEMENT_B(kthCol,k-1,k) = verdc; in Jac()
327 if (j != MY) SM_COLUMN_ELEMENT_B(kthCol,k+1,k) = verdc; in Jac()
/dports/science/cantera/cantera-2.5.1-611-gc4d6ecc15/ext/sundials/examples/cvode/C_openmp/
H A DcvAdvDiff_bnd_omp.c296 realtype *kthCol, hordc, horac, verdc; in Jac() local
316 #pragma omp parallel for collapse(2) default(shared) private(i, j, k, kthCol) num_threads(data->nth… in Jac()
320 kthCol = SUNBandMatrix_Column(J,k); in Jac()
324 SM_COLUMN_ELEMENT_B(kthCol,k,k) = -TWO*(verdc+hordc); in Jac()
325 if (i != 1) SM_COLUMN_ELEMENT_B(kthCol,k-MY,k) = hordc + horac; in Jac()
326 if (i != MX) SM_COLUMN_ELEMENT_B(kthCol,k+MY,k) = hordc - horac; in Jac()
327 if (j != 1) SM_COLUMN_ELEMENT_B(kthCol,k-1,k) = verdc; in Jac()
328 if (j != MY) SM_COLUMN_ELEMENT_B(kthCol,k+1,k) = verdc; in Jac()
/dports/math/sundials/sundials-5.7.0/examples/cvode/C_openmp/
H A DcvAdvDiff_bnd_omp.c296 realtype *kthCol, hordc, horac, verdc; in Jac() local
316 #pragma omp parallel for collapse(2) default(shared) private(i, j, k, kthCol) num_threads(data->nth… in Jac()
320 kthCol = SUNBandMatrix_Column(J,k); in Jac()
324 SM_COLUMN_ELEMENT_B(kthCol,k,k) = -TWO*(verdc+hordc); in Jac()
325 if (i != 1) SM_COLUMN_ELEMENT_B(kthCol,k-MY,k) = hordc + horac; in Jac()
326 if (i != MX) SM_COLUMN_ELEMENT_B(kthCol,k+MY,k) = hordc - horac; in Jac()
327 if (j != 1) SM_COLUMN_ELEMENT_B(kthCol,k-1,k) = verdc; in Jac()
328 if (j != MY) SM_COLUMN_ELEMENT_B(kthCol,k+1,k) = verdc; in Jac()
/dports/math/sundials/sundials-5.7.0/examples/cvodes/C_openmp/
H A DcvsAdvDiff_bnd_omp.c295 realtype *kthCol, hordc, horac, verdc; in Jac() local
315 #pragma omp parallel for collapse(2) default(shared) private(i, j, k, kthCol) num_threads(data->nth… in Jac()
319 kthCol = SUNBandMatrix_Column(J,k); in Jac()
323 SM_COLUMN_ELEMENT_B(kthCol,k,k) = -TWO*(verdc+hordc); in Jac()
324 if (i != 1) SM_COLUMN_ELEMENT_B(kthCol,k-MY,k) = hordc + horac; in Jac()
325 if (i != MX) SM_COLUMN_ELEMENT_B(kthCol,k+MY,k) = hordc - horac; in Jac()
326 if (j != 1) SM_COLUMN_ELEMENT_B(kthCol,k-1,k) = verdc; in Jac()
327 if (j != MY) SM_COLUMN_ELEMENT_B(kthCol,k+1,k) = verdc; in Jac()