| /petsc/src/ksp/ksp/utils/lmvm/ |
| H A D | lmbasis.c | 13 LMBasis basis; in LMBasisCreate() local
24 PetscCall(PetscNew(&basis)); in LMBasisCreate()
25 *basis_p = basis; in LMBasisCreate()
26 basis->m = m; in LMBasisCreate()
27 basis->k = 0; in LMBasisCreate()
28 basis->vecs = backing; in LMBasisCreate()
32 static PetscErrorCode LMBasisGetVec_Internal(LMBasis basis, PetscInt idx, PetscMemoryAccessMode mod… in LMBasisGetVec_Internal() argument
35 PetscAssertPointer(basis, 1); in LMBasisGetVec_Internal()
37 PetscValidLogicalCollectiveInt(basis->vecs, idx, 2); in LMBasisGetVec_Internal()
38 …basis->k, PetscObjectComm((PetscObject)basis->vecs), PETSC_ERR_ARG_OUTOFRANGE, "Asked for index %"… in LMBasisGetVec_Internal()
[all …]
|
| H A D | lmvmutils.c | 1031 LMBasis basis; in MatLMVMGetUpdatedBasis() local
1039 *basis_p = lmvm->basis[type]; in MatLMVMGetUpdatedBasis()
1050 *basis_p = lmvm->basis[type == LMBASIS_B0S ? LMBASIS_S : LMBASIS_Y]; in MatLMVMGetUpdatedBasis()
1054 LMBasis orig_basis = (type == LMBASIS_B0S) ? lmvm->basis[LMBASIS_S] : lmvm->basis[LMBASIS_Y]; in MatLMVMGetUpdatedBasis()
1058 …if (!lmvm->basis[type]) PetscCall(LMBasisCreate(MatLMVMBasisSizeOf(type) == LMBASIS_S ? lmvm->Xpre… in MatLMVMGetUpdatedBasis()
1059 basis = lmvm->basis[type]; in MatLMVMGetUpdatedBasis()
1060 …PetscCall(MatLMVMUpdateOpVecs(B, orig_basis, basis, (type == LMBASIS_B0S) ? MatLMVMApplyJ0Fwd : Ma… in MatLMVMGetUpdatedBasis()
1061 *basis_p = basis; in MatLMVMGetUpdatedBasis()
1071 …if (!lmvm->basis[type]) PetscCall(LMBasisCreate(MatLMVMBasisSizeOf(type) == LMBASIS_S ? lmvm->Xpre… in MatLMVMGetUpdatedBasis()
1072 basis = lmvm->basis[type]; in MatLMVMGetUpdatedBasis()
[all …]
|
| H A D | lmvmimpl.c | 62 for (PetscInt i = 0; i < LMBASIS_END; i++) PetscCall(LMBasisDestroy(&lmvm->basis[i])); in MatReset_LMVM()
70 for (PetscInt i = 0; i < LMBASIS_END; i++) PetscCall(LMBasisReset(lmvm->basis[i])); in MatReset_LMVM()
96 …if (!lmvm->basis[LMBASIS_S]) PetscCall(LMBasisCreate(lmvm->Xprev, lmvm->m, &lmvm->basis[LMBASIS_S]… in MatLMVMAllocateBases()
97 …if (!lmvm->basis[LMBASIS_Y]) PetscCall(LMBasisCreate(lmvm->Fprev, lmvm->m, &lmvm->basis[LMBASIS_Y]… in MatLMVMAllocateBases()
139 PetscCall(LMBasisGetNextVec(lmvm->basis[LMBASIS_S], &s_k)); in MatUpdateKernel_LMVM()
141 PetscCall(LMBasisRestoreNextVec(lmvm->basis[LMBASIS_S], &s_k)); in MatUpdateKernel_LMVM()
143 PetscCall(LMBasisGetNextVec(lmvm->basis[LMBASIS_Y], &y_k)); in MatUpdateKernel_LMVM()
145 PetscCall(LMBasisRestoreNextVec(lmvm->basis[LMBASIS_Y], &y_k)); in MatUpdateKernel_LMVM()
148 …PetscAssert(lmvm->k == lmvm->basis[LMBASIS_S]->k, PetscObjectComm((PetscObject)B), PETSC_ERR_PLIB,… in MatUpdateKernel_LMVM()
149 …PetscAssert(lmvm->k == lmvm->basis[LMBASIS_Y]->k, PetscObjectComm((PetscObject)B), PETSC_ERR_PLIB,… in MatUpdateKernel_LMVM()
[all …]
|
| /petsc/src/dm/dt/dualspace/impls/lagrange/tutorials/output/ |
| H A D | ex1_2.out | 5 Dual basis vector 0
8 Dual basis vector 1
11 Dual basis vector 2
14 Dual basis vector 3
17 Dual basis vector 4
20 Dual basis vector 5
23 Dual basis vector 6
26 Dual basis vector 7
29 Dual basis vector 8
32 Dual basis vector 9
[all …]
|
| H A D | ex1_4.out | 5 Dual basis vector 0
8 Dual basis vector 1
11 Dual basis vector 2
14 Dual basis vector 3
17 Dual basis vector 4
20 Dual basis vector 5
23 Dual basis vector 6
26 Dual basis vector 7
29 Dual basis vector 8
32 Dual basis vector 9
|
| H A D | ex1_6.out | 5 Dual basis vector 0
8 Dual basis vector 1
11 Dual basis vector 2
14 Dual basis vector 3
17 Dual basis vector 4
20 Dual basis vector 5
23 Dual basis vector 6
26 Dual basis vector 7
29 Dual basis vector 8
32 Dual basis vector 9
|
| H A D | ex1_7.out | 5 Dual basis vector 0
8 Dual basis vector 1
11 Dual basis vector 2
14 Dual basis vector 3
17 Dual basis vector 4
20 Dual basis vector 5
23 Dual basis vector 6
26 Dual basis vector 7
29 Dual basis vector 8
32 Dual basis vector 9
|
| H A D | ex1_5.out | 5 Dual basis vector 0
8 Dual basis vector 1
11 Dual basis vector 2
14 Dual basis vector 3
17 Dual basis vector 4
20 Dual basis vector 5
23 Dual basis vector 6
26 Dual basis vector 7
29 Dual basis vector 8
32 Dual basis vector 9
|
| H A D | ex1_3.out | 5 Dual basis vector 0
8 Dual basis vector 1
11 Dual basis vector 2
14 Dual basis vector 3
17 Dual basis vector 4
20 Dual basis vector 5
23 Dual basis vector 6
26 Dual basis vector 7
29 Dual basis vector 8
32 Dual basis vector 9
[all …]
|
| H A D | ex1_0.out | 5 Dual basis vector 0
8 Dual basis vector 1
11 Dual basis vector 2
14 Dual basis vector 3
17 Dual basis vector 4
20 Dual basis vector 5
|
| H A D | ex1_1.out | 5 Dual basis vector 0
8 Dual basis vector 1
11 Dual basis vector 2
14 Dual basis vector 3
|
| /petsc/src/dm/impls/swarm/ |
| H A D | swarmpic_da.c | 93 PetscReal **basis; in private_DMSwarmInsertPointsUsingCellDM_DA_Q1() local
131 PetscCall(PetscMalloc1(npoints_q, &basis)); in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
133 PetscCall(PetscMalloc1(npe, &basis[q])); in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
137 basis[q][0] = 0.5 * (1.0 - xi[dim * q + 0]); in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
138 basis[q][1] = 0.5 * (1.0 + xi[dim * q + 0]); in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
141 basis[q][0] = 0.25 * (1.0 - xi[dim * q + 0]) * (1.0 - xi[dim * q + 1]); in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
142 basis[q][1] = 0.25 * (1.0 + xi[dim * q + 0]) * (1.0 - xi[dim * q + 1]); in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
143 basis[q][2] = 0.25 * (1.0 + xi[dim * q + 0]) * (1.0 + xi[dim * q + 1]); in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
144 basis[q][3] = 0.25 * (1.0 - xi[dim * q + 0]) * (1.0 + xi[dim * q + 1]); in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
148 … basis[q][0] = 0.125 * (1.0 - xi[dim * q + 0]) * (1.0 - xi[dim * q + 1]) * (1.0 - xi[dim * q + 2]); in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
[all …]
|
| H A D | swarmpic_plex.c | 162 PetscReal **basis; in private_DMSwarmInsertPointsUsingCellDM_PLEX2D_Regular() local
201 PetscCall(PetscMalloc1(npoints_q, &basis)); in private_DMSwarmInsertPointsUsingCellDM_PLEX2D_Regular()
203 PetscCall(PetscMalloc1(npe, &basis[q])); in private_DMSwarmInsertPointsUsingCellDM_PLEX2D_Regular()
205 basis[q][0] = 1.0 - xi[dim * q + 0] - xi[dim * q + 1]; in private_DMSwarmInsertPointsUsingCellDM_PLEX2D_Regular()
206 basis[q][1] = xi[dim * q + 0]; in private_DMSwarmInsertPointsUsingCellDM_PLEX2D_Regular()
207 basis[q][2] = xi[dim * q + 1]; in private_DMSwarmInsertPointsUsingCellDM_PLEX2D_Regular()
233 …for (k = 0; k < npe; k++) swarm_coor[dim * pcnt + d] += basis[q][k] * PetscRealPart(elcoor[dim * k… in private_DMSwarmInsertPointsUsingCellDM_PLEX2D_Regular()
244 for (q = 0; q < npoints_q; q++) PetscCall(PetscFree(basis[q])); in private_DMSwarmInsertPointsUsingCellDM_PLEX2D_Regular()
245 PetscCall(PetscFree(basis)); in private_DMSwarmInsertPointsUsingCellDM_PLEX2D_Regular()
|
| /petsc/src/dm/dt/tests/output/ |
| H A D | ex16_0.out | 16 PetscSpace Object: vector basis space 1 MPI process
33 Vector basis to original basis:
39 PetscSpace Object: vector basis space 1 MPI process
56 Vector basis to original basis:
62 PetscSpace Object: vector basis space 1 MPI process
79 Vector basis to original basis:
85 PetscSpace Object: vector basis space 1 MPI process
102 Vector basis to original basis:
108 PetscSpace Object: vector basis space 1 MPI process
125 Vector basis to original basis:
[all …]
|
| /petsc/src/binding/petsc4py/test/ |
| H A D | test_nsp.py | 22 basis = [u1, u2]
23 nullsp = PETSc.NullSpace().create(False, basis, comm=PETSc.COMM_SELF)
24 self.basis = basis
28 self.basis = None
|
| /petsc/src/dm/dt/fv/interface/ceed/ |
| H A D | fvceed.c | 45 PetscErrorCode PetscFVGetCeedBasis(PetscFV fv, CeedBasis *basis) in PetscFVGetCeedBasis() argument
52 PetscAssertPointer(basis, 2); in PetscFVGetCeedBasis()
60 *basis = fv->ceedBasis; in PetscFVGetCeedBasis()
|
| /petsc/src/dm/dt/fe/interface/ceed/ |
| H A D | feceed.c | 45 PetscErrorCode PetscFEGetCeedBasis(PetscFE fe, CeedBasis *basis) in PetscFEGetCeedBasis() argument
53 PetscAssertPointer(basis, 2); in PetscFEGetCeedBasis()
63 *basis = fe->ceedBasis; in PetscFEGetCeedBasis()
|
| /petsc/doc/manualpages/MANSECHeaders/ |
| H A D | DMSwarm | 3 …. It also can conservatively project fields between the particle basis and an FEM basis on a `DMPL…
|
| /petsc/src/ksp/ksp/tutorials/output/ |
| H A D | ex71_bddc_elast_4lev.out | 19 Use change of basis on local edge nodes: 0
20 Use change of basis on local face nodes: 0
21 User defined change of basis matrix: 0
22 Has change of basis matrix: 0
40 Symmetric computation of primal basis functions: 1
132 Use change of basis on local edge nodes: 0
133 Use change of basis on local face nodes: 0
134 User defined change of basis matrix: 0
135 Has change of basis matrix: 0
153 Symmetric computation of primal basis functions: 1
[all …]
|
| H A D | ex71_bddc_elast_3lev.out | 19 Use change of basis on local edge nodes: 0
20 Use change of basis on local face nodes: 0
21 User defined change of basis matrix: 0
22 Has change of basis matrix: 0
40 Symmetric computation of primal basis functions: 1
132 Use change of basis on local edge nodes: 0
133 Use change of basis on local face nodes: 0
134 User defined change of basis matrix: 0
135 Has change of basis matrix: 0
153 Symmetric computation of primal basis functions: 1
|
| H A D | ex59_bddc_fetidp_ml_2.out | 19 Use change of basis on local edge nodes: 0
20 Use change of basis on local face nodes: 0
21 User defined change of basis matrix: 0
22 Has change of basis matrix: 0
40 Symmetric computation of primal basis functions: 1
130 Use change of basis on local edge nodes: 0
131 Use change of basis on local face nodes: 0
132 User defined change of basis matrix: 0
133 Has change of basis matrix: 0
151 Symmetric computation of primal basis functions: 1
|
| H A D | ex71_bddc_elast_3lev_alt.out | 19 Use change of basis on local edge nodes: 0
20 Use change of basis on local face nodes: 0
21 User defined change of basis matrix: 0
22 Has change of basis matrix: 0
40 Symmetric computation of primal basis functions: 1
132 Use change of basis on local edge nodes: 0
133 Use change of basis on local face nodes: 0
134 User defined change of basis matrix: 0
135 Has change of basis matrix: 0
153 Symmetric computation of primal basis functions: 1
|
| H A D | ex59_bddc_fetidp_ml_1.out | 19 Use change of basis on local edge nodes: 0
20 Use change of basis on local face nodes: 0
21 User defined change of basis matrix: 0
22 Has change of basis matrix: 0
40 Symmetric computation of primal basis functions: 1
130 Use change of basis on local edge nodes: 0
131 Use change of basis on local face nodes: 0
132 User defined change of basis matrix: 0
133 Has change of basis matrix: 0
151 Symmetric computation of primal basis functions: 1
|
| /petsc/src/dm/interface/ |
| H A D | dmceed.c | 156 PetscCall(PetscFEGetCeedBasis(fe, &sd->basis)); in DMCeedCreateGeometry()
157 PetscCall(CeedBasisGetNumQuadraturePoints(sd->basis, &Nq)); in DMCeedCreateGeometry()
180 PetscCallCEED(CeedOperatorSetField(sd->op, "x", sd->er, sd->basis, CEED_VECTOR_ACTIVE)); in DMCeedCreateGeometry()
181 PetscCallCEED(CeedOperatorSetField(sd->op, "dx", sd->er, sd->basis, CEED_VECTOR_ACTIVE)); in DMCeedCreateGeometry()
182 …PetscCallCEED(CeedOperatorSetField(sd->op, "weight", CEED_ELEMRESTRICTION_NONE, sd->basis, CEED_VE… in DMCeedCreateGeometry()
212 PetscCall(PetscFEGetCeedBasis(fe, &sd->basis)); in DMCeedCreate_Internal()
214 PetscCall(CeedBasisGetNumQuadraturePoints(sd->basis, &Nq)); in DMCeedCreate_Internal()
228 PetscCallCEED(CeedBasisGetNumQuadraturePoints(sd->geom->basis, &Nqx)); in DMCeedCreate_Internal()
243 PetscCallCEED(CeedOperatorSetField(sd->op, "u", sd->er, sd->basis, CEED_VECTOR_ACTIVE)); in DMCeedCreate_Internal()
244 PetscCallCEED(CeedOperatorSetField(sd->op, "du", sd->er, sd->basis, CEED_VECTOR_ACTIVE)); in DMCeedCreate_Internal()
[all …]
|
| /petsc/src/snes/tutorials/output/ |
| H A D | ex12_p4est_full_q2_nonconformal_parallel_bddcfas.out | 50 Use change of basis on local edge nodes: 0
51 Use change of basis on local face nodes: 0
52 User defined change of basis matrix: 0
53 Has change of basis matrix: 0
71 Symmetric computation of primal basis functions: 1
225 Use change of basis on local edge nodes: 0
226 Use change of basis on local face nodes: 0
227 User defined change of basis matrix: 0
228 Has change of basis matrix: 0
246 Symmetric computation of primal basis functions: 1
[all …]
|