1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
|
#ifndef _PCIMPL_H
#define _PCIMPL_H
#include "petscksp.h"
#include "petscpc.h"
typedef struct _PCOps *PCOps;
struct _PCOps {
PetscErrorCode (*setup)(PC);
PetscErrorCode (*apply)(PC,Vec,Vec);
PetscErrorCode (*applyrichardson)(PC,Vec,Vec,Vec,PetscReal,PetscReal,PetscReal,PetscInt,PetscTruth,PetscInt*,PCRichardsonConvergedReason*);
PetscErrorCode (*applyBA)(PC,PCSide,Vec,Vec,Vec);
PetscErrorCode (*applytranspose)(PC,Vec,Vec);
PetscErrorCode (*applyBAtranspose)(PC,PetscInt,Vec,Vec,Vec);
PetscErrorCode (*setfromoptions)(PC);
PetscErrorCode (*presolve)(PC,KSP,Vec,Vec);
PetscErrorCode (*postsolve)(PC,KSP,Vec,Vec);
PetscErrorCode (*getfactoredmatrix)(PC,Mat*);
PetscErrorCode (*applysymmetricleft)(PC,Vec,Vec);
PetscErrorCode (*applysymmetricright)(PC,Vec,Vec);
PetscErrorCode (*setuponblocks)(PC);
PetscErrorCode (*destroy)(PC);
PetscErrorCode (*view)(PC,PetscViewer);
};
/*
Preconditioner context
*/
struct _p_PC {
PETSCHEADER(struct _PCOps);
PetscInt setupcalled;
PetscInt setfromoptionscalled;
MatStructure flag;
Mat mat,pmat;
Vec diagonalscaleright,diagonalscaleleft; /* used for time integration scaling */
PetscTruth diagonalscale;
PetscTruth nonzero_guess; /* used by PCKSP, PCREDUNDANT and PCOPENMP */
PetscErrorCode (*modifysubmatrices)(PC,PetscInt,const IS[],const IS[],Mat[],void*); /* user provided routine */
void *modifysubmatricesP; /* context for user routine */
void *data;
};
extern PetscLogEvent PC_SetUp, PC_SetUpOnBlocks, PC_Apply, PC_ApplyCoarse, PC_ApplyMultiple, PC_ApplySymmetricLeft;
extern PetscLogEvent PC_ApplySymmetricRight, PC_ModifySubMatrices, PC_ApplyOnBlocks, PC_ApplyTransposeOnBlocks;
#endif
|
