| /petsc/src/vec/vec/tests/ |
| H A D | ex39.c | 16 Vec X, Y, F, X1, Y1, X2, Y2, F1, F2; in main() local 40 PetscCall(VecCreate(PETSC_COMM_SELF, &X1)); in main() 41 PetscCall(VecSetSizes(X1, N / 2, N / 2)); in main() 42 PetscCall(VecSetFromOptions(X1)); in main() 43 PetscCall(VecDuplicate(X1, &X2)); in main() 44 PetscCall(VecDuplicate(X1, &Y1)); in main() 45 PetscCall(VecDuplicate(X1, &Y2)); in main() 46 PetscCall(VecDuplicate(X1, &F1)); in main() 47 PetscCall(VecDuplicate(X1, &F2)); in main() 54 PetscCall(VecPlaceArray(X1, x)); in main() [all …]
|
| H A D | ex32.c | 15 PetscInt *X, *X1, *XR, *XSO, *W, *Y, *Z, *XP, *X1P; in main() local 38 PetscCall(PetscCalloc6(n, &X, n, &X1, n, &XR, n, &XSO, n, &Y, n, &Z)); in main() 75 PetscCall(PetscArraycpy(X1, order ? XSO : XR, n)); in main() 79 PetscCall(PetscIntSortSemiOrdered(n, X1)); in main() 89 …X1[i], PETSC_COMM_SELF, PETSC_ERR_PLIB, "PetscIntSortSemiOrdered() rep %" PetscInt_FMT " X1[%" Pet… in main() 91 …for (i = 0; i < n - 1; i++) PetscCheck(X1[i] <= X1[i + 1], PETSC_COMM_SELF, PETSC_ERR_PLIB, "Petsc… in main() 93 PetscCall(PetscArrayzero(X1, n)); in main() 110 PetscCall(PetscArraycpy(X1, order ? XSO : XR, n)); in main() 118 PetscCall(PetscIntSortSemiOrderedWithArray(n, X1, X1P)); in main() 138 …FMT "! randomSeed %" PetscInt64_FMT " orderedSeed %" PetscInt64_FMT, l, i, X1[i], i, X[i], seedr, … in main() [all …]
|
| /petsc/src/mat/tests/ |
| H A D | ex97.c | 51 static PetscErrorCode CheckMatrices(Mat A, Mat B, Vec left, Vec right, Vec X, Vec Y, Vec X1, Vec Y1) in CheckMatrices() argument 81 PetscCall(VecCopy(X1, rtmp[0])); in CheckMatrices() 82 PetscCall(VecCopy(X1, rtmp[1])); in CheckMatrices() 87 PetscCall(MatMultTransposeAdd(A, Y, X1, rtmp[0])); in CheckMatrices() 88 PetscCall(MatMultTransposeAdd(B, Y, X1, rtmp[1])); in CheckMatrices() 100 Vec left, right, X, Y, X1, Y1; in main() local 128 PetscCall(VecDuplicate(right, &X1)); in main() 138 PetscCall(VecSetRandom(X1, NULL)); in main() 145 PetscCall(VecSet(X1, 3.0)); in main() 148 PetscCall(CheckMatrices(Asub, Bsub, left, right, X, Y, X1, Y1)); in main() [all …]
|
| H A D | ex99.c | 51 static PetscErrorCode CheckMatrices(Mat A, Mat B, Vec left, Vec right, Vec X, Vec Y, Vec X1, Vec Y1) in CheckMatrices() argument 83 PetscCall(VecCopy(X1, rtmp[0])); in CheckMatrices() 84 PetscCall(VecCopy(X1, rtmp[1])); in CheckMatrices() 89 PetscCall(MatMultTransposeAdd(A, Y, X1, rtmp[0])); in CheckMatrices() 90 PetscCall(MatMultTransposeAdd(B, Y, X1, rtmp[1])); in CheckMatrices() 102 Vec left, right, X, Y, X1, Y1; in main() local 129 PetscCall(VecDuplicate(right, &X1)); in main() 139 PetscCall(VecSetRandom(X1, NULL)); in main() 146 PetscCall(VecSet(X1, 3.0)); in main() 149 PetscCall(CheckMatrices(Asub, Bsub, left, right, X, Y, X1, Y1)); in main() [all …]
|
| /petsc/src/ts/impls/implicit/alpha/ |
| H A D | alpha1.c | 23 Vec X0, Xa, X1; member 73 Vec X1 = X, V1 = th->V1; in TSAlpha_StageVecs() local 83 PetscCall(VecWAXPY(V1, -1.0, X0, X1)); in TSAlpha_StageVecs() 86 PetscCall(VecWAXPY(Xa, -1.0, X0, X1)); in TSAlpha_StageVecs() 118 Vec X0 = ts->vec_sol, X1, X2 = th->X1; in TSAlpha_Restart() local 122 PetscCall(VecDuplicate(X0, &X1)); in TSAlpha_Restart() 137 PetscCall(VecCopy(th->X0, X1)); in TSAlpha_Restart() 138 PetscCall(TSAlpha_SNESSolve(ts, NULL, X1)); in TSAlpha_Restart() 139 PetscCall(TSPostStage(ts, th->stage_time, 0, &X1)); in TSAlpha_Restart() 140 PetscCall(TSAdaptCheckStage(ts->adapt, ts, th->stage_time, X1, &stageok)); in TSAlpha_Restart() [all …]
|
| H A D | alpha2.c | 24 Vec X0, Xa, X1; member 75 static PetscErrorCode TSAlpha_ApplyPredictor(TS ts, Vec X1) in TSAlpha_ApplyPredictor() argument 81 if (th->predictor) PetscCall(th->predictor(ts, th->X0, th->V0, th->A0, X1, th->predictor_ctx)); in TSAlpha_ApplyPredictor() 82 else PetscCall(VecCopy(th->X0, X1)); in TSAlpha_ApplyPredictor() 107 Vec X1 = X, V1 = th->V1, A1 = th->A1; in TSAlpha_StageVecs() local 118 PetscCall(VecWAXPY(A1, -1.0, X0, X1)); in TSAlpha_StageVecs() 125 PetscCall(VecWAXPY(Xa, -1.0, X0, X1)); in TSAlpha_StageVecs() 160 Vec X0 = ts->vec_sol, X1, X2 = th->X1; in TSAlpha_Restart() local 165 PetscCall(VecDuplicate(X0, &X1)); in TSAlpha_Restart() 182 PetscCall(TSAlpha_ApplyPredictor(ts, X1)); in TSAlpha_Restart() [all …]
|
| /petsc/src/mat/tutorials/ |
| H A D | ex2.c | 14 Vec X, X1, X2, Y, Z, Z1, Z2; in main() local 93 PetscCall(VecCreateSeqWithArray(PETSC_COMM_SELF, 1, size1, x, &X1)); in main() 103 PetscCall(MatMult(A11, X1, Z1)); in main() 106 PetscCall(MatMultAdd(A21, X1, Z2, Z2)); in main() 133 PetscCall(MatMult(A11, X1, Z1)); in main() 136 PetscCall(MatMultAdd(A21, X1, Z2, Z2)); in main() 145 PetscCall(MatMultTranspose(A11, X1, Z1)); in main() 148 PetscCall(MatMultTransposeAdd(A12, X1, Z2, Z2)); in main() 167 PetscCall(VecDestroy(&X1)); in main()
|
| /petsc/src/snes/tutorials/ |
| H A D | ex73f90t.F90 | 137 subroutine InitialGuessLocal(solver, X1, ierr) argument 140 Vec:: X1 156 PetscCall(VecGetOwnershipRange(X1, low, high, ierr)) 167 PetscCall(VecSetValues(X1, ione, [row], [v], INSERT_VALUES, ierr)) 238 subroutine FormJacobianLocal(X1, jac, solver, add_nl_term, ierr) argument 241 Vec:: X1 264 PetscCall(VecGetOwnershipRange(X1, low, high, ierr)) 265 PetscCall(VecGetArrayRead(X1, lx_v, ierr)) 300 PetscCall(VecRestoreArrayRead(X1, lx_v, ierr)) 364 subroutine FormFunctionNLTerm(X1, F1, solver, ierr) argument [all …]
|
| H A D | ex59.c | 16 #define X1 5.0 macro 35 const PetscScalar b = X1; in main() 146 ff[n - 1] = d * d * (xx[n - 1] - X1); in FormFunction()
|
| /petsc/src/ts/tutorials/ |
| H A D | ex43.c | 53 PetscErrorCode Predictor(TS ts, Vec X0, Vec V0, Vec A0, Vec X1, PetscCtx ctx) in Predictor() argument 61 PetscCall(VecCopy(X0, X1)); in Predictor() 62 PetscCall(VecAXPY(X1, dt, V0)); in Predictor() 63 PetscCall(VecAXPY(X1, accel_fac, A0)); in Predictor()
|
| /petsc/src/ksp/ksp/impls/gmres/dgmres/ |
| H A D | dgmres.c | 81 X1 = NULL; in KSPSetUp_DGMRES() 281 PetscCall(PetscFree(X1)); in KSPDestroy_DGMRES() 794 if (!X1) { in KSPDGMRESApplyDeflation_DGMRES() 795 PetscCall(PetscMalloc1(bmax, &X1)); in KSPDGMRESApplyDeflation_DGMRES() 798 PetscCall(VecMDot(x, r, UU, X1)); in KSPDGMRESApplyDeflation_DGMRES() 801 PetscCall(PetscArraycpy(X2, X1, br)); in KSPDGMRESApplyDeflation_DGMRES() 805 PetscCallBLAS("LAPACKgetrs", LAPACKgetrs_("N", &br, &nrhs, TTF, &bmax, INVP, X1, &bmax, &info)); in KSPDGMRESApplyDeflation_DGMRES() 817 …fs", LAPACKgerfs_("N", &br, &nrhs, TT, &bmax, TTF, &bmax, INVP, X2, &bmax, X1, &bmax, &ferr, &berr… in KSPDGMRESApplyDeflation_DGMRES() 821 for (i = 0; i < r; i++) X2[i] = X1[i] / lambda - X2[i]; in KSPDGMRESApplyDeflation_DGMRES()
|
| H A D | dgmresimpl.h | 79 #define X1 dgmres->x1 macro
|
| /petsc/src/tao/constrained/impls/ipm/ |
| H A D | ipm.c | 971 static PetscErrorCode IPMGatherRHS(Tao tao, Vec RHS, Vec X1, Vec X2, Vec X3, Vec X4) in IPMGatherRHS() argument 980 if (X1) { in IPMGatherRHS() 981 PetscCall(VecScatterBegin(ipmP->rhs1, X1, RHS, INSERT_VALUES, SCATTER_FORWARD)); in IPMGatherRHS() 982 PetscCall(VecScatterEnd(ipmP->rhs1, X1, RHS, INSERT_VALUES, SCATTER_FORWARD)); in IPMGatherRHS() 1001 static PetscErrorCode IPMScatterStep(Tao tao, Vec STEP, Vec X1, Vec X2, Vec X3, Vec X4) in IPMScatterStep() argument 1011 if (X1) { in IPMScatterStep() 1012 PetscCall(VecScatterBegin(ipmP->step1, STEP, X1, INSERT_VALUES, SCATTER_FORWARD)); in IPMScatterStep() 1013 PetscCall(VecScatterEnd(ipmP->step1, STEP, X1, INSERT_VALUES, SCATTER_FORWARD)); in IPMScatterStep()
|
| /petsc/doc/changes/ |
| H A D | 216.md | 6 - Added support for Cray X1
|
| /petsc/src/sys/classes/draw/impls/null/ |
| H A D | drawnull.c | 60 static PetscErrorCode PetscDrawTriangle_Null(PetscDraw draw, PetscReal X1, PetscReal Y_1, PetscReal… in PetscDrawTriangle_Null() argument
|
| /petsc/doc/miscellaneous/ |
| H A D | acknowledgements.md | 25 - Richard Mills, the `MATAIJPERM` matrix format for the Cray X1; universal F90 array
|
| /petsc/src/sys/classes/draw/impls/x/ |
| H A D | xops.c | 167 static PetscErrorCode PetscDrawTriangle_X(PetscDraw draw, PetscReal X1, PetscReal Y_1, PetscReal X2… in PetscDrawTriangle_X() argument 175 pt[0].x = (short)XTRANS(draw, XiWin, X1); in PetscDrawTriangle_X() 184 x1 = XTRANS(draw, XiWin, X1); in PetscDrawTriangle_X()
|
| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex42.c | 1783 Vec f, X, X1; in solve_stokes_3d_coupled() local 2066 PetscCall(VecDuplicate(X, &X1)); in solve_stokes_3d_coupled() 2074 PetscCall(KSPSolve(ksp_S, f, X1)); in solve_stokes_3d_coupled() 2076 PetscCall(VecDestroy(&X1)); in solve_stokes_3d_coupled()
|
| /petsc/src/benchmarks/results/ |
| H A D | performance_cfd_2_10 | 450 ex10 on a X1 named phoenix with 1 processor, by kaushik Thu May 29 16:06:33 2003 [400/800 MHz 12.8G…
|
| /petsc/include/ |
| H A D | petscts.h | 1577 …EF typedef PetscErrorCode TSAlpha2PredictorFn(TS ts, Vec X0, Vec V0, Vec A0, Vec X1, PetscCtx ctx);
|