| /petsc/src/tao/constrained/impls/almm/ |
| H A D | almm.c | 45 if (tao->eq_constrained) PetscCall(VecDot(auglag->Ce, auglag->Ce, &auglag->cenorm)); in TaoSolve_ALMM()
85 PetscCall(VecAXPY(auglag->Ye, auglag->mu, auglag->Ce)); in TaoSolve_ALMM()
161 auglag->Ce = tao->constraints_equality; in TaoSetUp_ALMM()
164 PetscCall(VecDuplicate(auglag->Ce, &auglag->Ye)); in TaoSetUp_ALMM()
166 if (!auglag->Cework) PetscCall(VecDuplicate(auglag->Ce, &auglag->Cework)); in TaoSetUp_ALMM()
214 if (!auglag->C) auglag->C = auglag->Ce; in TaoSetUp_ALMM()
582 if (tao->eq_constrained) PetscCall(VecNorm(auglag->Ce, NORM_INFINITY, &auglag->cenorm)); in TaoALMMComputeOptimalityNorms_Private()
611 PetscCall(TaoComputeEqualityConstraints(tao, auglag->Px, auglag->Ce)); in TaoALMMEvaluateIterate_Private()
632 PetscCall(TaoALMMCombineDual_Private(tao, auglag->Ce, auglag->Ci, auglag->C)); in TaoALMMEvaluateIterate_Private()
652 PetscCall(VecWAXPY(auglag->Cework, 1.0 / auglag->mu, auglag->Ye, auglag->Ce)); in TaoALMMComputePHRLagAndGradient_Private()
[all …]
|
| H A D | almm.h | 13 Vec Px, LgradX, Ce, Ci, G; /* aliased vectors (do not destroy!) */ member
|
| /petsc/src/ksp/ksp/tests/ |
| H A D | ex21.c | 161 Mat B, C, Ce, Cee, Cer; in main() local
179 PetscCall(MatComputeOperator(C, MATAIJ, &Ce)); in main()
183 PetscCall(MatNormDifference(Ce, Cee, &err)); in main()
188 PetscCall(MatDestroy(&Ce)); in main()
|
| /petsc/src/mat/impls/dense/mpi/ |
| H A D | mpidense.h | 9 Mat Ae, Be, Ce; /* matrix in Elemental format */ member
|
| H A D | mpidense.c | 2420 PetscCall(MatDestroy(&ab->Ce)); in MatProductCtxDestroy_MatMatMult_MPIDense_MPIDense()
2437 if (ab->Ae && ab->Ce) { in MatMatMultNumeric_MPIDense_MPIDense()
2441 PetscCall(MatMatMultNumeric_Elemental(ab->Ae, ab->Be, ab->Ce)); in MatMatMultNumeric_MPIDense_MPIDense()
2442 PetscCall(MatConvert(ab->Ce, MATMPIDENSE, MAT_REUSE_MATRIX, &C)); in MatMatMultNumeric_MPIDense_MPIDense()
2562 PetscCall(MatCreate(PetscObjectComm((PetscObject)C), &ab->Ce)); in MatMatMultSymbolic_MPIDense_MPIDense()
2563 PetscCall(MatMatMultSymbolic_Elemental(ab->Ae, ab->Be, fill, ab->Ce)); in MatMatMultSymbolic_MPIDense_MPIDense()
2571 ab->Ce = NULL; in MatMatMultSymbolic_MPIDense_MPIDense()
|
| /petsc/src/mat/impls/elemental/ |
| H A D | matelem.cxx | 316 PetscErrorCode MatMatMultSymbolic_Elemental(Mat A, Mat B, PetscReal, Mat Ce) in MatMatMultSymbolic_Elemental() argument
319 PetscCall(MatSetSizes(Ce, A->rmap->n, B->cmap->n, PETSC_DECIDE, PETSC_DECIDE)); in MatMatMultSymbolic_Elemental()
320 PetscCall(MatSetType(Ce, MATELEMENTAL)); in MatMatMultSymbolic_Elemental()
321 PetscCall(MatSetUp(Ce)); in MatMatMultSymbolic_Elemental()
322 Ce->ops->matmultnumeric = MatMatMultNumeric_Elemental; in MatMatMultSymbolic_Elemental()
386 Mat Be, Ce; in MatMatMultNumeric_Elemental_MPIDense() local
390 PetscCall(MatMatMult(A, Be, MAT_INITIAL_MATRIX, PETSC_DETERMINE, &Ce)); in MatMatMultNumeric_Elemental_MPIDense()
391 PetscCall(MatConvert(Ce, MATMPIDENSE, MAT_REUSE_MATRIX, &C)); in MatMatMultNumeric_Elemental_MPIDense()
393 PetscCall(MatDestroy(&Ce)); in MatMatMultNumeric_Elemental_MPIDense()
|
| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex70.c | 347 PetscScalar Ce[NODES_PER_EL * P_DOFS * NODES_PER_EL * P_DOFS]; in AssembleStokes_A() local
375 PetscCall(PetscMemzero(Ce, sizeof(Ce))); in AssembleStokes_A()
381 BForm_Stabilisation(Ce, el_coords, prop_eta); in AssembleStokes_A()
388 …PetscCall(MatSetValuesLocal(A, NODES_PER_EL * P_DOFS, p_eqn, NODES_PER_EL * P_DOFS, p_eqn, Ce, ADD… in AssembleStokes_A()
410 PetscScalar Ce[NODES_PER_EL * P_DOFS * NODES_PER_EL * P_DOFS]; in AssembleStokes_PC() local
438 PetscCall(PetscMemzero(Ce, sizeof(Ce))); in AssembleStokes_PC()
443 BForm_ScaledMassMatrix(Ce, el_coords, prop_eta); in AssembleStokes_PC()
450 …PetscCall(MatSetValuesLocal(A, NODES_PER_EL * P_DOFS, p_eqn, NODES_PER_EL * P_DOFS, p_eqn, Ce, ADD… in AssembleStokes_PC()
|
| H A D | ex43.c | 715 PetscScalar Ce[NODES_PER_EL * P_DOFS * NODES_PER_EL * P_DOFS]; in AssembleA_Stokes() local
747 PetscCall(PetscMemzero(Ce, sizeof(Ce))); in AssembleA_Stokes()
753 FormStabilisationOperatorQ1(Ce, el_coords, prop_eta); in AssembleA_Stokes()
760 …atSetValuesStencil(A, NODES_PER_EL * P_DOFS, p_eqn, NODES_PER_EL * P_DOFS, p_eqn, Ce, ADD_VALUES)); in AssembleA_Stokes()
785 PetscScalar Ce[NODES_PER_EL * P_DOFS * NODES_PER_EL * P_DOFS]; in AssembleA_PCStokes() local
817 PetscCall(PetscMemzero(Ce, sizeof(Ce))); in AssembleA_PCStokes()
822 FormScaledMassMatrixOperatorQ1(Ce, el_coords, prop_eta); in AssembleA_PCStokes()
828 …atSetValuesStencil(A, NODES_PER_EL * P_DOFS, p_eqn, NODES_PER_EL * P_DOFS, p_eqn, Ce, ADD_VALUES)); in AssembleA_PCStokes()
|
| H A D | ex42.c | 908 PetscScalar Ce[NODES_PER_EL * P_DOFS * NODES_PER_EL * P_DOFS]; in AssembleA_Stokes() local
936 PetscCall(PetscMemzero(Ce, sizeof(Ce))); in AssembleA_Stokes()
944 FormStabilisationOperatorQ13D(Ce, el_coords, prop_eta); in AssembleA_Stokes()
971 …atSetValuesStencil(A, NODES_PER_EL * P_DOFS, p_eqn, NODES_PER_EL * P_DOFS, p_eqn, Ce, ADD_VALUES)); in AssembleA_Stokes()
994 PetscScalar Ce[NODES_PER_EL * P_DOFS * NODES_PER_EL * P_DOFS]; in AssembleA_PCStokes() local
1023 PetscCall(PetscMemzero(Ce, sizeof(Ce))); in AssembleA_PCStokes()
1029 FormScaledMassMatrixOperatorQ13D(Ce, el_coords, prop_eta); in AssembleA_PCStokes()
1056 …atSetValuesStencil(A, NODES_PER_EL * P_DOFS, p_eqn, NODES_PER_EL * P_DOFS, p_eqn, Ce, ADD_VALUES)); in AssembleA_PCStokes()
|