Lines Matching refs:At

763     const PetscReal At[3][3] = {  in TSARKIMEXRegisterAll()  local
770 …PetscCall(TSARKIMEXRegister(TSARKIMEX1BEE, 2, 3, &At[0][0], b, NULL, &A[0][0], b, NULL, bembedt, b… in TSARKIMEXRegisterAll()
777 const PetscReal At[2][2] = { in TSARKIMEXRegisterAll() local
784 …PetscCall(TSARKIMEXRegister(TSARKIMEXARS122, 2, 2, &At[0][0], b, NULL, &A[0][0], b, NULL, bembedt,… in TSARKIMEXRegisterAll()
791 const PetscReal At[2][2] = { in TSARKIMEXRegisterAll() local
798 …PetscCall(TSARKIMEXRegister(TSARKIMEXA2, 2, 2, &At[0][0], b, NULL, &A[0][0], b, NULL, bembedt, bem… in TSARKIMEXRegisterAll()
805 const PetscReal At[2][2] = { in TSARKIMEXRegisterAll() local
819 …PetscCall(TSARKIMEXRegister(TSARKIMEXL2, 2, 2, &At[0][0], b, NULL, &A[0][0], b, NULL, bembedt, bem… in TSARKIMEXRegisterAll()
827 const PetscReal At[3][3] = { in TSARKIMEXRegisterAll() local
838 …PetscCall(TSARKIMEXRegister(TSARKIMEX2C, 2, 3, &At[0][0], NULL, NULL, &A[0][0], NULL, NULL, bembed… in TSARKIMEXRegisterAll()
846 const PetscReal At[3][3] = { in TSARKIMEXRegisterAll() local
857 …PetscCall(TSARKIMEXRegister(TSARKIMEX2D, 2, 3, &At[0][0], NULL, NULL, &A[0][0], NULL, NULL, bembed… in TSARKIMEXRegisterAll()
865 const PetscReal At[3][3] = { in TSARKIMEXRegisterAll() local
876 …PetscCall(TSARKIMEXRegister(TSARKIMEX2E, 2, 3, &At[0][0], NULL, NULL, &A[0][0], NULL, NULL, bembed… in TSARKIMEXRegisterAll()
884 const PetscReal At[3][3] = { in TSARKIMEXRegisterAll() local
889 …PetscCall(TSARKIMEXRegister(TSARKIMEXPRSSP2, 2, 3, &At[0][0], NULL, NULL, &A[0][0], NULL, NULL, NU… in TSARKIMEXRegisterAll()
898 const PetscReal At[4][4] = { in TSARKIMEXRegisterAll() local
911 …PetscCall(TSARKIMEXRegister(TSARKIMEX3, 3, 4, &At[0][0], NULL, NULL, &A[0][0], NULL, NULL, bembedt… in TSARKIMEXRegisterAll()
921 const PetscReal At[5][5] = { in TSARKIMEXRegisterAll() local
928 …PetscCall(TSARKIMEXRegister(TSARKIMEXARS443, 3, 5, &At[0][0], NULL, NULL, &A[0][0], NULL, NULL, NU… in TSARKIMEXRegisterAll()
938 const PetscReal At[5][5] = { in TSARKIMEXRegisterAll() local
945 …PetscCall(TSARKIMEXRegister(TSARKIMEXBPR3, 3, 5, &At[0][0], NULL, NULL, &A[0][0], NULL, NULL, NULL… in TSARKIMEXRegisterAll()
956 const PetscReal At[6][6] = { in TSARKIMEXRegisterAll() local
973 …PetscCall(TSARKIMEXRegister(TSARKIMEX4, 4, 6, &At[0][0], NULL, NULL, &A[0][0], NULL, NULL, bembedt… in TSARKIMEXRegisterAll()
986 const PetscReal At[8][8] = { in TSARKIMEXRegisterAll() local
1007 …PetscCall(TSARKIMEXRegister(TSARKIMEX5, 5, 8, &At[0][0], NULL, NULL, &A[0][0], NULL, NULL, bembedt… in TSARKIMEXRegisterAll()
1032 PetscCall(PetscFree6(t->At, t->bt, t->ct, t->A, t->b, t->c)); in TSARKIMEXRegisterDestroy()
1104 PetscErrorCode TSARKIMEXRegister(TSARKIMEXType name, PetscInt order, PetscInt s, const PetscReal At in TSARKIMEXRegister() argument
1124 PetscCall(PetscMalloc6(s * s, &t->At, s, &t->bt, s, &t->ct, s * s, &t->A, s, &t->b, s, &t->c)); in TSARKIMEXRegister()
1125 PetscCall(PetscArraycpy(t->At, At, s * s)); in TSARKIMEXRegister()
1133 for (i = 0; i < s; i++) t->bt[i] = At[(s - 1) * s + i]; in TSARKIMEXRegister()
1144 for (j = 0, t->ct[i] = 0; j < s; j++) t->ct[i] += At[i * s + j]; in TSARKIMEXRegister()
1155 if (t->At[(s - 1) * s + i] != t->bt[i]) t->stiffly_accurate = PETSC_FALSE; in TSARKIMEXRegister()
1159 if (t->At[i] != 0.0) t->explicit_first_stage = PETSC_FALSE; in TSARKIMEXRegister()
1203 PetscErrorCode TSDIRKRegister(TSDIRKType name, PetscInt order, PetscInt s, const PetscReal At[], co… in TSDIRKRegister() argument
1206 …PetscCall(TSARKIMEXRegister(name, order, s, At, bt, ct, NULL, NULL, NULL, bembedt, NULL, pinterp, … in TSDIRKRegister()
1335 const PetscReal *At = tab->At, *A = tab->A, *ct = tab->ct, *c = tab->c; in TSStep_ARKIMEX() local
1439 if (At[i * s + i] == 0) { /* This stage is explicit */ in TSStep_ARKIMEX()
1442 for (j = 0; j < i; j++) w[j] = h * At[i * s + j]; in TSStep_ARKIMEX()
1451 ark->scoeff = 1. / At[i * s + i]; in TSStep_ARKIMEX()
1454 for (j = 0; j < i; j++) w[j] = h * At[i * s + j]; in TSStep_ARKIMEX()
1565 …const PetscReal *At = tab->At, *A = tab->A, *ct = tab->ct, *c = tab->c, *b = tab->b, *bt = tab->bt; in TSAdjointStep_ARKIMEX() local
1583 if (At[i * s + i] == 0) { // This stage is explicit in TSAdjointStep_ARKIMEX()
1586 …ark->scoeff = -1. / At[i * s + i]; // this makes shift=ark->scoeff/ts->time_step positive since ts… in TSAdjointStep_ARKIMEX()
1600 for (j = i + 1; j < s; j++) w[j - i - 1] = -At[j * s + i]; in TSAdjointStep_ARKIMEX()
1639 if (At[i * s + i] == 0) { // This stage is explicit in TSAdjointStep_ARKIMEX()
1646 PetscCall(VecScale(VecsDeltaLam[nadj * s + i], 1. / At[i * s + i])); in TSAdjointStep_ARKIMEX()
1657 …PetscCall(VecAXPY(ts->vecs_sensip[nadj], -At[i * s + i] * adjoint_time_step, VecsSensiPTemp[nadj])… in TSAdjointStep_ARKIMEX()
2138 PetscCall(PetscFormatRealArray(buf, sizeof(buf), "% 8.6f", tab->s, tab->At + i * tab->s)); in TSView_ARKIMEX()