Home
last modified time | relevance | path

Searched refs:nz (Results 1 – 25 of 197) sorted by relevance

12345678

/petsc/src/mat/tests/output/ !
H A Dex134_1.out7 [0] Local rows 4 nz 32 nz alloced 32 bs 2 mem.
8 [0] on-diagonal part: nz 16
9 [0] off-diagonal part: nz 16
10 [1] Local rows 4 nz 32 nz alloced 32 bs 2 mem.
11 [1] on-diagonal part: nz 16
12 [1] off-diagonal part: nz 16
47 [0] Local rows 4 nz 28 nz alloced 32 bs 2 mem.
48 [0] on-diagonal part: nz 12
49 [0] off-diagonal part: nz 16
50 [1] Local rows 4 nz 12 nz alloced 32 bs 2 mem.
[all …]
H A Dex135_1.out6 [0] Local rows 250 nz 500 nz alloced 500 bs 1 mem 18900.
7 [0] on-diagonal part: nz 250
8 [0] off-diagonal part: nz 250
9 [1] Local rows 250 nz 500 nz alloced 500 bs 1 mem 18900.
10 [1] on-diagonal part: nz 250
11 [1] off-diagonal part: nz 250
12 [2] Local rows 250 nz 500 nz alloced 500 bs 1 mem 18900.
13 [2] on-diagonal part: nz 250
14 [2] off-diagonal part: nz 250
15 [3] Local rows 250 nz 499 nz alloced 500 bs 1 mem 19398.
[all …]
/petsc/src/binding/petsc4py/demo/legacy/perftest/ !
H A DApp.f905 integer(kind=C_INT), intent(in) :: nx, ny, nz local
7 real(kind=C_DOUBLE), intent(in) :: x(nx,ny,nz), xdot(nx,ny,nz)
8 real(kind=C_DOUBLE), intent(inout) :: f(nx,ny,nz)
9 call formfunction_f(nx, ny, nz, h, t, x, xdot, f)
16 integer(kind=C_INT), intent(in) :: nx, ny, nz local
18 real(kind=C_DOUBLE), intent(inout) :: x(nx,ny,nz)
19 call forminitial_f(nx, ny, nz, h, t, x)
45 integer, intent(in) :: nx, ny, nz local
47 real(kind=8), intent(inout) :: x(nx,ny,nz)
54 integer, intent(in) :: nx, ny, nz local
[all …]
H A DApp.pyf6 subroutine formFunction(nx, ny, nz, h, t, x, xdot, f)
10 integer, intent(in), intent(hide) :: nz = shape(x,2)
12 real(kind=8), intent(in) :: x(nx,ny,nz), xdot(nx,ny,nz)
13 real(kind=8), intent(inout) :: f(nx,ny,nz)
16 subroutine formInitial(nx, ny, nz, h, t, x)
20 integer, intent(in), intent(hide) :: nz = shape(x,2)
22 real(kind=8), intent(inout) :: x(nx,ny,nz)
H A Ddriver.py19 def __init__(self, nx, ny, nz): argument
20 self.n = np.array([nx, ny, nz], dtype='i')
40 nx, ny, nz = self.n
43 U = x.getArray(readonly=1).reshape(nx,ny,nz, order='f')
46 Z = U[:,:,nz//2]
56 X, Y = mgrid[0:1:1j*nx,0:1:1j*nz]
67 X, Y = mgrid[0:1:1j*ny,0:1:1j*nz]
79 def run_test(nx,ny,nz,samples,plot=False): argument
90 x = PETSc.Vec().createSeq(nx*ny*nz)
92 app = AppCtx(nx, ny, nz)
[all …]
H A Ddriver.c11 PetscInt nx,ny,nz; member
23 formInitial(&app->nx,&app->ny,&app->nz,app->h,&t,x); in FormInitial()
41 formFunction(&app->nx,&app->ny,&app->nz,app->h,&t,x,xdot,f); in FormFunction()
49 PetscErrorCode RunTest(int nx, int ny, int nz, int loops, double *wt) in RunTest() argument
59 app->nz = nz; app->h[2] = 1./(nz-1); in RunTest()
62 PetscCall(VecSetSizes(x,nx*ny*nz,nx*ny*nz)); in RunTest()
125 int nx=n+1, ny=n+1, nz=n+1; in main() local
126 PetscCall(RunTest(nx,ny,nz,samples,&wt)); in main()
127 …Printf(PETSC_COMM_SELF,"Grid %3d x %3d x %3d -> %f seconds (%2d samples)\n",nx,ny,nz,wt,samples)); in main()
/petsc/src/mat/impls/sbaij/seq/ !
H A Drelax.h64 PetscInt mbs = a->mbs, i, j, nz; in MatMult_SeqSBAIJ_1_ushort() local
83 nz = ai[i + 1] - ai[i]; /* length of i_th row of A */ in MatMult_SeqSBAIJ_1_ushort()
84 if (!nz) continue; /* Move to the next row if the current row is empty */ in MatMult_SeqSBAIJ_1_ushort()
93 …PetscPrefetchBlock(ib + nz, nz, 0, PETSC_PREFETCH_HINT_NTA); /* Indices for the next row (assumes … in MatMult_SeqSBAIJ_1_ushort()
94 PetscPrefetchBlock(v + nz, nz, 0, PETSC_PREFETCH_HINT_NTA); /* Entries for the next row */ in MatMult_SeqSBAIJ_1_ushort()
96 for (j = jmin; j < nz; j++) { in MatMult_SeqSBAIJ_1_ushort()
103 for (j = jmin; j < nz; j++) { in MatMult_SeqSBAIJ_1_ushort()
111 v += nz; in MatMult_SeqSBAIJ_1_ushort()
112 ib += nz; in MatMult_SeqSBAIJ_1_ushort()
117 PetscCall(PetscLogFlops(2.0 * (2.0 * a->nz - nonzerorow) - nonzerorow)); in MatMult_SeqSBAIJ_1_ushort()
[all …]
H A Dsbaijfact2.c15 PetscInt nz, *vj, k, idx, k1; in MatSolve_SeqSBAIJ_N_inplace() local
39 nz = ai[k + 1] - ai[k]; in MatSolve_SeqSBAIJ_N_inplace()
42 while (nz--) { in MatSolve_SeqSBAIJ_N_inplace()
58 nz = ai[k + 1] - ai[k]; in MatSolve_SeqSBAIJ_N_inplace()
61 while (nz--) { in MatSolve_SeqSBAIJ_N_inplace()
76 PetscCall(PetscLogFlops(4.0 * bs2 * a->nz - (bs + 2.0 * bs2) * mbs)); in MatSolve_SeqSBAIJ_N_inplace()
94 PetscInt nz, k; in MatForwardSolve_SeqSBAIJ_N_NaturalOrdering() local
106 nz = ai[k + 1] - ai[k]; in MatForwardSolve_SeqSBAIJ_N_NaturalOrdering()
109 while (nz--) { in MatForwardSolve_SeqSBAIJ_N_NaturalOrdering()
126 PetscInt nz, k; in MatBackwardSolve_SeqSBAIJ_N_NaturalOrdering() local
[all …]
/petsc/src/mat/impls/baij/seq/ !
H A Dbaijsolvnat2.c19 PetscInt jdx, idt, idx, nz, i; in MatSolve_SeqBAIJ_2_NaturalOrdering_inplace() local
32 nz = diag[i] - ai[i]; in MatSolve_SeqBAIJ_2_NaturalOrdering_inplace()
36 while (nz--) { in MatSolve_SeqBAIJ_2_NaturalOrdering_inplace()
51 nz = ai[i + 1] - diag[i] - 1; in MatSolve_SeqBAIJ_2_NaturalOrdering_inplace()
55 while (nz--) { in MatSolve_SeqBAIJ_2_NaturalOrdering_inplace()
70 PetscCall(PetscLogFlops(2.0 * 4 * (a->nz) - 2.0 * A->cmap->n)); in MatSolve_SeqBAIJ_2_NaturalOrdering_inplace()
78 PetscInt i, k, nz, idx, idt, jdx; in MatSolve_SeqBAIJ_2_NaturalOrdering() local
93 nz = ai[i + 1] - ai[i]; in MatSolve_SeqBAIJ_2_NaturalOrdering()
97 PetscPrefetchBlock(vi + nz, nz, 0, PETSC_PREFETCH_HINT_NTA); in MatSolve_SeqBAIJ_2_NaturalOrdering()
98 PetscPrefetchBlock(v + 4 * nz, 4 * nz, 0, PETSC_PREFETCH_HINT_NTA); in MatSolve_SeqBAIJ_2_NaturalOrdering()
[all …]
H A Dbaijsolvnat1.c16 PetscInt jdx, idt, idx, nz, i; in MatSolve_SeqBAIJ_1_NaturalOrdering_inplace() local
28 nz = diag[i] - ai[i]; in MatSolve_SeqBAIJ_1_NaturalOrdering_inplace()
31 while (nz--) { in MatSolve_SeqBAIJ_1_NaturalOrdering_inplace()
43 nz = ai[i + 1] - diag[i] - 1; in MatSolve_SeqBAIJ_1_NaturalOrdering_inplace()
46 while (nz--) { in MatSolve_SeqBAIJ_1_NaturalOrdering_inplace()
57 PetscCall(PetscLogFlops(2.0 * (a->nz) - A->cmap->n)); in MatSolve_SeqBAIJ_1_NaturalOrdering_inplace()
68 PetscInt i, nz; in MatForwardSolve_SeqBAIJ_1_NaturalOrdering() local
81 nz = ai[i + 1] - ai[i]; in MatForwardSolve_SeqBAIJ_1_NaturalOrdering()
83 PetscSparseDenseMinusDot(sum, x, v, vi, nz); in MatForwardSolve_SeqBAIJ_1_NaturalOrdering()
84 v += nz; in MatForwardSolve_SeqBAIJ_1_NaturalOrdering()
[all …]
H A Dbaijsolvtrannat1.c8 PetscInt nz; in MatSolveTranspose_SeqBAIJ_1_NaturalOrdering() local
25 nz = adiag[i] - adiag[i + 1] - 1; in MatSolveTranspose_SeqBAIJ_1_NaturalOrdering()
27 s1 *= v[nz]; /* multiply by inverse of diagonal entry */ in MatSolveTranspose_SeqBAIJ_1_NaturalOrdering()
28 for (j = 0; j < nz; j++) tmp[vi[j]] -= s1 * v[j]; in MatSolveTranspose_SeqBAIJ_1_NaturalOrdering()
36 nz = ai[i + 1] - ai[i]; in MatSolveTranspose_SeqBAIJ_1_NaturalOrdering()
38 for (j = 0; j < nz; j++) tmp[vi[j]] -= s1 * v[j]; in MatSolveTranspose_SeqBAIJ_1_NaturalOrdering()
47 PetscCall(PetscLogFlops(2.0 * a->nz - A->cmap->n)); in MatSolveTranspose_SeqBAIJ_1_NaturalOrdering()
54 PetscInt i, nz; in MatSolveTranspose_SeqBAIJ_1_NaturalOrdering_inplace() local
69 nz = ai[i + 1] - diag[i] - 1; in MatSolveTranspose_SeqBAIJ_1_NaturalOrdering_inplace()
70 while (nz--) x[*vi++] -= (*v++) * s1; in MatSolveTranspose_SeqBAIJ_1_NaturalOrdering_inplace()
[all …]
H A Dbaijfact13.c16 PetscInt *ajtmpold, *ajtmp, nz, row, *ai = a->i, *aj = a->j; in MatILUFactorNumeric_SeqBAIJ_3_inplace() local
37 nz = bi[i + 1] - bi[i]; in MatILUFactorNumeric_SeqBAIJ_3_inplace()
39 for (j = 0; j < nz; j++) { in MatILUFactorNumeric_SeqBAIJ_3_inplace()
45 nz = ai[idx + 1] - ai[idx]; in MatILUFactorNumeric_SeqBAIJ_3_inplace()
48 for (j = 0; j < nz; j++) { in MatILUFactorNumeric_SeqBAIJ_3_inplace()
96 nz = bi[row + 1] - diag_offset[row] - 1; in MatILUFactorNumeric_SeqBAIJ_3_inplace()
98 for (j = 0; j < nz; j++) { in MatILUFactorNumeric_SeqBAIJ_3_inplace()
122 PetscCall(PetscLogFlops(54.0 * nz + 36.0)); in MatILUFactorNumeric_SeqBAIJ_3_inplace()
129 nz = bi[i + 1] - bi[i]; in MatILUFactorNumeric_SeqBAIJ_3_inplace()
130 for (j = 0; j < nz; j++) { in MatILUFactorNumeric_SeqBAIJ_3_inplace()
[all …]
H A Dbaijsolvtran1.c10 PetscInt nz; in MatSolveTranspose_SeqBAIJ_1() local
32 nz = adiag[i] - adiag[i + 1] - 1; in MatSolveTranspose_SeqBAIJ_1()
34 s1 *= v[nz]; /* multiply by inverse of diagonal entry */ in MatSolveTranspose_SeqBAIJ_1()
35 for (j = 0; j < nz; j++) tmp[vi[j]] -= s1 * v[j]; in MatSolveTranspose_SeqBAIJ_1()
43 nz = ai[i + 1] - ai[i]; in MatSolveTranspose_SeqBAIJ_1()
45 for (j = 0; j < nz; j++) tmp[vi[j]] -= s1 * v[j]; in MatSolveTranspose_SeqBAIJ_1()
56 PetscCall(PetscLogFlops(2.0 * a->nz - A->cmap->n)); in MatSolveTranspose_SeqBAIJ_1()
66 PetscInt i, nz; in MatSolveTranspose_SeqBAIJ_1_inplace() local
90 nz = ai[i + 1] - diag[i] - 1; in MatSolveTranspose_SeqBAIJ_1_inplace()
91 while (nz--) t[*vi++] -= (*v++) * s1; in MatSolveTranspose_SeqBAIJ_1_inplace()
[all …]
H A Dbaijfact.c32 PetscInt nzL, nz = bi[i + 1] - bi[i]; in MatLUFactorNumeric_SeqBAIJ_2() local
34 for (PetscInt j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_2()
37 nz = bdiag[i] - bdiag[i + 1]; in MatLUFactorNumeric_SeqBAIJ_2()
39 for (PetscInt j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_2()
42 nz = ai[r[i] + 1] - ai[r[i]]; in MatLUFactorNumeric_SeqBAIJ_2()
45 …for (PetscInt j = 0; j < nz; j++) PetscCall(PetscArraycpy(rtmp + bs2 * ic[ajtmp[j]], v + bs2 * j, … in MatLUFactorNumeric_SeqBAIJ_2()
68 nz = bdiag[row] - bdiag[row + 1] - 1; /* num of entries inU(row,:), excluding diag */ in MatLUFactorNumeric_SeqBAIJ_2()
69 for (PetscInt j = 0; j < nz; j++) { in MatLUFactorNumeric_SeqBAIJ_2()
76 PetscCall(PetscLogFlops(16.0 * nz + 12)); /* flops = 2*bs^3*nz + 2*bs^3 - bs2) */ in MatLUFactorNumeric_SeqBAIJ_2()
84 nz = bi[i + 1] - bi[i]; in MatLUFactorNumeric_SeqBAIJ_2()
[all …]
H A Dbaijsolvtrannat2.c6 PetscInt i, nz, idx, idt, oidx; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering_inplace() local
27 nz = ai[i + 1] - diag[i] - 1; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering_inplace()
28 while (nz--) { in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering_inplace()
42 nz = diag[i] - ai[i]; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering_inplace()
46 while (nz--) { in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering_inplace()
54 PetscCall(PetscLogFlops(2.0 * 4.0 * (a->nz) - 2.0 * A->cmap->n)); in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering_inplace()
62 PetscInt nz, idx, idt, j, i, oidx; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering() local
83 nz = diag[i] - diag[i + 1] - 1; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering()
84 for (j = 0; j > -nz; j--) { in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering()
98 nz = ai[i + 1] - ai[i]; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering()
[all …]
/petsc/src/mat/tests/ !
H A Dmmio.c20 int M, N, nz; in mm_read_unsymmetric_sparse() local
41 if (mm_read_mtx_crd_size(f, &M, &N, &nz) != 0) { in mm_read_unsymmetric_sparse()
48 *nz_ = nz; in mm_read_unsymmetric_sparse()
52 ia = (int *)malloc(nz * sizeof(int)); in mm_read_unsymmetric_sparse()
53 ja = (int *)malloc(nz * sizeof(int)); in mm_read_unsymmetric_sparse()
54 val = (double *)malloc(nz * sizeof(double)); in mm_read_unsymmetric_sparse()
64 for (i = 0; i < nz; i++) { in mm_read_unsymmetric_sparse()
140 int mm_write_mtx_crd_size(FILE *f, int M, int N, int nz) in mm_write_mtx_crd_size() argument
142 if (fprintf(f, "%d %d %d\n", M, N, nz) < 0) return MM_COULD_NOT_WRITE_FILE; in mm_write_mtx_crd_size()
146 int mm_read_mtx_crd_size(FILE *f, int *M, int *N, int *nz) in mm_read_mtx_crd_size() argument
[all …]
H A Dmmloader.c9 PetscInt i, j, nz, *rownz; in MatCreateFromMTX() local
25 …PetscCheck(mm_read_mtx_crd_size(file, &M, &N, &nz) == 0, PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Si… in MatCreateFromMTX()
28 PetscCall(PetscMalloc4(nz, &ia, nz, &ja, nz, &val, M, &rownz)); in MatCreateFromMTX()
34 for (i = 0; i < nz; i++) { in MatCreateFromMTX()
68 …for (j = 0; j < nz; j++) PetscCall(MatSetValues(*A, 1, &ia[j], 1, &ja[j], &val[j], INSERT_VALUES)); in MatCreateFromMTX()
73 …for (j = 0; j < nz; j++) PetscCall(MatSetValues(*A, 1, &ja[j], 1, &ia[j], &val[j], INSERT_VALUES)); in MatCreateFromMTX()
76 for (j = 0; j < nz; j++) { in MatCreateFromMTX()
/petsc/src/mat/impls/aij/seq/ !
H A Dfdaij.c26 PetscInt nz, brows, bcols; in MatFDColoringCreate_SeqXAIJ() local
29 nz = spA->nz; in MatFDColoringCreate_SeqXAIJ()
32 nz = spA->nz; in MatFDColoringCreate_SeqXAIJ()
36 mem = nz * (sizeof(PetscScalar) + sizeof(PetscInt)) + 3 * m * sizeof(PetscInt); in MatFDColoringCreate_SeqXAIJ()
62 PetscErrorCode MatFDColoringSetUpBlocked_AIJ_Private(Mat mat, MatFDColoring c, PetscInt nz) in MatFDColoringSetUpBlocked_AIJ_Private() argument
80 PetscCall(PetscMalloc1(nz, &Jentry_new)); in MatFDColoringSetUpBlocked_AIJ_Private()
99 nz = color_start[j]; in MatFDColoringSetUpBlocked_AIJ_Private()
101 if (Jentry[nz].row >= row_end) { in MatFDColoringSetUpBlocked_AIJ_Private()
102 color_start[j] = nz; in MatFDColoringSetUpBlocked_AIJ_Private()
105 Jentry_new[nz_new].row = Jentry[nz].row + j * mbs; /* index in dy-array */ in MatFDColoringSetUpBlocked_AIJ_Private()
[all …]
H A Dij.c33 PetscInt *work, *ia, *ja, *j, i, nz, row, col; in MatToSymmetricIJ_SeqAIJ() local
44 nz = ai[row + 1] - ai[row]; in MatToSymmetricIJ_SeqAIJ()
46 while (nz--) { in MatToSymmetricIJ_SeqAIJ()
66 nz = ia[m] + (!shiftin); in MatToSymmetricIJ_SeqAIJ()
67 PetscCall(PetscMalloc1(nz, &ja)); in MatToSymmetricIJ_SeqAIJ()
72 nz = ai[row + 1] - ai[row]; in MatToSymmetricIJ_SeqAIJ()
74 while (nz--) { in MatToSymmetricIJ_SeqAIJ()
H A Daijfact.c169 b->maxnz = b->nz = bdiag[0] + 1; in MatLUFactorSymbolic_SeqAIJ()
225 PetscInt i, j, k, nz, nzL, row, *pj; in MatLUFactorNumeric_SeqAIJ() local
250 nz = ai[i + 1] - ai[i]; in MatLUFactorNumeric_SeqAIJ()
251 for (j = 0; j < nz; j++) rs += PetscAbsScalar(v[j]); in MatLUFactorNumeric_SeqAIJ()
270 nz = bi[i + 1] - bi[i]; in MatLUFactorNumeric_SeqAIJ()
272 for (j = 0; j < nz; j++) rtmp[bjtmp[j]] = 0.0; in MatLUFactorNumeric_SeqAIJ()
275 nz = bdiag[i] - bdiag[i + 1]; in MatLUFactorNumeric_SeqAIJ()
277 for (j = 0; j < nz; j++) rtmp[bjtmp[j]] = 0.0; in MatLUFactorNumeric_SeqAIJ()
280 nz = ai[r[i] + 1] - ai[r[i]]; in MatLUFactorNumeric_SeqAIJ()
283 for (j = 0; j < nz; j++) rtmp[ics[ajtmp[j]]] = v[j]; in MatLUFactorNumeric_SeqAIJ()
[all …]
/petsc/src/mat/impls/aij/mpi/ !
H A Dfdmpiaij.c34 PetscInt k, cstart, cend, l, row, col, nz, spidx, i, j; in MatFDColoringApply_BAIJ() local
95 nz = 0; in MatFDColoringApply_BAIJ()
145 row = bs * Jentry2[nz].row; /* local row index */ in MatFDColoringApply_BAIJ()
146 valaddr = Jentry2[nz++].valaddr; in MatFDColoringApply_BAIJ()
158 row = bs * Jentry[nz].row; /* local row index */ in MatFDColoringApply_BAIJ()
159 col = bs * Jentry[nz].col; /* local column index */ in MatFDColoringApply_BAIJ()
160 valaddr = Jentry[nz++].valaddr; in MatFDColoringApply_BAIJ()
185 PetscInt k, cstart, cend, l, row, col, nz; in MatFDColoringApply_AIJ() local
245 nz = 0; in MatFDColoringApply_AIJ()
302 row = Jentry2[nz].row; /* local row index */ in MatFDColoringApply_AIJ()
[all …]
/petsc/src/ts/tests/ !
H A Dex3.c27 PetscInt nz; /* total number of grid points */ member
45 PetscInt i, m, nz, steps, max_steps, k, nphase = 1; in main() local
63 nz = num_z; in main()
64 m = nz - 2; in main()
65 appctx.nz = nz; in main()
69 appctx.max_probsz = nz; in main()
95 PetscCall(PetscMalloc1(nz + 1, &z)); in main()
96 for (i = 0; i < nz; i++) z[i] = (i) * ((zFinal - zInitial) / (nz - 1)); in main()
98 PetscCall(femA(&appctx, nz, z)); in main()
113 for (i = 0; i < nz - 2; i++) { in main()
[all …]
/petsc/src/ksp/ksp/tutorials/ !
H A Dex83f.F9015 PetscCount nz
36 nz = 3*n - 4
41 allocate (ia(n + 1), ja(nz), a(nz))
42 allocate (rows(nz), cols(nz))
70 ja(nz) = n - 1
71 rows(nz) = n - 1; cols(nz) = n - 1
72 a(nz) = 1.0
75 …PetscCallA(MatCreateSeqAIJFromTriple(PETSC_COMM_SELF, n, n, rows, cols, a, Jt, nz, PETSC_FALSE, ie…
82 PetscCallA(MatSetPreallocationCOO(Jr, nz, rows, cols, ierr))
H A Dex88f.F9012 PetscCount nz
29 nz = 3*n - 4
33 allocate (rows(nz), cols(nz), a(nz))
51 rows(nz) = n - 1; cols(nz) = n - 1
52 a(nz) = 1.0
57 PetscCallA(MatSetPreallocationCOO(J, nz, rows, cols, ierr))
H A Dex61f.F9079 integer :: i, j, ij, ij2, ii, jj, nz, ip, dx, dy, icase
128 nz = 0
142 nz = nz + 1
143 ilist(nz) = ij
144 jlist(nz) = ij2
150 alist(nz) = aij
157 print *, 'nz = ', nz
162 ilist(1:nz) = ilist(1:nz) - 1
163 jlist(1:nz) = jlist(1:nz) - 1
182 do ip = 1, nz
[all …]

12345678