| /petsc/src/ts/tests/output/ |
| H A D | ex2_be_adapt.out | 1 At t = 0.000000e+00 u = 1.000000e+00 1.000000e+00 1.000000e+00 2 At t = 0.000000e+00 errors = 0.000000e+00 0.000000e+00 0.000000e+00 3 At t = 1.000000e-03 u = 1.003010e+00 1.004014e+00 1.003010e+00 4 At t = 1.000000e-03 errors = 5.028445e-06 7.040158e-06 5.028445e-06 5 At t = 2.000000e-03 u = 1.006030e+00 1.008042e+00 1.006030e+00 6 At t = 2.000000e-03 errors = 1.009117e-05 1.412871e-05 1.009117e-05 7 At t = 4.000000e-03 u = 1.012111e+00 1.016155e+00 1.012111e+00 8 At t = 4.000000e-03 errors = 3.052675e-05 4.274366e-05 3.052675e-05 9 At t = 8.000000e-03 u = 1.024437e+00 1.032612e+00 1.024437e+00 10 At t = 8.000000e-03 errors = 1.139107e-04 1.595211e-04 1.139107e-04 [all …]
|
| H A D | ex2_euler.out | 1 At t = 0.000000e+00 u = 1.000000e+00 1.000000e+00 1.000000e+00 2 At t = 0.000000e+00 errors = 0.000000e+00 0.000000e+00 0.000000e+00 3 At t = 1.000000e-03 u = 1.003000e+00 1.004000e+00 1.003000e+00 4 At t = 1.000000e-03 errors = -5.005672e-06 -7.008007e-06 -5.005672e-06 5 At t = 2.000000e-03 u = 1.006010e+00 1.008014e+00 1.006010e+00 6 At t = 2.000000e-03 errors = -1.004541e-05 -1.406411e-05 -1.004541e-05 7 At t = 3.000000e-03 u = 1.009030e+00 1.012042e+00 1.009030e+00 8 At t = 3.000000e-03 errors = -1.511939e-05 -2.116855e-05 -1.511939e-05 9 At t = 4.000000e-03 u = 1.012060e+00 1.016084e+00 1.012060e+00 10 At t = 4.000000e-03 errors = -2.022779e-05 -2.832159e-05 -2.022779e-05 [all …]
|
| H A D | ex2_beuler.out | 1 At t = 0.000000e+00 u = 1.000000e+00 1.000000e+00 1.000000e+00 2 At t = 0.000000e+00 errors = 0.000000e+00 0.000000e+00 0.000000e+00 3 At t = 1.000000e-03 u = 1.003010e+00 1.004014e+00 1.003010e+00 4 At t = 1.000000e-03 errors = 5.028445e-06 7.040158e-06 5.028445e-06 5 At t = 2.000000e-03 u = 1.006030e+00 1.008042e+00 1.006030e+00 6 At t = 2.000000e-03 errors = 1.009117e-05 1.412871e-05 1.009117e-05 7 At t = 3.000000e-03 u = 1.009060e+00 1.012084e+00 1.009060e+00 8 At t = 3.000000e-03 errors = 1.518836e-05 2.126592e-05 1.518836e-05 9 At t = 4.000000e-03 u = 1.012101e+00 1.016141e+00 1.012101e+00 10 At t = 4.000000e-03 errors = 2.032017e-05 2.845202e-05 2.032017e-05 [all …]
|
| H A D | ex2_rk.out | 1 At t = 0.000000e+00 u = 1.000000e+00 1.000000e+00 1.000000e+00 2 At t = 0.000000e+00 errors = 0.000000e+00 0.000000e+00 0.000000e+00 4 At t = 1.000000e-03 u = 1.003005e+00 1.004007e+00 1.003005e+00 5 At t = 1.000000e-03 errors = -4.836798e-12 -6.837864e-12 -4.836798e-12 7 At t = 1.100000e-02 u = 1.033613e+00 1.044858e+00 1.033613e+00 8 At t = 1.100000e-02 errors = -4.883664e-08 -6.904508e-08 -4.883664e-08 10 At t = 6.696691e-02 u = 1.225076e+00 1.301734e+00 1.225076e+00 11 At t = 6.696691e-02 errors = -5.123787e-05 -7.244106e-05 -5.123787e-05 13 At t = 1.216158e-01 u = 1.450026e+00 1.605847e+00 1.450026e+00 14 At t = 1.216158e-01 errors = -1.180127e-04 -1.668553e-04 -1.180127e-04 [all …]
|
| H A D | ex4_8.out | 1 At t[0] = 0.00e+00 u= 1.00e+00 at the center 2 At t[10] = 1.00e-01 u= 4.72e-01 at the center 3 At t[20] = 2.00e-01 u= 2.24e-01 at the center 4 At t[30] = 3.00e-01 u= 1.09e-01 at the center 5 At t[40] = 4.00e-01 u= 5.37e-02 at the center 6 At t[50] = 5.00e-01 u= 2.70e-02 at the center 7 At t[60] = 6.00e-01 u= 1.38e-02 at the center 8 At t[70] = 7.00e-01 u= 7.15e-03 at the center 9 At t[80] = 8.00e-01 u= 3.74e-03 at the center 10 At t[90] = 9.00e-01 u= 1.98e-03 at the center [all …]
|
| H A D | ex2_sundials.out | 1 At t = 0.000000e+00 u = 1.000000e+00 1.000000e+00 1.000000e+00 2 At t = 0.000000e+00 errors = 0.000000e+00 0.000000e+00 0.000000e+00 3 At t = 1.000000e+00 u = 2.620596e+01 3.631673e+01 2.620596e+01 4 At t = 1.000000e+00 errors = 7.994325e-04 1.130509e-03 7.994325e-04
|
| H A D | ex4_3.out | 1 At t[0] = 0.00e+00 u= 1.00e+00 at the center 2 At t[10] = 1.00e+00 u= 8.42e-04 at the center
|
| H A D | ex4.out | 1 At t[0] = 0.00e+00 u= 1.00e+00 at the center 2 At t[10] = 1.00e+00 u= 4.67e-03 at the center
|
| H A D | ex4_7.out | 1 At t[0] = 0.00e+00 u= 1.00e+00 at the center 12 At t[10] = 1.00e+00 u= 4.67e-03 at the center
|
| /petsc/src/ts/tutorials/hamiltonian/output/ |
| H A D | ex1_2.out | 1 At time 0.000000, Energy = 81.92, Modified Energy = 81.92 2 At time 0.010000, Energy = 81.9203, Modified Energy = 82.1714 3 At time 0.020000, Energy = 81.9208, Modified Energy = 82.0648 4 At time 0.030000, Energy = 81.9207, Modified Energy = 81.7522 5 At time 0.040000, Energy = 81.9203, Modified Energy = 81.6796 6 At time 0.050000, Energy = 81.92, Modified Energy = 81.9506 7 At time 0.060000, Energy = 81.9203, Modified Energy = 82.1785 8 At time 0.070000, Energy = 81.9208, Modified Energy = 82.0383 9 At time 0.080000, Energy = 81.9207, Modified Energy = 81.7299 10 At time 0.090000, Energy = 81.9202, Modified Energy = 81.6933 [all …]
|
| H A D | ex1_1.out | 1 At time 0.000000, Energy = 81.92, Modified Energy = 81.92 2 At time 0.010000, Energy = 81.6695, Modified Energy = 81.92 3 At time 0.020000, Energy = 81.7775, Modified Energy = 81.92 4 At time 0.030000, Energy = 82.09, Modified Energy = 81.92 5 At time 0.040000, Energy = 82.1612, Modified Energy = 81.92 6 At time 0.050000, Energy = 81.8895, Modified Energy = 81.92 7 At time 0.060000, Energy = 81.6625, Modified Energy = 81.92 8 At time 0.070000, Energy = 81.8041, Modified Energy = 81.92 9 At time 0.080000, Energy = 82.1122, Modified Energy = 81.92 10 At time 0.090000, Energy = 82.1473, Modified Energy = 81.92 [all …]
|
| /petsc/share/petsc/matlab/ |
| H A D | generatePetscTestFiles.m | 17 [At,bt] = PetscBinaryRead('spd-real-int32-float64', 'indices','int32','precision','float64'); 18 if max(max(max(abs(At-A))),max(abs(bt-b))) ~= 0 25 [At,bt] = PetscBinaryRead('spd-real-int64-float64','indices','int64','precision','float64'); 26 if max(max(max(abs(At-A))),max(abs(bt-b))) ~= 0 34 [At,bt] = PetscBinaryRead('spd-real-int32-float32', 'indices','int32','precision','float32'); 35 if max(max(max(abs(At-A))),max(abs(bt-b))) ~= 0 42 [At,bt] = PetscBinaryRead('spd-real-int64-float32','indices','int64','precision','float32'); 43 if max(max(max(abs(At-A))),max(abs(bt-b))) ~= 0 58 [At,bt] = PetscBinaryRead('ns-real-int32-float64','indices','int32','precision','float64'); 59 if max(max(max(abs(At-A))),max(abs(bt-b))) ~= 0 [all …]
|
| /petsc/src/mat/tests/ |
| H A D | ex80.c | 15 Mat A, At; in main() local 108 PetscCall(MatCreate(PETSC_COMM_WORLD, &At)); in main() 109 PetscCall(MatSetSizes(At, 4, 4, 16, 16)); in main() 110 PetscCall(MatSetType(At, MATMPIADJ)); in main() 113 PetscCall(MatSetValues(At, 1, &row, ia[i + 1] - ia[i], ja + ia[i], NULL, INSERT_VALUES)); in main() 115 PetscCall(MatAssemblyBegin(At, MAT_FINAL_ASSEMBLY)); in main() 116 PetscCall(MatAssemblyEnd(At, MAT_FINAL_ASSEMBLY)); in main() 117 PetscCall(MatEqual(A, At, &equal)); in main() 119 PetscCall(MatDestroy(&At)); in main()
|
| H A D | ex123.c | 7 Mat A, At, AAt, T = NULL; in main() local 108 PetscCall(MatTranspose(T, MAT_INITIAL_MATRIX, &At)); in main() 110 PetscCall(MatMatMult(T, At, MAT_INITIAL_MATRIX, PETSC_DETERMINE, &AAt)); in main() 113 PetscCall(MatMatMult(At, T, MAT_INITIAL_MATRIX, PETSC_DETERMINE, &AAt)); in main() 117 PetscCall(MatDestroy(&At)); in main() 158 PetscCall(MatTranspose(T, MAT_INITIAL_MATRIX, &At)); in main() 160 PetscCall(MatMatMult(T, At, MAT_INITIAL_MATRIX, PETSC_DETERMINE, &AAt)); in main() 163 PetscCall(MatMatMult(At, T, MAT_INITIAL_MATRIX, PETSC_DETERMINE, &AAt)); in main() 167 PetscCall(MatDestroy(&At)); in main() 224 PetscCall(MatTranspose(T, MAT_INITIAL_MATRIX, &At)); in main() [all …]
|
| H A D | ex222.c | 31 Mat At; in main() local 32 PetscCall(MatTranspose(A, MAT_INITIAL_MATRIX, &At)); in main() 33 PetscCall(MatAXPY(A, 1.0, At, DIFFERENT_NONZERO_PATTERN)); in main() 34 PetscCall(MatDestroy(&At)); in main()
|
| H A D | ex267.c | 120 Mat A, At, Aht; in main() local 153 PetscCall(MatCreateTranspose(A, &At)); in main() 156 PetscCall(TestMatrix("LU T", At, nrhs, inplace, PETSC_FALSE)); in main() 159 PetscCall(TestMatrix("Chol T", At, nrhs, inplace, PETSC_TRUE)); in main() 173 PetscCall(TestMatrix("LU T nonsym", At, nrhs, inplace, PETSC_FALSE)); in main() 177 PetscCall(MatDestroy(&At)); in main()
|
| H A D | ex53.c | 7 Mat A, B, C, At, Bt; in main() local 162 PetscCall(MatTranspose(A, MAT_INITIAL_MATRIX, &At)); in main() 166 PetscCall(MatMult(At, xx, s1)); in main() 173 PetscCall(MatDestroy(&At)); in main()
|
| /petsc/src/mat/impls/aij/seq/ |
| H A D | symtranspose.c | 13 Mat At; in MatTransposeSymbolic_SeqAIJ() local 43 …PetscCall(MatCreateSeqAIJWithArrays(PetscObjectComm((PetscObject)A), an, am, ati, atj, NULL, &At)); in MatTransposeSymbolic_SeqAIJ() 44 PetscCall(MatSetBlockSizes(At, A->cmap->bs, A->rmap->bs)); in MatTransposeSymbolic_SeqAIJ() 45 PetscCall(MatSetType(At, ((PetscObject)A)->type_name)); in MatTransposeSymbolic_SeqAIJ() 46 at = (Mat_SeqAIJ *)At->data; in MatTransposeSymbolic_SeqAIJ() 51 *B = At; in MatTransposeSymbolic_SeqAIJ() 58 Mat At; in MatTranspose_SeqAIJ() local 93 At = *B; in MatTranspose_SeqAIJ() 124 … PetscCall(MatCreateSeqAIJWithArrays(PetscObjectComm((PetscObject)A), an, am, ati, atj, ata, &At)); in MatTranspose_SeqAIJ() 125 PetscCall(MatSetBlockSizes(At, A->cmap->bs, A->rmap->bs)); in MatTranspose_SeqAIJ() [all …]
|
| H A D | matptap.c | 42 atb->At = Pt; in MatProductSymbolic_PtAP_SeqAIJ_SeqAIJ() 286 PetscCall(MatTranspose(P, MAT_REUSE_MATRIX, &atb->At)); in MatPtAPNumeric_SeqAIJ_SeqAIJ() 290 PetscCall((*C->ops->matmatmultnumeric)(atb->At, A, P, C)); in MatPtAPNumeric_SeqAIJ_SeqAIJ()
|
| /petsc/src/ts/impls/arkimex/ |
| H A D | arkimex.c | 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() [all …]
|
| H A D | arkimex.h | 11 PetscReal *At, *bt, *ct; /* Stiff tableau */ member
|
| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex74.c | 88 PetscScalar *A, *B, *At, *b, *zvals, one = 1.0; in main() local 157 PetscCall(PetscMalloc2(nstages, &At, PetscSqr(nstages), &B)); in main() 159 At[s] = 0; in main() 161 At[s] += A[s + nstages * t]; /* Row sums of */ in main() 188 PetscCall(MatCreateKAIJ(Identity, nstages, 1, NULL, At, &R)); in main() 242 PetscCall(PetscFree2(At, B)); in main()
|
| /petsc/src/mat/utils/ |
| H A D | multequal.c | 116 …atMatMultEqual_Private(Mat A, Mat B, Mat C, PetscInt n, PetscBool *flg, PetscBool At, PetscBool Bt) in MatMatMultEqual_Private() argument 136 PetscValidLogicalCollectiveBool(A, At, 6); in MatMatMultEqual_Private() 141 if (At) { in MatMatMultEqual_Private() 154 sop = sops[(At ? 1 : 0) + 2 * (Bt ? 1 : 0)]; in MatMatMultEqual_Private() 163 if (At) { in MatMatMultEqual_Private() 179 if (At) { in MatMatMultEqual_Private()
|
| /petsc/src/ts/tutorials/ |
| H A D | ex19.c | 122 …At[3][3] = {{0, 0, 0}, {0.12132034355964257320, 0.29289321881345247560, 0}, {0.2071067811865475244… in RegisterMyARK2() local 123 …PetscCall(TSARKIMEXRegister("myark2", 2, 3, &At[0][0], NULL, NULL, &A[0][0], NULL, NULL, bembedt, … in RegisterMyARK2()
|
| H A D | ex16.c | 149 …At[3][3] = {{0, 0, 0}, {0.12132034355964257320, 0.29289321881345247560, 0}, {0.2071067811865475244… in RegisterMyARK2() local 150 …PetscCall(TSARKIMEXRegister("myark2", 2, 3, &At[0][0], NULL, NULL, &A[0][0], NULL, NULL, bembedt, … in RegisterMyARK2()
|