| #ifndef _EPSIMPL |
| #define _EPSIMPL |
| #include "slepceps.h" |
| #include <slepceps.h> |
| extern PetscFList EPSList; |
| extern PetscLogEvent EPS_SetUp, EPS_Solve, EPS_Dense; |
| #ifndef _IPIMPL |
| #define _IPIMPL |
| #include "slepcip.h" |
| #include <slepcip.h> |
| extern PetscLogEvent IP_InnerProduct,IP_Orthogonalize,IP_ApplyMatrix; |
| #ifndef _SVDIMPL |
| #define _SVDIMPL |
| #include "slepcsvd.h" |
| #include "slepcip.h" |
| #include <slepcsvd.h> |
| #include <slepcip.h> |
| extern PetscFList SVDList; |
| extern PetscLogEvent SVD_SetUp, SVD_Solve, SVD_Dense; |
| #ifndef _QEPIMPL |
| #define _QEPIMPL |
| #include "slepcqep.h" |
| #include <slepcqep.h> |
| extern PetscFList QEPList; |
| extern PetscLogEvent QEP_SetUp, QEP_Solve, QEP_Dense; |
| #ifndef _SLEPCIMPL |
| #define _SLEPCIMPL |
| #include "slepcsys.h" |
| #include <slepcsys.h> |
| /* Private functions that are shared by several classes */ |
| #ifndef _STIMPL |
| #define _STIMPL |
| #include "slepceps.h" |
| #include <slepceps.h> |
| extern PetscLogEvent ST_SetUp, ST_Apply, ST_ApplyB, ST_ApplyTranspose; |
| extern PetscFList STList; |
| #include "petscsys.h" |
| #include "petscvec.h" |
| #include "petscmat.h" |
| /* |
| slepcimpl.h contains definitions common to all SLEPc objects |
| */ |
| #include "private/slepcimpl.h" |
| PETSC_EXTERN_CXX_BEGIN |
| /* |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/fortranimpl.h" |
| #include "slepcst.h" |
| #include <private/fortranimpl.h> |
| #include <slepcst.h> |
| #ifdef PETSC_HAVE_FORTRAN_CAPS |
| #define stsettype_ STSETTYPE |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/stimpl.h" /*I "slepcst.h" I*/ |
| #include <private/stimpl.h> /*I "slepcst.h" I*/ |
| #undef __FUNCT__ |
| #define __FUNCT__ "STApply" |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/stimpl.h" /*I "slepcst.h" I*/ |
| #include <private/stimpl.h> /*I "slepcst.h" I*/ |
| EXTERN_C_BEGIN |
| extern PetscErrorCode STCreate_Shell(ST); |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/stimpl.h" /*I "slepcst.h" I*/ |
| #include <private/stimpl.h> /*I "slepcst.h" I*/ |
| PetscClassId ST_CLASSID = 0; |
| PetscLogEvent ST_SetUp = 0, ST_Apply = 0, ST_ApplyTranspose = 0; |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/stimpl.h" /*I "slepcst.h" I*/ |
| #include "petscsys.h" |
| #include <private/stimpl.h> /*I "slepcst.h" I*/ |
| /* |
| Contains the list of registered ST routines |
| */ |
| #include "private/stimpl.h" |
| #include <private/stimpl.h> |
| #undef __FUNCT__ |
| #define __FUNCT__ "STMatShellMult" |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/stimpl.h" /*I "slepcst.h" I*/ |
| #include <private/stimpl.h> /*I "slepcst.h" I*/ |
| #undef __FUNCT__ |
| #define __FUNCT__ "STAssociatedKSPSolve" |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/stimpl.h" /*I "slepcst.h" I*/ |
| #include <private/stimpl.h> /*I "slepcst.h" I*/ |
| typedef struct { |
| Vec w2; |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/stimpl.h" /*I "slepcst.h" I*/ |
| #include <private/stimpl.h> /*I "slepcst.h" I*/ |
| typedef struct { |
| PetscScalar nu; |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/fortranimpl.h" |
| #include "slepcst.h" |
| #include <private/fortranimpl.h> |
| #include <slepcst.h> |
| #ifdef PETSC_HAVE_FORTRAN_CAPS |
| #define stshellgetcontext_ STSHELLGETCONTEXT |
| */ |
| #include "private/stimpl.h" /*I "slepcst.h" I*/ |
| #include "slepceps.h" |
| EXTERN_C_BEGIN |
| typedef struct { |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/stimpl.h" /*I "slepcst.h" I*/ |
| #include <private/stimpl.h> /*I "slepcst.h" I*/ |
| PetscErrorCode STDestroy_Precond(ST st); |
| PetscErrorCode STSetFromOptions_Precond(ST st); |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/stimpl.h" /*I "slepcst.h" I*/ |
| #include <private/stimpl.h> /*I "slepcst.h" I*/ |
| #undef __FUNCT__ |
| #define __FUNCT__ "STApply_Sinvert" |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/stimpl.h" /*I "slepcst.h" I*/ |
| #include <private/stimpl.h> /*I "slepcst.h" I*/ |
| #undef __FUNCT__ |
| #define __FUNCT__ "STApply_Shift" |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/epsimpl.h" /*I "slepceps.h" I*/ |
| #include "slepcblaslapack.h" |
| #include <private/epsimpl.h> /*I "slepceps.h" I*/ |
| #include <slepcblaslapack.h> |
| #undef __FUNCT__ |
| #define __FUNCT__ "EPSDenseNHEP" |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/fortranimpl.h" |
| #include "slepceps.h" |
| #include "private/epsimpl.h" |
| #include <private/fortranimpl.h> |
| #include <private/epsimpl.h> /*I "slepceps.h" I*/ |
| #ifdef PETSC_HAVE_FORTRAN_CAPS |
| #define epsview_ EPSVIEW |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/epsimpl.h" /*I "slepceps.h" I*/ |
| #include <private/epsimpl.h> /*I "slepceps.h" I*/ |
| #undef __FUNCT__ |
| #define __FUNCT__ "EPSAllocateSolution" |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/epsimpl.h" /*I "slepceps.h" I*/ |
| #include <private/epsimpl.h> /*I "slepceps.h" I*/ |
| #undef __FUNCT__ |
| #define __FUNCT__ "EPSSetUp" |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/epsimpl.h" /*I "slepceps.h" I*/ |
| #include <private/epsimpl.h> /*I "slepceps.h" I*/ |
| PetscFList EPSList = 0; |
| PetscClassId EPS_CLASSID = 0; |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/epsimpl.h" /*I "slepceps.h" I*/ |
| #include "slepcblaslapack.h" |
| #include <private/epsimpl.h> /*I "slepceps.h" I*/ |
| #include <slepcblaslapack.h> |
| #undef __FUNCT__ |
| #define __FUNCT__ "EPSDestroy_Default" |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/epsimpl.h" /*I "slepceps.h" I*/ |
| #include <private/epsimpl.h> /*I "slepceps.h" I*/ |
| #undef __FUNCT__ |
| #define __FUNCT__ "EPSSetFromOptions" |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/epsimpl.h" /*I "slepceps.h" I*/ |
| #include <private/epsimpl.h> /*I "slepceps.h" I*/ |
| EXTERN_C_BEGIN |
| extern PetscErrorCode EPSCreate_POWER(EPS); |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/epsimpl.h" /*I "slepceps.h" I*/ |
| #include <private/epsimpl.h> /*I "slepceps.h" I*/ |
| #undef __FUNCT__ |
| #define __FUNCT__ "EPSMonitorSet" |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/epsimpl.h" /*I "slepceps.h" I*/ |
| #include <private/epsimpl.h> /*I "slepceps.h" I*/ |
| typedef struct { |
| /* old values of eps */ |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "../src/eps/impls/external/trlan/trlanp.h" |
| #include <../src/eps/impls/external/trlan/trlanp.h> |
| PetscErrorCode EPSSolve_TRLAN(EPS); |
| #if !defined(__TRLANP_H) |
| #define __TRLANP_H |
| #include "private/epsimpl.h" |
| #include <private/epsimpl.h> /*I "slepceps.h" I*/ |
| typedef struct { |
| PetscBLASInt maxlan; |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "../src/eps/impls/external/arpack/arpackp.h" |
| #include "private/stimpl.h" |
| #include <../src/eps/impls/external/arpack/arpackp.h> |
| #include <private/stimpl.h> /*I "slepcst.h" I*/ |
| PetscErrorCode EPSSolve_ARPACK(EPS); |
| #if !defined(__ARPACKP_H) |
| #define __ARPACKP_H |
| #include "private/epsimpl.h" /*I "slepceps.h" I*/ |
| #include <private/epsimpl.h> /*I "slepceps.h" I*/ |
| typedef struct { |
| PetscBool *select; |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "../src/eps/impls/external/blzpack/blzpackp.h" /*I "slepceps.h" I*/ |
| #include "private/stimpl.h" |
| #include <../src/eps/impls/external/blzpack/blzpackp.h> |
| #include <private/stimpl.h> /*I "slepcst.h" I*/ |
| PetscErrorCode EPSSolve_BLZPACK(EPS); |
| #if !defined(__BLZPACKP_H) |
| #define __BLZPACKP_H |
| #include "private/epsimpl.h" |
| #include <private/epsimpl.h> /*I "slepceps.h" I*/ |
| typedef struct { |
| PetscBLASInt block_size; /* block size */ |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "petscsys.h" |
| #include "petscvec.h" |
| #include "petscmat.h" |
| #include <petscsys.h> |
| #include <petscvec.h> |
| #include <assert.h> |
| #include <stdlib.h> |
| #include "petscblaslapack.h" |
| #include "interpreter.h" |
| #include "temp_multivector.h" |
| #define BlopexInt PetscInt |
| #include <interpreter.h> |
| #include <temp_multivector.h> |
| #include "slepc-interface.h" |
| static void* mv_TempMultiVectorCreateFromPETScVector( void* ii_, int n, void* sample ) |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/stimpl.h" |
| #include "private/epsimpl.h" |
| #include <private/stimpl.h> /*I "slepcst.h" I*/ |
| #include <private/epsimpl.h> /*I "slepceps.h" I*/ |
| #define BlopexInt PetscInt |
| #include "slepc-interface.h" |
| #include "lobpcg.h" |
| #include "interpreter.h" |
| #include "multivector.h" |
| #include "temp_multivector.h" |
| #include <lobpcg.h> |
| #include <interpreter.h> |
| #include <multivector.h> |
| #include <temp_multivector.h> |
| PetscErrorCode EPSSolve_BLOPEX(EPS); |
| PetscFunctionBegin; |
| #ifdef PETSC_USE_COMPLEX |
| info = lobpcg_solve_complex(blopex->eigenvectors,eps,OperatorAMultiVector, |
| //#ifdef PETSC_USE_COMPLEX |
| // info = lobpcg_solve_complex(blopex->eigenvectors,eps,OperatorAMultiVector, |
| // eps->isgeneralized?eps:PETSC_NULL,eps->isgeneralized?OperatorBMultiVector:PETSC_NULL, |
| // eps,Precond_FnMultiVector,blopex->Y, |
| // blopex->blap_fn,blopex->tol,eps->max_it,0,&its, |
| // eps->eigr,PETSC_NULL,0,eps->errest,PETSC_NULL,0); |
| //#else |
| info = lobpcg_solve(blopex->eigenvectors,eps,OperatorAMultiVector, |
| eps->isgeneralized?eps:PETSC_NULL,eps->isgeneralized?OperatorBMultiVector:PETSC_NULL, |
| eps,Precond_FnMultiVector,blopex->Y, |
| blopex->blap_fn,blopex->tol,eps->max_it,0,&its, |
| (komplex*)eps->eigr,PETSC_NULL,0,eps->errest,PETSC_NULL,0); |
| #else |
| info = lobpcg_solve_double(blopex->eigenvectors,eps,OperatorAMultiVector, |
| eps->isgeneralized?eps:PETSC_NULL,eps->isgeneralized?OperatorBMultiVector:PETSC_NULL, |
| eps,Precond_FnMultiVector,blopex->Y, |
| blopex->blap_fn,blopex->tol,eps->max_it,0,&its, |
| eps->eigr,PETSC_NULL,0,eps->errest,PETSC_NULL,0); |
| #endif |
| //#endif |
| if (info>0) SETERRQ1(((PetscObject)eps)->comm,PETSC_ERR_LIB,"Error in blopex (code=%d)",info); |
| eps->its = its; |
| #ifndef SLEPC_INTERFACE_HEADER |
| #define SLEPC_INTERFACE_HEADER |
| #include "lobpcg.h" |
| #include <lobpcg.h> |
| #define BlopexInt PetscInt |
| typedef struct {double real, imag;} komplex; |
| #ifdef BLOPEX_DIR |
| #include "petsc-interface.h" |
| #else |
| #include "../src/contrib/blopex/petsc-interface/petsc-interface.h" |
| #include <../src/contrib/blopex/petsc-interface/petsc-interface.h> |
| #endif |
| extern int |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/fortranimpl.h" |
| #include "slepceps.h" |
| #include <private/fortranimpl.h> |
| #include <slepceps.h> |
| #ifdef PETSC_HAVE_FORTRAN_CAPS |
| #define epsprimmegetmethod_ EPSPRIMMEGETMETHOD |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "petsc.h" |
| #include "private/epsimpl.h" /*I "slepceps.h" I*/ |
| #include "private/stimpl.h" |
| #include <petscsys.h> |
| #include <private/epsimpl.h> /*I "slepceps.h" I*/ |
| #include <private/stimpl.h> /*I "slepcst.h" I*/ |
| PetscErrorCode EPSSolve_PRIMME(EPS); |
| EXTERN_C_BEGIN |
| #include "primme.h" |
| #include <primme.h> |
| EXTERN_C_END |
| typedef struct { |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/epsimpl.h" /*I "slepceps.h" I*/ |
| #include "slepcblaslapack.h" |
| #include <private/epsimpl.h> /*I "slepceps.h" I*/ |
| #include <slepcblaslapack.h> |
| PetscErrorCode EPSSolve_SUBSPACE(EPS); |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/epsimpl.h" /*I "slepceps.h" I*/ |
| #include "slepcblaslapack.h" |
| #include <private/epsimpl.h> /*I "slepceps.h" I*/ |
| #include <slepcblaslapack.h> |
| PetscErrorCode EPSSolve_DSITRLANCZOS(EPS); |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/epsimpl.h" /*I "slepceps.h" I*/ |
| #include "private/slepcimpl.h" |
| #include "slepcblaslapack.h" |
| #include <private/epsimpl.h> /*I "slepceps.h" I*/ |
| #include <private/slepcimpl.h> /*I "slepcsys.h" I*/ |
| #include <slepcblaslapack.h> |
| #undef __FUNCT__ |
| #define __FUNCT__ "EPSBasicArnoldi" |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/epsimpl.h" /*I "slepceps.h" I*/ |
| #include "slepcblaslapack.h" |
| #include <private/epsimpl.h> /*I "slepceps.h" I*/ |
| #include <slepcblaslapack.h> |
| PetscErrorCode EPSSolve_ARNOLDI(EPS); |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/epsimpl.h" /*I "slepceps.h" I*/ |
| #include <private/epsimpl.h> /*I "slepceps.h" I*/ |
| PetscErrorCode EPSSolve_LANCZOS(EPS); |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/epsimpl.h" /*I "slepceps.h" I*/ |
| #include "slepcblaslapack.h" |
| #include <private/epsimpl.h> /*I "slepceps.h" I*/ |
| #include <slepcblaslapack.h> |
| #undef __FUNCT__ |
| #define __FUNCT__ "ArrowTridFlip" |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/epsimpl.h" /*I "slepceps.h" I*/ |
| #include "slepcblaslapack.h" |
| #include <private/epsimpl.h> /*I "slepceps.h" I*/ |
| #include <slepcblaslapack.h> |
| #undef __FUNCT__ |
| #define __FUNCT__ "EPSTranslateHarmonic" |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/epsimpl.h" /*I "slepceps.h" I*/ |
| #include "slepcblaslapack.h" |
| #include <private/epsimpl.h> /*I "slepceps.h" I*/ |
| #include <slepcblaslapack.h> |
| PetscErrorCode EPSSolve_KRYLOVSCHUR_DEFAULT(EPS); |
| extern PetscErrorCode EPSSolve_KRYLOVSCHUR_HARMONIC(EPS); |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "slepc.h" /*I "slepc.h" I*/ |
| #include "slepceps.h" /*I "slepc.h" I*/ |
| #include "slepcblaslapack.h" |
| #include "davidson.h" |
| PetscLogEvent SLEPC_SlepcDenseMatProd = 0; |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "slepc.h" |
| #include "private/epsimpl.h" |
| #include "davidson.h" |
| #include "slepcblaslapack.h" |
| #include <slepcblaslapack.h> |
| PetscErrorCode dvd_calcpairs_proj_qz(dvdDashboard *d); |
| PetscErrorCode dvd_calcpairs_qz_start(dvdDashboard *d); |
| */ |
| #include "davidson.h" |
| #include "slepcblaslapack.h" |
| #include "slepcvec.h" |
| #include <slepcvec.h> |
| PetscErrorCode dvd_improvex_PfuncV(dvdDashboard *d, void *funcV, Vec *D, |
| PetscInt max_size_D, PetscInt r_s, PetscInt r_e, Vec *auxV, |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/epsimpl.h" |
| #include "private/stimpl.h" |
| #include "davidson.h" |
| #include "slepcblaslapack.h" |
| PetscErrorCode EPSView_DAVIDSON(EPS eps,PetscViewer viewer); |
| options. |
| */ |
| #include "petsc.h" |
| #include "private/epsimpl.h" |
| #include "private/stimpl.h" |
| #include <private/epsimpl.h> /*I "slepceps.h" I*/ |
| #include <private/stimpl.h> /*I "slepcst.h" I*/ |
| #include <slepcblaslapack.h> |
| typedef struct _dvdFunctionList { |
| PetscErrorCode (*f)(void*); |
| void *d; |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/epsimpl.h" /*I "slepceps.h" I*/ |
| #include "private/stimpl.h" |
| #include "../src/eps/impls/davidson/common/davidson.h" |
| #include "slepcblaslapack.h" |
| #include <../src/eps/impls/davidson/common/davidson.h> |
| PetscErrorCode EPSSetUp_GD(EPS eps); |
| PetscErrorCode EPSDestroy_GD(EPS eps); |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/epsimpl.h" /*I "slepceps.h" I*/ |
| #include "private/stimpl.h" |
| #include "../src/eps/impls/davidson/common/davidson.h" |
| #include "slepcblaslapack.h" |
| #include <private/epsimpl.h> /*I "slepceps.h" I*/ |
| #include <private/stimpl.h> /*I "slepcst.h" I*/ |
| #include <../src/eps/impls/davidson/common/davidson.h> |
| #include <slepcblaslapack.h> |
| PetscErrorCode EPSSetUp_JD(EPS eps); |
| PetscErrorCode EPSDestroy_JD(EPS eps); |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/epsimpl.h" |
| #include "slepcblaslapack.h" |
| #include <private/epsimpl.h> /*I "slepceps.h" I*/ |
| #include <slepcblaslapack.h> |
| typedef struct { |
| Mat OP,A,B; |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/epsimpl.h" /*I "slepceps.h" I*/ |
| #include "slepcblaslapack.h" |
| #include <private/epsimpl.h> /*I "slepceps.h" I*/ |
| #include <slepcblaslapack.h> |
| PetscErrorCode EPSSolve_POWER(EPS); |
| PetscErrorCode EPSSolve_TS_POWER(EPS); |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/fortranimpl.h" |
| #include "slepcip.h" |
| #include <private/fortranimpl.h> |
| #include <slepcip.h> |
| #ifdef PETSC_HAVE_FORTRAN_CAPS |
| #define ipcreate_ IPCREATE |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/ipimpl.h" /*I "slepcip.h" I*/ |
| #include <private/ipimpl.h> /*I "slepcip.h" I*/ |
| #undef __FUNCT__ |
| #define __FUNCT__ "IPNorm" |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/ipimpl.h" /*I "slepcip.h" I*/ |
| #include <private/ipimpl.h> /*I "slepcip.h" I*/ |
| PetscClassId IP_CLASSID = 0; |
| PetscLogEvent IP_InnerProduct = 0, IP_Orthogonalize = 0, IP_ApplyMatrix = 0; |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/ipimpl.h" /*I "slepcip.h" I*/ |
| #include "slepcblaslapack.h" |
| #include <private/ipimpl.h> /*I "slepcip.h" I*/ |
| #include <slepcblaslapack.h> |
| /* |
| IPOrthogonalizeMGS1 - Compute one step of Modified Gram-Schmidt |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/ipimpl.h" /*I "slepcip.h" I*/ |
| #include <private/ipimpl.h> /*I "slepcip.h" I*/ |
| #undef __FUNCT__ |
| #define __FUNCT__ "IPSetBilinearForm" |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "petscmat.h" |
| #include "private/vecimpl.h" /*I "vec.h" I*/ |
| #include <private/vecimpl.h> /*I "petsvec.h" I*/ |
| #ifdef __WITH_MPI__ |
| #define __SUF__(A) A##_MPI |
| " -n <n>, where <n> = number of grid subdivisions in x dimension.\n" |
| " -m <m>, where <m> = number of grid subdivisions in y dimension.\n\n"; |
| #include "slepceps.h" |
| #include <slepceps.h> |
| #undef __FUNCT__ |
| #define __FUNCT__ "main" |
| "The command line options are:\n" |
| " -n <n>, where <n> = number of grid subdivisions in both x and y dimensions.\n\n"; |
| #include "slepceps.h" |
| #include "petscblaslapack.h" |
| #include <slepceps.h> |
| #include <petscblaslapack.h> |
| /* |
| User-defined routines |
| "The command line options are:\n" |
| " -file <filename>, where <filename> = matrix file in PETSc binary form.\n\n"; |
| #include "slepceps.h" |
| #include <slepceps.h> |
| #undef __FUNCT__ |
| #define __FUNCT__ "main" |
| "The command line options are:\n" |
| " -m <m>, where <m> = number of grid subdivisions in each dimension.\n\n"; |
| #include "slepceps.h" |
| #include <slepceps.h> |
| /* |
| User-defined routines |
| " -f1 <filename>, where <filename> = matrix (A) file in PETSc binary form.\n" |
| " -f2 <filename>, where <filename> = matrix (B) file in PETSc binary form.\n\n"; |
| #include "slepceps.h" |
| #include <slepceps.h> |
| #undef __FUNCT__ |
| #define __FUNCT__ "main" |
| "The command line options are:\n" |
| " -n <n>, where <n> = matrix dimension.\n\n"; |
| #include "slepcsvd.h" |
| #include <slepcsvd.h> |
| /* |
| This example computes the singular values of an nxn Grcar matrix, |
| " -alpha <alpha>, -beta <beta>, -delta1 <delta1>, -delta2 <delta2>,\n" |
| " where <alpha> <beta> <delta1> <delta2> = model parameters.\n\n"; |
| #include "slepceps.h" |
| #include <slepceps.h> |
| /* |
| This example computes the eigenvalues with largest real part of the |
| program main |
| implicit none |
| #include "finclude/petscsys.h" |
| #include "finclude/petscvec.h" |
| #include "finclude/petscmat.h" |
| #include "finclude/slepcsys.h" |
| #include "finclude/slepceps.h" |
| #include <finclude/petscsys.h> |
| #include <finclude/petscvec.h> |
| #include <finclude/petscmat.h> |
| #include <finclude/slepcsys.h> |
| #include <finclude/slepceps.h> |
| ! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| ! Declarations |
| program main |
| implicit none |
| #include "finclude/petscsys.h" |
| #include "finclude/petscvec.h" |
| #include "finclude/petscmat.h" |
| #include "finclude/slepcsys.h" |
| #include "finclude/slepceps.h" |
| #include <finclude/petscsys.h> |
| #include <finclude/petscvec.h> |
| #include <finclude/petscmat.h> |
| #include <finclude/slepcsys.h> |
| #include <finclude/slepceps.h> |
| ! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| ! Declarations |
| subroutine MatIsing_Mult(A,x,y,ierr) |
| implicit none |
| #include "finclude/petscsys.h" |
| #include "finclude/petscvec.h" |
| #include "finclude/petscmat.h" |
| #include <finclude/petscsys.h> |
| #include <finclude/petscvec.h> |
| #include <finclude/petscmat.h> |
| Mat A |
| Vec x,y |
| "The command line options are:\n" |
| " -n <n>, where <n> = number of grid subdivisions = matrix dimension.\n\n"; |
| #include "slepceps.h" |
| #include <slepceps.h> |
| /* Define context for user-provided spectral transformation */ |
| typedef struct { |
| " -n <n>, where <n> = number of grid subdivisions in x dimension.\n" |
| " -m <m>, where <m> = number of grid subdivisions in y dimension.\n\n"; |
| #include "slepceps.h" |
| #include <slepceps.h> |
| #undef __FUNCT__ |
| #define __FUNCT__ "main" |
| ! |
| program main |
| #include "finclude/slepcqepdef.h" |
| #include <finclude/slepcqepdef.h> |
| use slepcqep |
| implicit none |
| ! For usage without modules, uncomment the following lines and remove |
| ! the previous lines between 'program main' and 'implicit none' |
| ! |
| !#include "finclude/petsc.h" |
| !#include "finclude/slepc.h" |
| !#include <finclude/petsc.h> |
| !#include <finclude/slepc.h> |
| ! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| ! Declarations |
| "The command line options are:\n" |
| " -m <m>, where <m> = number of grid subdivisions in each dimension.\n\n"; |
| #include "slepceps.h" |
| #include <slepceps.h> |
| /* |
| User-defined routines |
| " -m <m>, where <m> = number of grid subdivisions in y dimension.\n" |
| " -nulldim <k>, where <k> = dimension of the nullspace of B.\n\n"; |
| #include "slepceps.h" |
| #include <slepceps.h> |
| #undef __FUNCT__ |
| #define __FUNCT__ "main" |
| "The command line options are:\n" |
| " -file <filename>, where <filename> = matrix file in PETSc binary form.\n\n"; |
| #include "slepcsvd.h" |
| #include <slepcsvd.h> |
| #undef __FUNCT__ |
| #define __FUNCT__ "main" |
| " -n <n>, where <n> = matrix dimension.\n" |
| " -mu <mu>, where <mu> = subdiagonal value.\n\n"; |
| #include "slepcsvd.h" |
| #include <slepcsvd.h> |
| #undef __FUNCT__ |
| #define __FUNCT__ "main" |
| " -n <n>, where <n> = number of grid subdivisions in x dimension.\n" |
| " -m <m>, where <m> = number of grid subdivisions in y dimension.\n\n"; |
| #include "slepcqep.h" |
| #include <slepcqep.h> |
| #undef __FUNCT__ |
| #define __FUNCT__ "main" |
| " -C <filename>, where <filename> = matrix (C) file in PETSc binary form.\n" |
| " -K <filename>, where <filename> = matrix (K) file in PETSc binary form.\n\n"; |
| #include "slepcqep.h" |
| #include <slepcqep.h> |
| #undef __FUNCT__ |
| #define __FUNCT__ "main" |
| "The command line options are:\n" |
| " -m <m>, where <m> = number of grid subdivisions in each dimension.\n\n"; |
| #include "slepceps.h" |
| #include <slepceps.h> |
| /* |
| User-defined routines |
| ! |
| program main |
| #include "finclude/slepcepsdef.h" |
| #include <finclude/slepcepsdef.h> |
| use slepceps |
| implicit none |
| ! For usage without modules, uncomment the following lines and remove |
| ! the previous lines between 'program main' and 'implicit none' |
| ! |
| !#include "finclude/petsc.h" |
| !#include "finclude/slepc.h" |
| !#include <finclude/petsc.h> |
| !#include <finclude/slepc.h> |
| ! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| ! Declarations |
| program main |
| implicit none |
| #include "finclude/petscsys.h" |
| #include "finclude/petscvec.h" |
| #include "finclude/petscmat.h" |
| #include "finclude/slepcsys.h" |
| #include "finclude/slepcsvd.h" |
| #include <finclude/petscsys.h> |
| #include <finclude/petscvec.h> |
| #include <finclude/petscmat.h> |
| #include <finclude/slepcsys.h> |
| #include <finclude/slepcsvd.h> |
| ! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| ! Declarations |
| "The command line options are:\n" |
| " -n <n>, where <n> = number of grid subdivisions = matrix dimension.\n\n"; |
| #include "slepceps.h" |
| #include <slepceps.h> |
| #undef __FUNCT__ |
| #define __FUNCT__ "main" |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "private/fortranimpl.h" |
| #include <private/fortranimpl.h> |
| #ifdef PETSC_HAVE_FORTRAN_CAPS |
| #define slepcinitializefortran_ SLEPCINITIALIZEFORTRAN |
| */ |
| #include "private/fortranimpl.h" |
| #include "slepcsys.h" |
| #include "private/epsimpl.h" |
| #include "private/stimpl.h" |
| #include "private/svdimpl.h" |
| #include "private/ipimpl.h" |
| #include "private/qepimpl.h" |
| #include <stdlib.h> |
| extern PetscBool SlepcBeganPetsc; |
| extern PetscBool SlepcInitializeCalled; |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "slepcsys.h" /*I "slepcsys.h" I*/ |
| #include "petscblaslapack.h" |
| #include <slepcsys.h> |
| #include <petscblaslapack.h> |
| #include <stdlib.h> |
| PetscLogEvent SLEPC_UpdateVectors = 0, SLEPC_VecMAXPBY = 0; |
| - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
| */ |
| #include "slepcsys.h" /*I "slepcsys.h" I*/ |
| #include "private/epsimpl.h" |
| #include "private/stimpl.h" |
| #include "private/svdimpl.h" |
| #include "private/qepimpl.h" |
| #include "private/ipimpl.h" |
| #include <slepcsys.h> |
| #include <private/epsimpl.h> /*I "slepceps.h" I*/ |
| #include <private/stimpl.h> /*I "slepcst.h" I*/ |
| #include <private/svdimpl.h> /*I "slepcsvd.h" I*/ |
| #include <private/qepimpl.h> /*I "slepcqep.h" I*/ |
| #include <private/ipimpl.h> /*I "slepcip.h" I*/ |
| #include <stdlib.h> |
| #undef __FUNCT__ |