Subversion Repositories slepc-dev

/*

   - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

   SLEPc - Scalable Library for Eigenvalue Problem Computations

   Copyright (c) 2002-2010, Universidad Politecnica de Valencia, Spain

   This file is part of SLEPc.

   SLEPc is free software: you can redistribute it and/or modify it under  the

   terms of version 3 of the GNU Lesser General Public License as published by

   the Free Software Foundation.

   SLEPc  is  distributed in the hope that it will be useful, but WITHOUT  ANY

   WARRANTY;  without even the implied warranty of MERCHANTABILITY or  FITNESS

   FOR  A  PARTICULAR PURPOSE. See the GNU Lesser General Public  License  for

   more details.

   You  should have received a copy of the GNU Lesser General  Public  License

   along with SLEPc. If not, see <http://www.gnu.org/licenses/>.

   - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

*/

#ifndef _SVDIMPL

#define _SVDIMPL

#include <slepcsvd.h>

#include <slepcip.h>

extern PetscFList SVDList;

extern PetscLogEvent SVD_SetUp, SVD_Solve, SVD_Dense;

typedef struct _SVDOps *SVDOps;

struct _SVDOps {

  PetscErrorCode  (*solve)(SVD);

  PetscErrorCode  (*setup)(SVD);

  PetscErrorCode  (*setfromoptions)(SVD);

  PetscErrorCode  (*publishoptions)(SVD);

  PetscErrorCode  (*destroy)(SVD);

  PetscErrorCode  (*view)(SVD,PetscViewer);

};

/*

     Maximum number of monitors you can run with a single SVD

*/

#define MAXSVDMONITORS 5 

/*

   Defines the SVD data structure.

*/

struct _p_SVD {

  PETSCHEADER(struct _SVDOps);

  Mat              OP;          /* problem matrix */

  Mat              A;           /* problem matrix (m>n) */

  Mat              AT;          /* transposed matrix */

  SVDTransposeMode transmode;   /* transpose mode */

  PetscReal        *sigma;      /* singular values */

  PetscInt         *perm;       /* permutation for singular value ordering */

  Vec              *U,*V;       /* left and right singular vectors */

  Vec              *IS;         /* placeholder for references to user-provided initial space */

  PetscInt         n;           /* maximun size of descomposition */

  SVDWhich         which;       /* which singular values are computed */

  PetscInt         nconv;       /* number of converged values */

  PetscInt         nsv;         /* number of requested values */

  PetscInt         ncv;         /* basis size */

  PetscInt         mpd;         /* maximum dimension of projected problem */

  PetscInt         nini;        /* number of initial vectors (negative means not copied yet) */

  PetscInt         its;         /* iteration counter */

  PetscInt         max_it;      /* max iterations */

  PetscReal        tol;         /* tolerance */

  PetscReal        *errest;     /* error estimates */

  PetscRandom      rand;        /* random number generator */

  void             *data;       /* placeholder for misc stuff associated

                                   with a particular solver */

  PetscInt         setupcalled;

  SVDConvergedReason reason;

  IP               ip;

  PetscBool        trackall;

  PetscErrorCode  (*monitor[MAXSVDMONITORS])(SVD,PetscInt,PetscInt,PetscReal*,PetscReal*,PetscInt,void*);

  PetscErrorCode  (*monitordestroy[MAXSVDMONITORS])(void*);

  void            *monitorcontext[MAXSVDMONITORS];

  PetscInt        numbermonitors;

  PetscInt        matvecs;

};

extern PetscErrorCode SVDRegisterAll(const char *);

extern PetscErrorCode SVDInitializePackage(const char*);

extern PetscErrorCode SVDFinalizePackage(void);

#define SVDMonitor(svd,it,nconv,sigma,errest,nest) \

        { PetscErrorCode _ierr; PetscInt _i,_im = svd->numbermonitors; \

          for ( _i=0; _i<_im; _i++ ) {\

            _ierr=(*svd->monitor[_i])(svd,it,nconv,sigma,errest,nest,svd->monitorcontext[_i]);\

            CHKERRQ(_ierr); \

          } \

        }

/* context for SVDMonitorConverged */

typedef struct {

  PetscViewerASCIIMonitor viewer;

  PetscInt oldnconv;

} SVDMONITOR_CONV;

extern PetscErrorCode SVDMonitorDestroy_Converged(SVDMONITOR_CONV*);

extern PetscErrorCode SVDDestroy_Default(SVD);

extern PetscErrorCode SVDMatMult(SVD,PetscBool,Vec,Vec);

extern PetscErrorCode SVDMatGetVecs(SVD,Vec*,Vec*);

extern PetscErrorCode SVDMatGetSize(SVD,PetscInt*,PetscInt*);

extern PetscErrorCode SVDMatGetLocalSize(SVD,PetscInt*,PetscInt*);

extern PetscErrorCode SVDTwoSideLanczos(SVD,PetscReal*,PetscReal*,Vec*,Vec,Vec*,PetscInt,PetscInt,PetscScalar*);

#endif