| /petsc/src/mat/impls/sbaij/mpi/ |
| H A D | mmsbaij.c | 10 PetscInt i, j, *aj = B->j, ec = 0, *garray, *sgarray; in MatSetUpMultiply_MPISBAIJ() local
34 if (!data) PetscCall(PetscHMapISet(gid1_lid1, gid1, ++ec)); in MatSetUpMultiply_MPISBAIJ()
38 PetscCall(PetscMalloc1(ec, &garray)); in MatSetUpMultiply_MPISBAIJ()
48 PetscCall(PetscSortInt(ec, garray)); in MatSetUpMultiply_MPISBAIJ()
50 for (i = 0; i < ec; i++) PetscCall(PetscHMapISet(gid1_lid1, garray[i] + 1, i + 1)); in MatSetUpMultiply_MPISBAIJ()
61 PetscCall(PetscMalloc2(2 * ec, &sgarray, ec, &ec_owner)); in MatSetUpMultiply_MPISBAIJ()
62 for (i = j = 0; i < ec; i++) { in MatSetUpMultiply_MPISBAIJ()
72 if (!indices[aj[B->i[i] + j]]) ec++; in MatSetUpMultiply_MPISBAIJ()
78 PetscCall(PetscMalloc1(ec, &garray)); in MatSetUpMultiply_MPISBAIJ()
79 PetscCall(PetscMalloc2(2 * ec, &sgarray, ec, &ec_owner)); in MatSetUpMultiply_MPISBAIJ()
[all …]
|
| /petsc/src/mat/impls/baij/mpi/ |
| H A D | mmbaij.c | 11 PetscInt i, j, *aj = B->j, ec = 0, *garray; in MatSetUpMultiply_MPIBAIJ() local
33 PetscCall(PetscHMapISet(gid1_lid1, gid1, ++ec)); in MatSetUpMultiply_MPIBAIJ()
38 PetscCall(PetscMalloc1(ec, &garray)); in MatSetUpMultiply_MPIBAIJ()
48 PetscCall(PetscSortInt(ec, garray)); in MatSetUpMultiply_MPIBAIJ()
50 for (i = 0; i < ec; i++) PetscCall(PetscHMapISet(gid1_lid1, garray[i] + 1, i + 1)); in MatSetUpMultiply_MPIBAIJ()
60 B->nbs = ec; in MatSetUpMultiply_MPIBAIJ()
62 …etscLayoutCreateFromSizes(PetscObjectComm((PetscObject)baij->B), ec * mat->rmap->bs, ec * mat->rma… in MatSetUpMultiply_MPIBAIJ()
70 if (!indices[aj[B->i[i] + j]]) ec++; in MatSetUpMultiply_MPIBAIJ()
76 PetscCall(PetscMalloc1(ec, &garray)); in MatSetUpMultiply_MPIBAIJ()
77 ec = 0; in MatSetUpMultiply_MPIBAIJ()
[all …]
|
| /petsc/src/mat/impls/aij/mpi/ |
| H A D | mmaij.c | 13 PetscInt ec = 0; /* Number of nonzero external columns */ in MatSetUpMultiply_MPIAIJ() local
36 PetscCall(PetscHMapISet(gid1_lid1, gid1, ++ec)); in MatSetUpMultiply_MPIAIJ()
41 PetscCall(PetscMalloc1(ec, &garray)); in MatSetUpMultiply_MPIAIJ()
51 PetscCall(PetscSortInt(ec, garray)); /* sort, and rebuild */ in MatSetUpMultiply_MPIAIJ()
53 for (i = 0; i < ec; i++) PetscCall(PetscHMapISet(gid1_lid1, garray[i] + 1, i + 1)); in MatSetUpMultiply_MPIAIJ()
64 …PetscCall(PetscLayoutCreateFromSizes(PetscObjectComm((PetscObject)aij->B), ec, ec, 1, &aij->B->cma… in MatSetUpMultiply_MPIAIJ()
72 if (!indices[aj[B->i[i] + j]]) ec++; in MatSetUpMultiply_MPIAIJ()
78 PetscCall(PetscMalloc1(ec, &garray)); in MatSetUpMultiply_MPIAIJ()
79 ec = 0; in MatSetUpMultiply_MPIAIJ()
81 if (indices[i]) garray[ec++] = i; in MatSetUpMultiply_MPIAIJ()
[all …]
|
| /petsc/src/mat/impls/sell/mpi/ |
| H A D | mmsell.c | 77 PetscInt i, j, *bcolidx = B->colidx, ec = 0, *garray, totalslices; in MatSetUpMultiply_MPISELL() local
104 if (!data) PetscCall(PetscHMapISet(gid1_lid1, gid1, ++ec)); in MatSetUpMultiply_MPISELL()
110 PetscCall(PetscMalloc1(ec, &garray)); in MatSetUpMultiply_MPISELL()
120 PetscCall(PetscSortInt(ec, garray)); /* sort, and rebuild */ in MatSetUpMultiply_MPISELL()
122 for (i = 0; i < ec; i++) PetscCall(PetscHMapISet(gid1_lid1, garray[i] + 1, i + 1)); in MatSetUpMultiply_MPISELL()
137 …PetscCall(PetscLayoutCreateFromSizes(PetscObjectComm((PetscObject)sell->B), ec, ec, 1, &sell->B->c… in MatSetUpMultiply_MPISELL()
147 if (!indices[bcolidx[j]]) ec++; in MatSetUpMultiply_MPISELL()
154 PetscCall(PetscMalloc1(ec, &garray)); in MatSetUpMultiply_MPISELL()
155 ec = 0; in MatSetUpMultiply_MPISELL()
157 if (indices[i]) garray[ec++] = i; in MatSetUpMultiply_MPISELL()
[all …]
|
| /petsc/src/snes/tutorials/ |
| H A D | ex16.c | 416 …***x, CoordField ***c, PetscInt i, PetscInt j, PetscInt k, Field *ex, CoordField *ec, AppCtx *user) in GatherElementData() argument
431 for (m = 0; m < 3; m++) ec[idx][m] = c[k + kk][j + jj][i + ii][m]; in GatherElementData()
437 void QuadraturePointGeometricJacobian(CoordField *ec, PetscInt qi, PetscInt qj, PetscInt qk, PetscS… in QuadraturePointGeometricJacobian() argument
447 J[0] += grad[3 * bidx + 0] * ec[idx][0]; in QuadraturePointGeometricJacobian()
448 J[1] += grad[3 * bidx + 1] * ec[idx][0]; in QuadraturePointGeometricJacobian()
449 J[2] += grad[3 * bidx + 2] * ec[idx][0]; in QuadraturePointGeometricJacobian()
450 J[3] += grad[3 * bidx + 0] * ec[idx][1]; in QuadraturePointGeometricJacobian()
451 J[4] += grad[3 * bidx + 1] * ec[idx][1]; in QuadraturePointGeometricJacobian()
452 J[5] += grad[3 * bidx + 2] * ec[idx][1]; in QuadraturePointGeometricJacobian()
453 J[6] += grad[3 * bidx + 0] * ec[idx][2]; in QuadraturePointGeometricJacobian()
[all …]
|
| /petsc/lib/petsc/bin/maint/ |
| H A D | xclude | 99 petsc-dist/include/ec.h
117 petsc-dist/src/ec
|
| /petsc/src/mat/impls/aij/seq/ |
| H A D | aij.c | 3567 PetscInt gid, lid, i, ec, nz = aij->nz; in MatSeqAIJCompactOutExtraColumns_SeqAIJ() local
3575 ec = 0; in MatSeqAIJCompactOutExtraColumns_SeqAIJ()
3581 PetscCall(PetscHMapISet(gid1_lid1, gid1, ++ec)); in MatSeqAIJCompactOutExtraColumns_SeqAIJ()
3585 PetscCall(PetscMalloc1(ec, &garray)); in MatSeqAIJCompactOutExtraColumns_SeqAIJ()
3595 PetscCall(PetscSortInt(ec, garray)); /* sort, and rebuild */ in MatSeqAIJCompactOutExtraColumns_SeqAIJ()
3597 for (i = 0; i < ec; i++) PetscCall(PetscHMapISet(gid1_lid1, garray[i] + 1, i + 1)); in MatSeqAIJCompactOutExtraColumns_SeqAIJ()
3607 PetscCall(PetscLayoutCreateFromSizes(PetscObjectComm((PetscObject)mat), ec, ec, 1, &mat->cmap)); in MatSeqAIJCompactOutExtraColumns_SeqAIJ()
|
| /petsc/share/petsc/datafiles/meshes/ |
| H A D | testcase3D.cas | 6214 (udf/ec-rate-fcn "none")
|