Subversion Repositories slepc-dev

Compare Revisions

Ignore whitespace Rev 2343 → Rev 2344

/trunk/src/eps/interface/mem.c
72,9 → 72,7
ierr = PetscFree(eps->errest);CHKERRQ(ierr);
ierr = PetscFree(eps->errest_left);CHKERRQ(ierr);
ierr = SlepcVecDestroyVecs(eps->allocated_ncv,&eps->V);CHKERRQ(ierr);
if (eps->W) {
ierr = SlepcVecDestroyVecs(eps->allocated_ncv,&eps->W);CHKERRQ(ierr);
}
ierr = SlepcVecDestroyVecs(eps->allocated_ncv,&eps->W);CHKERRQ(ierr);
eps->allocated_ncv = 0;
}
PetscFunctionReturn(0);
/trunk/src/svd/interface/svdbasic.c
281,9 → 281,7
ierr = PetscFree((*svd)->sigma);CHKERRQ(ierr);
ierr = PetscFree((*svd)->perm);CHKERRQ(ierr);
ierr = PetscFree((*svd)->errest);CHKERRQ(ierr);
if ((*svd)->U) {
ierr = SlepcVecDestroyVecs((*svd)->n,&(*svd)->U);CHKERRQ(ierr);
}
ierr = SlepcVecDestroyVecs((*svd)->n,&(*svd)->U);CHKERRQ(ierr);
ierr = SlepcVecDestroyVecs((*svd)->n,&(*svd)->V);CHKERRQ(ierr);
}
ierr = SVDMonitorCancel(*svd);CHKERRQ(ierr);
/trunk/src/svd/impls/trlanczos/trlanczos.c
56,12 → 56,9
}
if (!svd->mpd) svd->mpd = svd->ncv;
if (svd->ncv>svd->nsv+svd->mpd) SETERRQ(((PetscObject)svd)->comm,1,"The value of ncv must not be larger than nev+mpd");
if (!svd->max_it)
svd->max_it = PetscMax(N/svd->ncv,100);
if (!svd->max_it) svd->max_it = PetscMax(N/svd->ncv,100);
if (svd->ncv!=svd->n) {
if (svd->U) {
ierr = SlepcVecDestroyVecs(svd->n,&svd->U);CHKERRQ(ierr);
}
ierr = SlepcVecDestroyVecs(svd->n,&svd->U);CHKERRQ(ierr);
ierr = SVDMatGetLocalSize(svd,&nloc,PETSC_NULL);CHKERRQ(ierr);
ierr = VecCreateMPIWithArray(((PetscObject)svd)->comm,nloc,PETSC_DECIDE,PETSC_NULL,&t);CHKERRQ(ierr);
ierr = SlepcVecDuplicateVecs(t,svd->ncv,&svd->U);CHKERRQ(ierr);
/trunk/src/svd/impls/cyclic/cyclic.c
173,9 → 173,7
ierr = EPSGetTolerances(cyclic->eps,&svd->tol,&svd->max_it);CHKERRQ(ierr);
 
if (svd->ncv != svd->n) {
if (svd->U) {
ierr = SlepcVecDestroyVecs(svd->n,&svd->U);CHKERRQ(ierr);
}
ierr = SlepcVecDestroyVecs(svd->n,&svd->U);CHKERRQ(ierr);
ierr = SVDMatGetLocalSize(svd,&nloc,PETSC_NULL);CHKERRQ(ierr);
ierr = VecCreateMPIWithArray(((PetscObject)svd)->comm,nloc,PETSC_DECIDE,PETSC_NULL,&t);CHKERRQ(ierr);
ierr = SlepcVecDuplicateVecs(t,svd->ncv,&svd->U);CHKERRQ(ierr);
/trunk/src/svd/impls/lanczos/gklanczos.c
57,11 → 57,8
}
if (!svd->mpd) svd->mpd = svd->ncv;
if (svd->ncv>svd->nsv+svd->mpd) SETERRQ(((PetscObject)svd)->comm,1,"The value of ncv must not be larger than nev+mpd");
if (!svd->max_it)
svd->max_it = PetscMax(N/svd->ncv,100);
if (svd->U) {
ierr = SlepcVecDestroyVecs(svd->n,&svd->U);CHKERRQ(ierr);
}
if (!svd->max_it) svd->max_it = PetscMax(N/svd->ncv,100);
ierr = SlepcVecDestroyVecs(svd->n,&svd->U);CHKERRQ(ierr);
if (!lanczos->oneside) {
ierr = SVDMatGetLocalSize(svd,&nloc,PETSC_NULL);CHKERRQ(ierr);
ierr = VecCreateMPIWithArray(((PetscObject)svd)->comm,nloc,PETSC_DECIDE,PETSC_NULL,&t);CHKERRQ(ierr);
/trunk/src/svd/impls/lapack/svdlapack.c
38,9 → 38,7
if (svd->mpd) PetscInfo(svd,"Warning: parameter mpd ignored\n");
svd->max_it = 1;
if (svd->ncv!=svd->n) {
if (svd->U) {
ierr = SlepcVecDestroyVecs(svd->n,&svd->U);CHKERRQ(ierr);
}
ierr = SlepcVecDestroyVecs(svd->n,&svd->U);CHKERRQ(ierr);
ierr = SVDMatGetLocalSize(svd,&nloc,PETSC_NULL);CHKERRQ(ierr);
ierr = VecCreateMPIWithArray(((PetscObject)svd)->comm,nloc,PETSC_DECIDE,PETSC_NULL,&t);CHKERRQ(ierr);
ierr = SlepcVecDuplicateVecs(t,svd->ncv,&svd->U);CHKERRQ(ierr);
/trunk/src/vec/contiguous.c
124,6 → 124,7
 
PetscFunctionBegin;
PetscValidPointer(V,2);
if (!*V) PetscFunctionReturn(0);
if (m <= 0) SETERRQ1(PETSC_COMM_SELF,PETSC_ERR_ARG_OUTOFRANGE,"m must be > 0: m = %D",m);
SlepcValidVecsContiguous(*V,m,2);
for (i=0;i<m;i++) {