Subversion Repositories slepc-dev

Rev

Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 2818 Rev 2821
Line 743... Line 743... 
      tmp = perm[j]; perm[j] = perm[j+1]; perm[j+1] = tmp; j--;
 
      tmp = perm[j]; perm[j] = perm[j+1]; perm[j+1] = tmp; j--;
 
      if (j>=0) {
 
      if (j>=0) {
 
        ierr = (*comp_func)(re,0.0,eig[perm[j]],0.0,&result,comp_ctx);CHKERRQ(ierr);
 
        ierr = (*comp_func)(re,0.0,eig[perm[j]],0.0,&result,comp_ctx);CHKERRQ(ierr);
 
      }
 
      }
 
    }
 
    }
 
 
 
  }
 
 
 
  PetscFunctionReturn(0);
 
 
 
}
 
 
 
 
 
 
 
#undef __FUNCT__  
 
 
 
#define __FUNCT__ "PSCopyMatrix_Private"
 
 
 
/*
 
 
 
  PSCopyMatrix_Private - Copies the trailing block of a matrix (from
 
 
 
  rows/columns l to n).
 
 
 
*/
 
 
 
PetscErrorCode PSCopyMatrix_Private(PS ps,PSMatType dst,PSMatType src)
 
 
 
{
 
 
 
  PetscErrorCode ierr;
 
 
 
  PetscInt    j,m,off,ld;
 
 
 
  PetscScalar *S,*D;
 
 
 
 
 
 
 
  PetscFunctionBegin;
 
 
 
  ld  = ps->ld;
 
 
 
  m   = ps->n-ps->l;
 
 
 
  off = ps->l+ps->l*ld;
 
 
 
  S   = ps->mat[src];
 
 
 
  D   = ps->mat[dst];
 
 
 
  for (j=0;j<m;j++) {
 
 
 
    ierr = PetscMemcpy(D+off+j*ld,S+off+j*ld,m*sizeof(PetscScalar));CHKERRQ(ierr);
 
  }
 
  }
 
  PetscFunctionReturn(0);
 
  PetscFunctionReturn(0);
 
}
 
}
 
 
 
 
 
#undef __FUNCT__  
 
#undef __FUNCT__