Subversion Repositories slepc-dev

Rev

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

Rev 1925 Rev 1940
Line 258... Line 258... 
    svd->its++;
 
    svd->its++;
 
 
 
 
 
    /* inner loop */
 
    /* inner loop */
 
    nv = PetscMin(svd->nconv+svd->mpd,svd->n);
 
    nv = PetscMin(svd->nconv+svd->mpd,svd->n);
 
    if (lanczos->oneside) {
 
    if (lanczos->oneside) {
 
      if (orthog == IP_MGS_ORTH) {
 
      if (orthog == IP_ORTH_MGS) {
 
        ierr = SVDOneSideTRLanczosMGS(svd,alpha,beta,b+svd->nconv,svd->V,v,svd->U,svd->nconv,l,nv,swork);CHKERRQ(ierr);
 
        ierr = SVDOneSideTRLanczosMGS(svd,alpha,beta,b+svd->nconv,svd->V,v,svd->U,svd->nconv,l,nv,swork);CHKERRQ(ierr);
 
      } else {
 
      } else {
 
        ierr = SVDOneSideTRLanczosCGS(svd,alpha,beta,b+svd->nconv,svd->V,v,svd->U,svd->nconv,l,nv,swork);CHKERRQ(ierr);
 
        ierr = SVDOneSideTRLanczosCGS(svd,alpha,beta,b+svd->nconv,svd->V,v,svd->U,svd->nconv,l,nv,swork);CHKERRQ(ierr);
 
      }
 
      }
 
    } else {
 
    } else {