Actual source code: ipimpl.h
1: /*
2: - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
3: SLEPc - Scalable Library for Eigenvalue Problem Computations
4: Copyright (c) 2002-2011, Universitat Politecnica de Valencia, Spain
6: This file is part of SLEPc.
7:
8: SLEPc is free software: you can redistribute it and/or modify it under the
9: terms of version 3 of the GNU Lesser General Public License as published by
10: the Free Software Foundation.
12: SLEPc is distributed in the hope that it will be useful, but WITHOUT ANY
13: WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
14: FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License for
15: more details.
17: You should have received a copy of the GNU Lesser General Public License
18: along with SLEPc. If not, see <http://www.gnu.org/licenses/>.
19: - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
20: */
22: #ifndef _IPIMPL
23: #define _IPIMPL
25: #include <slepcip.h>
27: extern PetscLogEvent IP_InnerProduct,IP_Orthogonalize,IP_ApplyMatrix;
29: typedef struct _IPOps *IPOps;
31: struct _IPOps {
32: PetscErrorCode (*normbegin)(IP,Vec,PetscReal*);
33: PetscErrorCode (*normend)(IP,Vec,PetscReal*);
34: PetscErrorCode (*innerproductbegin)(IP,Vec,Vec,PetscScalar*);
35: PetscErrorCode (*innerproductend)(IP,Vec,Vec,PetscScalar*);
36: PetscErrorCode (*minnerproductbegin)(IP,Vec,PetscInt,const Vec[],PetscScalar*);
37: PetscErrorCode (*minnerproductend)(IP,Vec,PetscInt,const Vec[],PetscScalar*);
38: };
40: struct _p_IP {
41: PETSCHEADER(struct _IPOps);
42: IPOrthogType orthog_type; /* which orthogonalization to use */
43: IPOrthogRefineType orthog_ref; /* refinement method */
44: PetscReal orthog_eta; /* refinement threshold */
45: Mat matrix;
46: PetscInt innerproducts;
48: /*------------------------- Cache Bx product -------------------*/
49: PetscInt xid;
50: PetscInt xstate;
51: Vec Bx;
52: };
54: extern PetscErrorCode IPSetDefaultType_Private(IP);
55: extern PetscErrorCode IPApplyMatrix_Private(IP,Vec);
56: extern PetscErrorCode IPOrthogonalizeCGS1(IP,PetscInt,Vec*,PetscInt,PetscBool*,Vec*,Vec,PetscScalar*,PetscReal*,PetscReal*);
58: #endif