| /petsc/src/tao/unconstrained/tutorials/ |
| H A D | adjointreads.m | 9 %Grad --- gradient,
25 % plot(xg,Grad,'bs','Markersize',10); drawnow
37 % plot(xg,Grad,'go-','Markersize',10); drawnow
123 % plot(xg,Grad,'k*-');
125 % plot(xg,Grad,'ro-');
134 % legend('Grad at it=0','Grad at it=1')
142 % plot(xg,Grad,'k*-','LineWidth',2,'Markersize',10);
153 % errgrad=max(abs(gradj./Mass-Grad))
|
| /petsc/src/tao/pde_constrained/tutorials/ |
| H A D | elliptic.c | 29 Mat Grad; member
263 PetscCall(MatPtAP(user->Diag, user->Grad, MAT_REUSE_MATRIX, 1.0, &user->DSG)); in FormJacobianState()
413 PetscCall(MatMult(user->Grad, user->y, user->Twork)); in DesignMatMult()
417 PetscCall(MatMultTranspose(user->Grad, user->Swork, Y)); in DesignMatMult()
423 PetscCall(MatMult(user->Grad, user->suby, user->Twork)); in DesignMatMult()
425 PetscCall(MatMultTranspose(user->Grad, user->Twork, user->subq)); in DesignMatMult()
455 PetscCall(MatMult(user->Grad, user->subq, user->Swork)); in DesignMatMultTranspose()
458 PetscCall(MatMult(user->Grad, user->suby, user->Twork)); in DesignMatMultTranspose()
487 PetscCall(MatMult(user->Grad, user->y, user->Swork)); in FormConstraints()
489 PetscCall(MatMultTranspose(user->Grad, user->Swork, C)); in FormConstraints()
[all …]
|
| H A D | parabolic.c | 30 Mat Grad; member
343 … PetscCall(MatMatMult(user->Divwork, user->Grad, MAT_INITIAL_MATRIX, PETSC_DETERMINE, &user->DSG)); in FormJacobianState()
405 PetscCall(MatMult(user->Grad, X, user->Swork)); in StateMatBlockMult()
440 PetscCall(MatMult(user->Grad, user->yi[i], user->Twork)); in DesignMatMult()
472 PetscCall(MatMult(user->Grad, user->yiwork[i], user->Swork)); in DesignMatMultTranspose()
475 PetscCall(MatMult(user->Grad, user->yi[i], user->Twork)); in DesignMatMultTranspose()
779 PetscCall(MatCreate(PETSC_COMM_WORLD, &user->Grad)); in ParabolicInitialize()
780 PetscCall(MatSetSizes(user->Grad, PETSC_DECIDE, PETSC_DECIDE, m, n)); in ParabolicInitialize()
781 PetscCall(MatSetFromOptions(user->Grad)); in ParabolicInitialize()
782 PetscCall(MatMPIAIJSetPreallocation(user->Grad, 2, NULL, 2, NULL)); in ParabolicInitialize()
[all …]
|
| H A D | hyperbolic.c | 26 Mat Grad, Gradxy[2]; member
400 PetscCall(MatMult(user->Grad, X, user->uiwork[i])); in StateMatBlockMultTranspose()
439 PetscCall(MatMult(user->Grad, user->yiwork[i], user->uiwork[i])); in DesignMatMultTranspose()
698 PetscCall(MatCreate(PETSC_COMM_WORLD, &user->Grad)); in HyperbolicInitialize()
699 PetscCall(MatSetSizes(user->Grad, PETSC_DECIDE, PETSC_DECIDE, 2 * n, n)); in HyperbolicInitialize()
700 PetscCall(MatSetFromOptions(user->Grad)); in HyperbolicInitialize()
701 PetscCall(MatMPIAIJSetPreallocation(user->Grad, 3, NULL, 3, NULL)); in HyperbolicInitialize()
702 PetscCall(MatSeqAIJSetPreallocation(user->Grad, 3, NULL)); in HyperbolicInitialize()
703 PetscCall(MatGetOwnershipRange(user->Grad, &istart, &iend)); in HyperbolicInitialize()
709 PetscCall(MatSetValues(user->Grad, 1, &i, 1, &j, &half_hinv, INSERT_VALUES)); in HyperbolicInitialize()
[all …]
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexfem.c | 1039 …uesRiemann(DM dm, PetscReal time, Vec faceGeometry, Vec cellGeometry, Vec Grad, PetscInt field, Pe… in DMPlexInsertBoundaryValuesRiemann() argument
1063 if (Grad) { in DMPlexInsertBoundaryValuesRiemann()
1067 PetscCall(VecGetDM(Grad, &dmGrad)); in DMPlexInsertBoundaryValuesRiemann()
1068 PetscCall(VecGetArrayRead(Grad, &grad)); in DMPlexInsertBoundaryValuesRiemann()
1091 if (Grad) { in DMPlexInsertBoundaryValuesRiemann()
1124 if (Grad) { in DMPlexInsertBoundaryValuesRiemann()
1126 PetscCall(VecRestoreArrayRead(Grad, &grad)); in DMPlexInsertBoundaryValuesRiemann()
|
| /petsc/doc/manual/ |
| H A D | ksp.md | 1294 others. For Grad-Div and Curl-Curl operators, you may want to try the
|
| H A D | tao.md | 3379 PetscCall(PetscViewerASCIIPrintf(viewer,"Grad. steps: %d\n",cg->ngradsteps));
|