| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex21.c | 27 Mat A, Ae = NULL, Ad = NULL; in main() local 94 PetscCall(MatComputeOperator(A, MATDENSE, &Ae)); in main() 95 PetscCall(MatAXPY(Ae, -1.0, Ad, SAME_NONZERO_PATTERN)); in main() 96 PetscCall(MatGetDiagonal(Ae, d)); in main() 97 PetscCall(MatViewFromOptions(Ae, NULL, "-A_view")); in main() 98 PetscCall(MatViewFromOptions(Ae, NULL, "-D_view")); in main() 99 PetscCall(MatNorm(Ae, NORM_FROBENIUS, &norm)); in main() 103 PetscCall(MatDestroy(&Ae)); in main() 133 PetscCall(MatComputeOperator(A, MATDENSE, &Ae)); in main() 134 PetscCall(MatAXPY(Ae, -1.0, Ad, SAME_NONZERO_PATTERN)); in main() [all …]
|
| H A D | ex42.c | 905 PetscScalar Ae[NODES_PER_EL * U_DOFS * NODES_PER_EL * U_DOFS]; in AssembleA_Stokes() local 933 PetscCall(PetscMemzero(Ae, sizeof(Ae))); in AssembleA_Stokes() 939 FormStressOperatorQ13D(Ae, el_coords, prop_eta); in AssembleA_Stokes() 951 _ZERO_ROWCOL_i(Ae, 3 * n); in AssembleA_Stokes() 957 _ZERO_ROWCOL_i(Ae, 3 * n + 1); in AssembleA_Stokes() 963 _ZERO_ROWCOL_i(Ae, 3 * n + 2); in AssembleA_Stokes() 968 …atSetValuesStencil(A, NODES_PER_EL * U_DOFS, u_eqn, NODES_PER_EL * U_DOFS, u_eqn, Ae, ADD_VALUES)); in AssembleA_Stokes() 991 PetscScalar Ae[NODES_PER_EL * U_DOFS * NODES_PER_EL * U_DOFS]; in AssembleA_PCStokes() local 1020 PetscCall(PetscMemzero(Ae, sizeof(Ae))); in AssembleA_PCStokes() 1026 FormStressOperatorQ13D(Ae, el_coords, prop_eta); in AssembleA_PCStokes() [all …]
|
| H A D | ex70.c | 344 PetscScalar Ae[NODES_PER_EL * U_DOFS * NODES_PER_EL * U_DOFS]; in AssembleStokes_A() local 372 PetscCall(PetscMemzero(Ae, sizeof(Ae))); in AssembleStokes_A() 378 BForm_DivT(Ae, el_coords, prop_eta); in AssembleStokes_A() 385 …PetscCall(MatSetValuesLocal(A, NODES_PER_EL * U_DOFS, u_eqn, NODES_PER_EL * U_DOFS, u_eqn, Ae, ADD… in AssembleStokes_A() 407 PetscScalar Ae[NODES_PER_EL * U_DOFS * NODES_PER_EL * U_DOFS]; in AssembleStokes_PC() local 435 PetscCall(PetscMemzero(Ae, sizeof(Ae))); in AssembleStokes_PC() 441 BForm_DivT(Ae, el_coords, prop_eta); in AssembleStokes_PC() 447 …PetscCall(MatSetValuesLocal(A, NODES_PER_EL * U_DOFS, u_eqn, NODES_PER_EL * U_DOFS, u_eqn, Ae, ADD… in AssembleStokes_PC()
|
| H A D | ex43.c | 712 PetscScalar Ae[NODES_PER_EL * U_DOFS * NODES_PER_EL * U_DOFS]; in AssembleA_Stokes() local 744 PetscCall(PetscMemzero(Ae, sizeof(Ae))); in AssembleA_Stokes() 750 FormStressOperatorQ1(Ae, el_coords, prop_eta); in AssembleA_Stokes() 757 …atSetValuesStencil(A, NODES_PER_EL * U_DOFS, u_eqn, NODES_PER_EL * U_DOFS, u_eqn, Ae, ADD_VALUES)); in AssembleA_Stokes() 782 PetscScalar Ae[NODES_PER_EL * U_DOFS * NODES_PER_EL * U_DOFS]; in AssembleA_PCStokes() local 814 PetscCall(PetscMemzero(Ae, sizeof(Ae))); in AssembleA_PCStokes() 820 FormStressOperatorQ1(Ae, el_coords, prop_eta); in AssembleA_PCStokes() 826 …atSetValuesStencil(A, NODES_PER_EL * U_DOFS, u_eqn, NODES_PER_EL * U_DOFS, u_eqn, Ae, ADD_VALUES)); in AssembleA_PCStokes()
|
| H A D | ex49.c | 574 PetscScalar Ae[NODES_PER_EL * U_DOFS * NODES_PER_EL * U_DOFS]; in AssembleA_Elasticity() local 604 PetscCall(PetscMemzero(Ae, sizeof(Ae))); in AssembleA_Elasticity() 607 FormStressOperatorQ1(Ae, el_coords, prop_E, prop_nu); in AssembleA_Elasticity() 611 …atSetValuesStencil(A, NODES_PER_EL * U_DOFS, u_eqn, NODES_PER_EL * U_DOFS, u_eqn, Ae, ADD_VALUES)); in AssembleA_Elasticity()
|
| /petsc/src/mat/tests/ |
| H A D | ex222.c | 7 Mat A, Ae, Aet; in main() local 49 PetscCall(MatComputeOperator(A, etype, &Ae)); in main() 50 PetscCall(PetscObjectSetName((PetscObject)Ae, "Explicit matrix")); in main() 51 PetscCall(MatViewFromOptions(Ae, NULL, "-view_expl")); in main() 60 PetscCall(MatAXPY(A2, -1.0, Ae, DIFFERENT_NONZERO_PATTERN)); in main() 71 PetscCall(MatDestroy(&Ae)); in main()
|
| H A D | ex174.cxx | 12 Mat A, Ae, B, Be; in main() local 63 Ae = A; in main() 86 PetscCall(MatConvert(A, MATELEMENTAL, MAT_INITIAL_MATRIX, &Ae)); in main() 90 PetscCall(MatMultEqual(A, Ae, 5, &flg)); in main() 97 PetscCall(MatDestroy(&Ae)); in main()
|
| H A D | ex244.cxx | 11 Mat A, Ae, B, Be; in main() local 62 Ae = A; in main() 85 PetscCall(MatConvert(A, MATSCALAPACK, MAT_INITIAL_MATRIX, &Ae)); in main() 89 PetscCall(MatMultEqual(A, Ae, 5, &flg)); in main() 96 PetscCall(MatDestroy(&Ae)); in main()
|
| H A D | ex125.c | 47 Mat A, Ae, RHS = NULL, RHS1 = NULL, C, F, X; in main() local 189 PetscCall(MatComputeOperator(A, MATAIJ, &Ae)); in main() 192 Ae = A; in main() 199 PetscCall(MatViewFromOptions(Ae, NULL, "-A_view_expl")); in main() 399 PetscCall(MatMatMult(Ae, C, MAT_INITIAL_MATRIX, 2.0, &RHS)); in main() 401 PetscCall(MatMatMult(Ae, C, MAT_REUSE_MATRIX, 2.0, &RHS)); in main() 447 PetscCall(MatTransposeMatMult(Ae, C, MAT_INITIAL_MATRIX, 2.0, &RHS1)); in main() 449 PetscCall(MatTransposeMatMult(Ae, C, MAT_REUSE_MATRIX, 2.0, &RHS1)); in main() 487 PetscCall(MatMult(Ae, x, b)); in main() 497 PetscCall(MatMult(Ae, x, u)); /* u = A*x */ in main() [all …]
|
| /petsc/src/tao/constrained/tutorials/ |
| H A D | ex1.c | 42 Mat Ae, Ai, H; member 87 …PetscCall(TaoSetJacobianEqualityRoutine(tao, user.Ae, user.Ae, FormEqualityJacobian, (void *)&user… in main() 130 PetscCall(FormEqualityJacobian(tao, user.x, user.Ae, user.Ae, (void *)&user)); in main() 131 PetscCall(MatCreateVecs(user.Ae, NULL, &CE)); in main() 135 PetscCall(MatView(user.Ae, PETSC_VIEWER_STDOUT_WORLD)); in main() 215 PetscCall(MatCreate(PETSC_COMM_WORLD, &user->Ae)); in InitializeProblem() 216 PetscCall(MatSetSizes(user->Ae, neloc, nloc, user->ne, user->n)); in InitializeProblem() 217 PetscCall(MatSetFromOptions(user->Ae)); in InitializeProblem() 218 PetscCall(MatSetUp(user->Ae)); in InitializeProblem() 236 if (!user->noeqflag) PetscCall(MatDestroy(&user->Ae)); in DestroyProblem()
|
| /petsc/src/tao/constrained/impls/almm/ |
| H A D | almm.h | 12 Mat Ae, Ai; /* aliased constraint Jacobians (do not destroy!) */ member
|
| H A D | almm.c | 162 auglag->Ae = tao->jacobian_equality; in TaoSetUp_ALMM() 612 PetscCall(TaoComputeJacobianEquality(tao, auglag->Px, auglag->Ae, auglag->Ae)); in TaoALMMEvaluateIterate_Private() 656 PetscCall(MatMultTransposeAdd(auglag->Ae, auglag->Cework, auglag->LgradX, auglag->LgradX)); in TaoALMMComputePHRLagAndGradient_Private() 695 PetscCall(MatMultTransposeAdd(auglag->Ae, auglag->Ye, auglag->LgradX, auglag->LgradX)); in TaoALMMComputeAugLagAndGradient_Private() 697 PetscCall(MatMultTranspose(auglag->Ae, auglag->Ce, auglag->Xwork)); in TaoALMMComputeAugLagAndGradient_Private()
|
| /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 | 2421 PetscCall(MatDestroy(&ab->Ae)); in MatProductCtxDestroy_MatMatMult_MPIDense_MPIDense() 2437 if (ab->Ae && ab->Ce) { in MatMatMultNumeric_MPIDense_MPIDense() 2439 PetscCall(MatConvert_MPIDense_Elemental(A, MATELEMENTAL, MAT_REUSE_MATRIX, &ab->Ae)); in MatMatMultNumeric_MPIDense_MPIDense() 2441 PetscCall(MatMatMultNumeric_Elemental(ab->Ae, ab->Be, ab->Ce)); in MatMatMultNumeric_MPIDense_MPIDense() 2549 PetscCall(MatCreate(PetscObjectComm((PetscObject)A), &ab->Ae)); in MatMatMultSymbolic_MPIDense_MPIDense() 2550 PetscCall(MatSetSizes(ab->Ae, PETSC_DECIDE, PETSC_DECIDE, A->rmap->N, A->cmap->N)); in MatMatMultSymbolic_MPIDense_MPIDense() 2551 PetscCall(MatSetType(ab->Ae, MATELEMENTAL)); in MatMatMultSymbolic_MPIDense_MPIDense() 2552 PetscCall(MatSetUp(ab->Ae)); in MatMatMultSymbolic_MPIDense_MPIDense() 2553 PetscCall(MatSetOption(ab->Ae, MAT_ROW_ORIENTED, PETSC_FALSE)); in MatMatMultSymbolic_MPIDense_MPIDense() 2563 PetscCall(MatMatMultSymbolic_Elemental(ab->Ae, ab->Be, fill, ab->Ce)); in MatMatMultSymbolic_MPIDense_MPIDense() [all …]
|
| /petsc/src/mat/impls/h2opus/cuda/ |
| H A D | math2opus.cu | 955 Mat E, Ae; in MatAssemblyEnd_H2OPUS() local 960 Ae = a->sampler->GetSamplingMat(); in MatAssemblyEnd_H2OPUS() 963 PetscCall(MatAXPY(E, -1.0, Ae, DIFFERENT_NONZERO_PATTERN)); in MatAssemblyEnd_H2OPUS() 972 PetscCall(MatComputeOperator(Ae, MATAIJ, &eAe)); in MatAssemblyEnd_H2OPUS() 987 PetscCall(MatGetOperation(Ae, MATOP_NORM, &normfunc)); in MatAssemblyEnd_H2OPUS() 988 PetscCall(MatSetOperation(Ae, MATOP_NORM, (PetscErrorCodeFn *)MatNorm_H2OPUS)); in MatAssemblyEnd_H2OPUS() 989 PetscCall(MatNorm(Ae, NORM_1, &n1A)); in MatAssemblyEnd_H2OPUS() 990 PetscCall(MatNorm(Ae, NORM_INFINITY, &niA)); in MatAssemblyEnd_H2OPUS() 991 PetscCall(MatNorm(Ae, NORM_2, &n2A)); in MatAssemblyEnd_H2OPUS() 995 PetscCall(MatSetOperation(Ae, MATOP_NORM, normfunc)); in MatAssemblyEnd_H2OPUS()
|
| /petsc/src/mat/impls/elemental/ |
| H A D | matelem.cxx | 1152 Mat Adense, Ae; in MatLoad_Elemental() local 1160 PetscCall(MatConvert(Adense, MATELEMENTAL, MAT_INITIAL_MATRIX, &Ae)); in MatLoad_Elemental() 1162 PetscCall(MatHeaderReplace(newMat, &Ae)); in MatLoad_Elemental()
|