Subversion Repositories slepc-dev

Rev

Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 1422 Rev 1504
Line 154... Line 154... 
#define __FUNCT__ "SVDSolve_CYCLIC"
 
#define __FUNCT__ "SVDSolve_CYCLIC"
 
PetscErrorCode SVDSolve_CYCLIC(SVD svd)
 
PetscErrorCode SVDSolve_CYCLIC(SVD svd)
 
{
 
{
 
  PetscErrorCode ierr;
 
  PetscErrorCode ierr;
 
  SVD_CYCLIC     *cyclic = (SVD_CYCLIC *)svd->data;
 
  SVD_CYCLIC     *cyclic = (SVD_CYCLIC *)svd->data;
 
  int            i,j;
 
  PetscInt       i,j,M,m,idx,start,end;
 
  PetscInt       M,m,idx,start,end;
 
 
 
  PetscScalar    sigma,*px;
 
  PetscScalar    sigma,*px;
 
  Vec            x;
 
  Vec            x;
 
  IS             isU,isV;
 
  IS             isU,isV;
 
  VecScatter     vsU,vsV;
 
  VecScatter     vsU,vsV;
 
 
 
 
Line 229... Line 228... 
  PetscFunctionReturn(0);
 
  PetscFunctionReturn(0);
 
}
 
}
 
 
 
 
 
#undef __FUNCT__  
 
#undef __FUNCT__  
 
#define __FUNCT__ "SVDMonitor_CYCLIC"
 
#define __FUNCT__ "SVDMonitor_CYCLIC"
 
PetscErrorCode SVDMonitor_CYCLIC(EPS eps,int its,int nconv,PetscScalar *eigr,PetscScalar *eigi,PetscReal *errest,int nest,void *ctx)
 
PetscErrorCode SVDMonitor_CYCLIC(EPS eps,PetscInt its,PetscInt nconv,PetscScalar *eigr,PetscScalar *eigi,PetscReal *errest,PetscInt nest,void *ctx)
 
{
 
{
 
  int        i,j;
 
  PetscInt   i,j;
 
  SVD        svd = (SVD)ctx;
 
  SVD        svd = (SVD)ctx;
 
 
 
 
 
  PetscFunctionBegin;
 
  PetscFunctionBegin;
 
  nconv = 0;
 
  nconv = 0;
 
  for (i=0,j=0;i<nest;i++) {
 
  for (i=0,j=0;i<nest;i++) {