Searched refs:I_s (Results 1 – 2 of 2) sorted by relevance
| /petsc/src/ts/impls/implicit/irk/ |
| H A D | irk.c | 21 PetscScalar *A_inv, *A_inv_rowsum, *I_s; member
67 …scReal *binterp, const PetscScalar *A_inv, const PetscScalar *A_inv_rowsum, const PetscScalar *I_s) in TSIRKTableauCreate() argument
74 …Malloc3(PetscSqr(nstages), &tab->A, PetscSqr(nstages), &tab->A_inv, PetscSqr(nstages), &tab->I_s)); in TSIRKTableauCreate()
83 if (I_s) PetscCall(PetscArraycpy(tab->I_s, I_s, PetscSqr(nstages))); in TSIRKTableauCreate()
92 PetscScalar *gauss_A, *gauss_A_inv, *gauss_A_inv_rowsum, *I_s; in TSIRKCreate_Gauss() local
100 …&gauss_A, PetscSqr(nstages), &gauss_A_inv, nstages, &gauss_A_inv_rowsum, PetscSqr(nstages), &I_s)); in TSIRKCreate_Gauss()
150 I_s[i + nstages * j] = 1. * (i == j); in TSIRKCreate_Gauss()
153 …auCreate(ts, nstages, gauss_A_real, gauss_b, gauss_c, NULL, gauss_A_inv, gauss_A_inv_rowsum, I_s)); in TSIRKCreate_Gauss()
155 PetscCall(PetscFree4(gauss_A, gauss_A_inv, gauss_A_inv_rowsum, I_s)); in TSIRKCreate_Gauss()
405 PetscCall(PetscFree3(tab->A, tab->A_inv, tab->I_s)); in TSIRKTableauReset()
[all …]
|
| /petsc/src/tao/tutorials/ |
| H A D | ex4.c | 51 PetscInt Istart, Iend, i, j, Ii, gridN, I_n, I_s, I_e, I_w; in CreateMatrix() local
76 I_s = i * gridN + j - 1; in CreateMatrix()
77 if (j - 1 < 0) I_s = -1; in CreateMatrix()
84 PetscCall(MatSetValue(ctx->F, Ii, I_s, -1., INSERT_VALUES)); in CreateMatrix()
|