/*
|
/*
|
SVD routines for setting solver options.
|
SVD routines for setting solver options.
|
|
|
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
SLEPc - Scalable Library for Eigenvalue Problem Computations
|
SLEPc - Scalable Library for Eigenvalue Problem Computations
|
Copyright (c) 2002-2009, Universidad Politecnica de Valencia, Spain
|
Copyright (c) 2002-2010, Universidad Politecnica de Valencia, Spain
|
|
|
This file is part of SLEPc.
|
This file is part of SLEPc.
|
|
|
SLEPc is free software: you can redistribute it and/or modify it under the
|
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
|
terms of version 3 of the GNU Lesser General Public License as published by
|
the Free Software Foundation.
|
the Free Software Foundation.
|
|
|
SLEPc is distributed in the hope that it will be useful, but WITHOUT ANY
|
SLEPc is distributed in the hope that it will be useful, but WITHOUT ANY
|
WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
|
WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
|
FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License for
|
FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License for
|
more details.
|
more details.
|
|
|
You should have received a copy of the GNU Lesser General Public License
|
You should have received a copy of the GNU Lesser General Public License
|
along with SLEPc. If not, see <http://www.gnu.org/licenses/>.
|
along with SLEPc. If not, see <http://www.gnu.org/licenses/>.
|
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
*/
|
*/
|
|
|
#include "private/svdimpl.h" /*I "slepcsvd.h" I*/
|
#include "private/svdimpl.h" /*I "slepcsvd.h" I*/
|
|
|
#undef __FUNCT__
|
#undef __FUNCT__
|
#define __FUNCT__ "SVDSetTransposeMode"
|
#define __FUNCT__ "SVDSetTransposeMode"
|
/*@
|
/*@
|
SVDSetTransposeMode - Sets how to handle the transpose of the matrix
|
SVDSetTransposeMode - Sets how to handle the transpose of the matrix
|
associated with the singular value problem.
|
associated with the singular value problem.
|
|
|
Collective on SVD
|
Collective on SVD
|
|
|
Input Parameters:
|
Input Parameters:
|
+ svd - the singular value solver context
|
+ svd - the singular value solver context
|
- mode - how to compute the transpose, one of SVD_TRANSPOSE_EXPLICIT
|
- mode - how to compute the transpose, one of SVD_TRANSPOSE_EXPLICIT
|
or SVD_TRANSPOSE_IMPLICIT (see notes below)
|
or SVD_TRANSPOSE_IMPLICIT (see notes below)
|
|
|
Options Database Key:
|
Options Database Key:
|
. -svd_transpose_mode <mode> - Indicates the mode flag, where <mode>
|
. -svd_transpose_mode <mode> - Indicates the mode flag, where <mode>
|
is one of 'explicit' or 'implicit'.
|
is one of 'explicit' or 'implicit'.
|
|
|
Notes:
|
Notes:
|
In the SVD_TRANSPOSE_EXPLICIT mode, the transpose of the matrix is
|
In the SVD_TRANSPOSE_EXPLICIT mode, the transpose of the matrix is
|
explicitly built.
|
explicitly built.
|
|
|
The option SVD_TRANSPOSE_IMPLICIT does not build the transpose, but
|
The option SVD_TRANSPOSE_IMPLICIT does not build the transpose, but
|
handles it implicitly via MatMultTranspose() operations. This is
|
handles it implicitly via MatMultTranspose() operations. This is
|
likely to be more inefficient than SVD_TRANSPOSE_EXPLICIT, both in
|
likely to be more inefficient than SVD_TRANSPOSE_EXPLICIT, both in
|
sequential and in parallel, but requires less storage.
|
sequential and in parallel, but requires less storage.
|
|
|
The default is SVD_TRANSPOSE_EXPLICIT if the matrix has defined the
|
The default is SVD_TRANSPOSE_EXPLICIT if the matrix has defined the
|
MatTranspose operation, and SVD_TRANSPOSE_IMPLICIT otherwise.
|
MatTranspose operation, and SVD_TRANSPOSE_IMPLICIT otherwise.
|
|
|
Level: advanced
|
Level: advanced
|
|
|
.seealso: SVDGetTransposeMode(), SVDSolve(), SVDSetOperator(),
|
.seealso: SVDGetTransposeMode(), SVDSolve(), SVDSetOperator(),
|
SVDGetOperator(), SVDTransposeMode
|
SVDGetOperator(), SVDTransposeMode
|
@*/
|
@*/
|
PetscErrorCode SVDSetTransposeMode(SVD svd,SVDTransposeMode mode)
|
PetscErrorCode SVDSetTransposeMode(SVD svd,SVDTransposeMode mode)
|
{
|
{
|
PetscFunctionBegin;
|
PetscFunctionBegin;
|
PetscValidHeaderSpecific(svd,SVD_COOKIE,1);
|
PetscValidHeaderSpecific(svd,SVD_COOKIE,1);
|
if (mode == PETSC_DEFAULT || mode == PETSC_DECIDE) mode = (SVDTransposeMode)PETSC_DECIDE;
|
if (mode == PETSC_DEFAULT || mode == PETSC_DECIDE) mode = (SVDTransposeMode)PETSC_DECIDE;
|
else switch (mode) {
|
else switch (mode) {
|
case SVD_TRANSPOSE_EXPLICIT:
|
case SVD_TRANSPOSE_EXPLICIT:
|
case SVD_TRANSPOSE_IMPLICIT:
|
case SVD_TRANSPOSE_IMPLICIT:
|
svd->transmode = mode;
|
svd->transmode = mode;
|
svd->setupcalled = 0;
|
svd->setupcalled = 0;
|
break;
|
break;
|
default:
|
default:
|
SETERRQ(PETSC_ERR_ARG_OUTOFRANGE,"Invalid transpose mode");
|
SETERRQ(PETSC_ERR_ARG_OUTOFRANGE,"Invalid transpose mode");
|
}
|
}
|
PetscFunctionReturn(0);
|
PetscFunctionReturn(0);
|
}
|
}
|
|
|
#undef __FUNCT__
|
#undef __FUNCT__
|
#define __FUNCT__ "SVDGetTransposeMode"
|
#define __FUNCT__ "SVDGetTransposeMode"
|
/*@C
|
/*@C
|
SVDGetTransposeMode - Gets the mode used to compute the transpose
|
SVDGetTransposeMode - Gets the mode used to compute the transpose
|
of the matrix associated with the singular value problem.
|
of the matrix associated with the singular value problem.
|
|
|
Not collective
|
Not collective
|
|
|
Input Parameter:
|
Input Parameter:
|
. svd - the singular value solver context
|
. svd - the singular value solver context
|
|
|
Output paramter:
|
Output paramter:
|
. mode - how to compute the transpose, one of SVD_TRANSPOSE_EXPLICIT
|
. mode - how to compute the transpose, one of SVD_TRANSPOSE_EXPLICIT
|
or SVD_TRANSPOSE_IMPLICIT
|
or SVD_TRANSPOSE_IMPLICIT
|
|
|
Level: advanced
|
Level: advanced
|
|
|
.seealso: SVDSetTransposeMode(), SVDSolve(), SVDSetOperator(),
|
.seealso: SVDSetTransposeMode(), SVDSolve(), SVDSetOperator(),
|
SVDGetOperator(), SVDTransposeMode
|
SVDGetOperator(), SVDTransposeMode
|
@*/
|
@*/
|
PetscErrorCode SVDGetTransposeMode(SVD svd,SVDTransposeMode *mode)
|
PetscErrorCode SVDGetTransposeMode(SVD svd,SVDTransposeMode *mode)
|
{
|
{
|
PetscFunctionBegin;
|
PetscFunctionBegin;
|
PetscValidHeaderSpecific(svd,SVD_COOKIE,1);
|
PetscValidHeaderSpecific(svd,SVD_COOKIE,1);
|
PetscValidPointer(mode,2);
|
PetscValidPointer(mode,2);
|
*mode = svd->transmode;
|
*mode = svd->transmode;
|
PetscFunctionReturn(0);
|
PetscFunctionReturn(0);
|
}
|
}
|
|
|
#undef __FUNCT__
|
#undef __FUNCT__
|
#define __FUNCT__ "SVDSetTolerances"
|
#define __FUNCT__ "SVDSetTolerances"
|
/*@
|
/*@
|
SVDSetTolerances - Sets the tolerance and maximum
|
SVDSetTolerances - Sets the tolerance and maximum
|
iteration count used by the default SVD convergence testers.
|
iteration count used by the default SVD convergence testers.
|
|
|
Collective on SVD
|
Collective on SVD
|
|
|
Input Parameters:
|
Input Parameters:
|
+ svd - the singular value solver context
|
+ svd - the singular value solver context
|
. tol - the convergence tolerance
|
. tol - the convergence tolerance
|
- maxits - maximum number of iterations to use
|
- maxits - maximum number of iterations to use
|
|
|
Options Database Keys:
|
Options Database Keys:
|
+ -svd_tol <tol> - Sets the convergence tolerance
|
+ -svd_tol <tol> - Sets the convergence tolerance
|
- -svd_max_it <maxits> - Sets the maximum number of iterations allowed
|
- -svd_max_it <maxits> - Sets the maximum number of iterations allowed
|
(use PETSC_DECIDE to compute an educated guess based on basis and matrix sizes)
|
(use PETSC_DECIDE to compute an educated guess based on basis and matrix sizes)
|
|
|
Notes:
|
Notes:
|
Use PETSC_IGNORE to retain the previous value of any parameter.
|
Use PETSC_IGNORE to retain the previous value of any parameter.
|
|
|
Level: intermediate
|
Level: intermediate
|
|
|
.seealso: SVDGetTolerances()
|
.seealso: SVDGetTolerances()
|
@*/
|
@*/
|
PetscErrorCode SVDSetTolerances(SVD svd,PetscReal tol,PetscInt maxits)
|
PetscErrorCode SVDSetTolerances(SVD svd,PetscReal tol,PetscInt maxits)
|
{
|
{
|
PetscFunctionBegin;
|
PetscFunctionBegin;
|
PetscValidHeaderSpecific(svd,SVD_COOKIE,1);
|
PetscValidHeaderSpecific(svd,SVD_COOKIE,1);
|
if (tol != PETSC_IGNORE) {
|
if (tol != PETSC_IGNORE) {
|
if (tol == PETSC_DEFAULT) {
|
if (tol == PETSC_DEFAULT) {
|
tol = 1e-7;
|
tol = 1e-7;
|
} else {
|
} else {
|
if (tol < 0.0) SETERRQ(PETSC_ERR_ARG_OUTOFRANGE,"Illegal value of tol. Must be > 0");
|
if (tol < 0.0) SETERRQ(PETSC_ERR_ARG_OUTOFRANGE,"Illegal value of tol. Must be > 0");
|
svd->tol = tol;
|
svd->tol = tol;
|
}
|
}
|
}
|
}
|
if (maxits != PETSC_IGNORE) {
|
if (maxits != PETSC_IGNORE) {
|
if (maxits == PETSC_DEFAULT || maxits == PETSC_DECIDE) {
|
if (maxits == PETSC_DEFAULT || maxits == PETSC_DECIDE) {
|
svd->max_it = 0;
|
svd->max_it = 0;
|
svd->setupcalled = 0;
|
svd->setupcalled = 0;
|
} else {
|
} else {
|
if (maxits < 0) SETERRQ(PETSC_ERR_ARG_OUTOFRANGE,"Illegal value of maxits. Must be > 0");
|
if (maxits < 0) SETERRQ(PETSC_ERR_ARG_OUTOFRANGE,"Illegal value of maxits. Must be > 0");
|
svd->max_it = maxits;
|
svd->max_it = maxits;
|
}
|
}
|
}
|
}
|
PetscFunctionReturn(0);
|
PetscFunctionReturn(0);
|
}
|
}
|
|
|
#undef __FUNCT__
|
#undef __FUNCT__
|
#define __FUNCT__ "SVDGetTolerances"
|
#define __FUNCT__ "SVDGetTolerances"
|
/*@
|
/*@
|
SVDGetTolerances - Gets the tolerance and maximum
|
SVDGetTolerances - Gets the tolerance and maximum
|
iteration count used by the default SVD convergence tests.
|
iteration count used by the default SVD convergence tests.
|
|
|
Not Collective
|
Not Collective
|
|
|
Input Parameter:
|
Input Parameter:
|
. svd - the singular value solver context
|
. svd - the singular value solver context
|
|
|
Output Parameters:
|
Output Parameters:
|
+ tol - the convergence tolerance
|
+ tol - the convergence tolerance
|
- maxits - maximum number of iterations
|
- maxits - maximum number of iterations
|
|
|
Notes:
|
Notes:
|
The user can specify PETSC_NULL for any parameter that is not needed.
|
The user can specify PETSC_NULL for any parameter that is not needed.
|
|
|
Level: intermediate
|
Level: intermediate
|
|
|
.seealso: SVDSetTolerances()
|
.seealso: SVDSetTolerances()
|
@*/
|
@*/
|
PetscErrorCode SVDGetTolerances(SVD svd,PetscReal *tol,PetscInt *maxits)
|
PetscErrorCode SVDGetTolerances(SVD svd,PetscReal *tol,PetscInt *maxits)
|
{
|
{
|
PetscFunctionBegin;
|
PetscFunctionBegin;
|
PetscValidHeaderSpecific(svd,SVD_COOKIE,1);
|
PetscValidHeaderSpecific(svd,SVD_COOKIE,1);
|
if (tol) *tol = svd->tol;
|
if (tol) *tol = svd->tol;
|
if (maxits) *maxits = svd->max_it;
|
if (maxits) *maxits = svd->max_it;
|
PetscFunctionReturn(0);
|
PetscFunctionReturn(0);
|
}
|
}
|
|
|
#undef __FUNCT__
|
#undef __FUNCT__
|
#define __FUNCT__ "SVDSetDimensions"
|
#define __FUNCT__ "SVDSetDimensions"
|
/*@
|
/*@
|
SVDSetDimensions - Sets the number of singular values to compute
|
SVDSetDimensions - Sets the number of singular values to compute
|
and the dimension of the subspace.
|
and the dimension of the subspace.
|
|
|
Collective on SVD
|
Collective on SVD
|
|
|
Input Parameters:
|
Input Parameters:
|
+ svd - the singular value solver context
|
+ svd - the singular value solver context
|
. nsv - number of singular values to compute
|
. nsv - number of singular values to compute
|
. ncv - the maximum dimension of the subspace to be used by the solver
|
. ncv - the maximum dimension of the subspace to be used by the solver
|
- mpd - the maximum dimension allowed for the projected problem
|
- mpd - the maximum dimension allowed for the projected problem
|
|
|
Options Database Keys:
|
Options Database Keys:
|
+ -svd_nsv <nsv> - Sets the number of singular values
|
+ -svd_nsv <nsv> - Sets the number of singular values
|
. -svd_ncv <ncv> - Sets the dimension of the subspace
|
. -svd_ncv <ncv> - Sets the dimension of the subspace
|
- -svd_mpd <mpd> - Sets the maximum projected dimension
|
- -svd_mpd <mpd> - Sets the maximum projected dimension
|
|
|
Notes:
|
Notes:
|
Use PETSC_IGNORE to retain the previous value of any parameter.
|
Use PETSC_IGNORE to retain the previous value of any parameter.
|
|
|
Use PETSC_DECIDE for ncv and mpd to assign a reasonably good value, which is
|
Use PETSC_DECIDE for ncv and mpd to assign a reasonably good value, which is
|
dependent on the solution method and the number of singular values required.
|
dependent on the solution method and the number of singular values required.
|
|
|
The parameters ncv and mpd are intimately related, so that the user is advised
|
The parameters ncv and mpd are intimately related, so that the user is advised
|
to set one of them at most. Normal usage is the following:
|
to set one of them at most. Normal usage is the following:
|
+ - In cases where nsv is small, the user sets ncv (a reasonable default is 2*nsv).
|
+ - In cases where nsv is small, the user sets ncv (a reasonable default is 2*nsv).
|
- - In cases where nsv is large, the user sets mpd.
|
- - In cases where nsv is large, the user sets mpd.
|
|
|
The value of ncv should always be between nsv and (nsv+mpd), typically
|
The value of ncv should always be between nsv and (nsv+mpd), typically
|
ncv=nsv+mpd. If nev is not too large, mpd=nsv is a reasonable choice, otherwise
|
ncv=nsv+mpd. If nev is not too large, mpd=nsv is a reasonable choice, otherwise
|
a smaller value should be used.
|
a smaller value should be used.
|
|
|
Level: intermediate
|
Level: intermediate
|
|
|
.seealso: SVDGetDimensions()
|
.seealso: SVDGetDimensions()
|
@*/
|
@*/
|
PetscErrorCode SVDSetDimensions(SVD svd,PetscInt nsv,PetscInt ncv,PetscInt mpd)
|
PetscErrorCode SVDSetDimensions(SVD svd,PetscInt nsv,PetscInt ncv,PetscInt mpd)
|
{
|
{
|
PetscFunctionBegin;
|
PetscFunctionBegin;
|
PetscValidHeaderSpecific(svd,SVD_COOKIE,1);
|
PetscValidHeaderSpecific(svd,SVD_COOKIE,1);
|
|
|
if (nsv != PETSC_IGNORE) {
|
if (nsv != PETSC_IGNORE) {
|
if (nsv<1) SETERRQ(PETSC_ERR_ARG_OUTOFRANGE,"Illegal value of nsv. Must be > 0");
|
if (nsv<1) SETERRQ(PETSC_ERR_ARG_OUTOFRANGE,"Illegal value of nsv. Must be > 0");
|
svd->nsv = nsv;
|
svd->nsv = nsv;
|
}
|
}
|
if (ncv != PETSC_IGNORE) {
|
if (ncv != PETSC_IGNORE) {
|
if (ncv == PETSC_DEFAULT || ncv == PETSC_DECIDE) {
|
if (ncv == PETSC_DEFAULT || ncv == PETSC_DECIDE) {
|
svd->ncv = 0;
|
svd->ncv = 0;
|
} else {
|
} else {
|
if (ncv<1) SETERRQ(PETSC_ERR_ARG_OUTOFRANGE,"Illegal value of ncv. Must be > 0");
|
if (ncv<1) SETERRQ(PETSC_ERR_ARG_OUTOFRANGE,"Illegal value of ncv. Must be > 0");
|
svd->ncv = ncv;
|
svd->ncv = ncv;
|
}
|
}
|
svd->setupcalled = 0;
|
svd->setupcalled = 0;
|
}
|
}
|
if( mpd != PETSC_IGNORE ) {
|
if( mpd != PETSC_IGNORE ) {
|
if (mpd == PETSC_DECIDE || mpd == PETSC_DEFAULT) {
|
if (mpd == PETSC_DECIDE || mpd == PETSC_DEFAULT) {
|
svd->mpd = 0;
|
svd->mpd = 0;
|
} else {
|
} else {
|
if (mpd<1) SETERRQ(PETSC_ERR_ARG_OUTOFRANGE,"Illegal value of mpd. Must be > 0");
|
if (mpd<1) SETERRQ(PETSC_ERR_ARG_OUTOFRANGE,"Illegal value of mpd. Must be > 0");
|
svd->mpd = mpd;
|
svd->mpd = mpd;
|
}
|
}
|
}
|
}
|
PetscFunctionReturn(0);
|
PetscFunctionReturn(0);
|
/* context for QEPMonitorConverged */
|
/* context for QEPMonitorConverged */
|
typedef struct {
|
typedef struct {
|
PetscViewerASCIIMonitor viewer;
|
PetscViewerASCIIMonitor viewer;
|
PetscInt oldnconv;
|
PetscInt oldnconv;
|
} QEPMONITOR_CONV;
|
} QEPMONITOR_CONV;
|
|
|
}
|
}
|
|
|
#undef __FUNCT__
|
#undef __FUNCT__
|
#define __FUNCT__ "SVDGetDimensions"
|
#define __FUNCT__ "SVDGetDimensions"
|
/*@
|
/*@
|
SVDGetDimensions - Gets the number of singular values to compute
|
SVDGetDimensions - Gets the number of singular values to compute
|
and the dimension of the subspace.
|
and the dimension of the subspace.
|
|
|
Not Collective
|
Not Collective
|
|
|
Input Parameter:
|
Input Parameter:
|
. svd - the singular value context
|
. svd - the singular value context
|
|
|
Output Parameters:
|
Output Parameters:
|
+ nsv - number of singular values to compute
|
+ nsv - number of singular values to compute
|
. ncv - the maximum dimension of the subspace to be used by the solver
|
. ncv - the maximum dimension of the subspace to be used by the solver
|
- mpd - the maximum dimension allowed for the projected problem
|
- mpd - the maximum dimension allowed for the projected problem
|
|
|
Notes:
|
Notes:
|
The user can specify PETSC_NULL for any parameter that is not needed.
|
The user can specify PETSC_NULL for any parameter that is not needed.
|
|
|
Level: intermediate
|
Level: intermediate
|
|
|
.seealso: SVDSetDimensions()
|
.seealso: SVDSetDimensions()
|
@*/
|
@*/
|
PetscErrorCode SVDGetDimensions(SVD svd,PetscInt *nsv,PetscInt *ncv,PetscInt *mpd)
|
PetscErrorCode SVDGetDimensions(SVD svd,PetscInt *nsv,PetscInt *ncv,PetscInt *mpd)
|
{
|
{
|
PetscFunctionBegin;
|
PetscFunctionBegin;
|
PetscValidHeaderSpecific(svd,SVD_COOKIE,1);
|
PetscValidHeaderSpecific(svd,SVD_COOKIE,1);
|
if (nsv) *nsv = svd->nsv;
|
if (nsv) *nsv = svd->nsv;
|
if (ncv) *ncv = svd->ncv;
|
if (ncv) *ncv = svd->ncv;
|
if (mpd) *mpd = svd->mpd;
|
if (mpd) *mpd = svd->mpd;
|
PetscFunctionReturn(0);
|
PetscFunctionReturn(0);
|
}
|
}
|
|
|
#undef __FUNCT__
|
#undef __FUNCT__
|
#define __FUNCT__ "SVDSetWhichSingularTriplets"
|
#define __FUNCT__ "SVDSetWhichSingularTriplets"
|
/*@
|
/*@
|
SVDSetWhichSingularTriplets - Specifies which singular triplets are
|
SVDSetWhichSingularTriplets - Specifies which singular triplets are
|
to be sought.
|
to be sought.
|
|
|
Collective on SVD
|
Collective on SVD
|
|
|
Input Parameter:
|
Input Parameter:
|
. svd - singular value solver context obtained from SVDCreate()
|
. svd - singular value solver context obtained from SVDCreate()
|
|
|
Output Parameter:
|
Output Parameter:
|
. which - which singular triplets are to be sought
|
. which - which singular triplets are to be sought
|
|
|
Possible values:
|
Possible values:
|
The parameter 'which' can have one of these values:
|
The parameter 'which' can have one of these values:
|
|
|
+ SVD_LARGEST - largest singular values
|
+ SVD_LARGEST - largest singular values
|
- SVD_SMALLEST - smallest singular values
|
- SVD_SMALLEST - smallest singular values
|
|
|
Options Database Keys:
|
Options Database Keys:
|
+ -svd_largest - Sets largest singular values
|
+ -svd_largest - Sets largest singular values
|
- -svd_smallest - Sets smallest singular values
|
- -svd_smallest - Sets smallest singular values
|
|
|
Level: intermediate
|
Level: intermediate
|
|
|
.seealso: SVDGetWhichSingularTriplets(), SVDWhich
|
.seealso: SVDGetWhichSingularTriplets(), SVDWhich
|
@*/
|
@*/
|
PetscErrorCode SVDSetWhichSingularTriplets(SVD svd,SVDWhich which)
|
PetscErrorCode SVDSetWhichSingularTriplets(SVD svd,SVDWhich which)
|
{
|
{
|
PetscFunctionBegin;
|
PetscFunctionBegin;
|
PetscValidHeaderSpecific(svd,SVD_COOKIE,1);
|
PetscValidHeaderSpecific(svd,SVD_COOKIE,1);
|
switch (which) {
|
switch (which) {
|
case SVD_LARGEST:
|
case SVD_LARGEST:
|
case SVD_SMALLEST:
|
case SVD_SMALLEST:
|
if (svd->which != which) {
|
if (svd->which != which) {
|
svd->setupcalled = 0;
|
svd->setupcalled = 0;
|
svd->which = which;
|
svd->which = which;
|
}
|
}
|
break;
|
break;
|
default:
|
default:
|
SETERRQ(PETSC_ERR_ARG_OUTOFRANGE,"Invalid 'which' parameter");
|
SETERRQ(PETSC_ERR_ARG_OUTOFRANGE,"Invalid 'which' parameter");
|
}
|
}
|
PetscFunctionReturn(0);
|
PetscFunctionReturn(0);
|
}
|
}
|
|
|
#undef __FUNCT__
|
#undef __FUNCT__
|
#define __FUNCT__ "SVDGetWhichSingularTriplets"
|
#define __FUNCT__ "SVDGetWhichSingularTriplets"
|
/*@C
|
/*@C
|
SVDGetWhichSingularTriplets - Returns which singular triplets are
|
SVDGetWhichSingularTriplets - Returns which singular triplets are
|
to be sought.
|
to be sought.
|
|
|
Not Collective
|
Not Collective
|
|
|
Input Parameter:
|
Input Parameter:
|
. svd - singular value solver context obtained from SVDCreate()
|
. svd - singular value solver context obtained from SVDCreate()
|
|
|
Output Parameter:
|
Output Parameter:
|
. which - which singular triplets are to be sought
|
. which - which singular triplets are to be sought
|
|
|
Notes:
|
Notes:
|
See SVDSetWhichSingularTriplets() for possible values of which
|
See SVDSetWhichSingularTriplets() for possible values of which
|
|
|
Level: intermediate
|
Level: intermediate
|
|
|
.seealso: SVDSetWhichSingularTriplets(), SVDWhich
|
.seealso: SVDSetWhichSingularTriplets(), SVDWhich
|
@*/
|
@*/
|
PetscErrorCode SVDGetWhichSingularTriplets(SVD svd,SVDWhich *which)
|
PetscErrorCode SVDGetWhichSingularTriplets(SVD svd,SVDWhich *which)
|
{
|
{
|
PetscFunctionBegin;
|
PetscFunctionBegin;
|
PetscValidHeaderSpecific(svd,SVD_COOKIE,1);
|
PetscValidHeaderSpecific(svd,SVD_COOKIE,1);
|
PetscValidPointer(which,2);
|
PetscValidPointer(which,2);
|
*which = svd->which;
|
*which = svd->which;
|
PetscFunctionReturn(0);
|
PetscFunctionReturn(0);
|
}
|
}
|
|
|
#undef __FUNCT__
|
#undef __FUNCT__
|
#define __FUNCT__ "SVDSetFromOptions"
|
#define __FUNCT__ "SVDSetFromOptions"
|
/*@
|
/*@
|
SVDSetFromOptions - Sets SVD options from the options database.
|
SVDSetFromOptions - Sets SVD options from the options database.
|
This routine must be called before SVDSetUp() if the user is to be
|
This routine must be called before SVDSetUp() if the user is to be
|
allowed to set the solver type.
|
allowed to set the solver type.
|
|
|
Collective on SVD
|
Collective on SVD
|
|
|
Input Parameters:
|
Input Parameters:
|
. svd - the singular value solver context
|
. svd - the singular value solver context
|
|
|
Notes:
|
Notes:
|
To see all options, run your program with the -help option.
|
To see all options, run your program with the -help option.
|
|
|
Level: beginner
|
Level: beginner
|
|
|
.seealso:
|
.seealso:
|
@*/
|
@*/
|
PetscErrorCode SVDSetFromOptions(SVD svd)
|
PetscErrorCode SVDSetFromOptions(SVD svd)
|
{
|
{
|
PetscErrorCode ierr;
|
PetscErrorCode ierr;
|
char type[256],monfilename[PETSC_MAX_PATH_LEN];
|
char type[256],monfilename[PETSC_MAX_PATH_LEN];
|
PetscTruth flg;
|
PetscTruth flg;
|
const char *mode_list[2] = { "explicit", "implicit" };
|
const char *mode_list[2] = { "explicit", "implicit" };
|
PetscInt i,j,k;
|
PetscInt i,j,k;
|
PetscReal r;
|
PetscReal r;
|
PetscViewerASCIIMonitor monviewer;
|
PetscViewerASCIIMonitor monviewer;
|
SVDMONITOR_CONV *ctx;
|
SVDMONITOR_CONV *ctx;
|
|
|
PetscFunctionBegin;
|
PetscFunctionBegin;
|
PetscValidHeaderSpecific(svd,SVD_COOKIE,1);
|
PetscValidHeaderSpecific(svd,SVD_COOKIE,1);
|
svd->setupcalled = 0;
|
svd->setupcalled = 0;
|
ierr = PetscOptionsBegin(((PetscObject)svd)->comm,((PetscObject)svd)->prefix,"Singular Value Solver (SVD) Options","SVD");CHKERRQ(ierr);
|
ierr = PetscOptionsBegin(((PetscObject)svd)->comm,((PetscObject)svd)->prefix,"Singular Value Solver (SVD) Options","SVD");CHKERRQ(ierr);
|
|
|
ierr = PetscOptionsList("-svd_type","Singular Value Solver method","SVDSetType",SVDList,(char*)(((PetscObject)svd)->type_name?((PetscObject)svd)->type_name:SVDCROSS),type,256,&flg);CHKERRQ(ierr);
|
ierr = PetscOptionsList("-svd_type","Singular Value Solver method","SVDSetType",SVDList,(char*)(((PetscObject)svd)->type_name?((PetscObject)svd)->type_name:SVDCROSS),type,256,&flg);CHKERRQ(ierr);
|
if (flg) {
|
if (flg) {
|
ierr = SVDSetType(svd,type);CHKERRQ(ierr);
|
ierr = SVDSetType(svd,type);CHKERRQ(ierr);
|
} else if (!((PetscObject)svd)->type_name) {
|
} else if (!((PetscObject)svd)->type_name) {
|
ierr = SVDSetType(svd,SVDCROSS);CHKERRQ(ierr);
|
ierr = SVDSetType(svd,SVDCROSS);CHKERRQ(ierr);
|
}
|
}
|
|
|
ierr = PetscOptionsName("-svd_view","Print detailed information on solver used","SVDView",0);CHKERRQ(ierr);
|
ierr = PetscOptionsName("-svd_view","Print detailed information on solver used","SVDView",0);CHKERRQ(ierr);
|
|
|
ierr = PetscOptionsEList("-svd_transpose_mode","Transpose SVD mode","SVDSetTransposeMode",mode_list,2,svd->transmode == PETSC_DECIDE ? "decide" : mode_list[svd->transmode],&i,&flg);CHKERRQ(ierr);
|
ierr = PetscOptionsEList("-svd_transpose_mode","Transpose SVD mode","SVDSetTransposeMode",mode_list,2,svd->transmode == PETSC_DECIDE ? "decide" : mode_list[svd->transmode],&i,&flg);CHKERRQ(ierr);
|
if (flg) {
|
if (flg) {
|
ierr = SVDSetTransposeMode(svd,(SVDTransposeMode)i);CHKERRQ(ierr);
|
ierr = SVDSetTransposeMode(svd,(SVDTransposeMode)i);CHKERRQ(ierr);
|
}
|
}
|
|
|
r = i = PETSC_IGNORE;
|
r = i = PETSC_IGNORE;
|
ierr = PetscOptionsInt("-svd_max_it","Maximum number of iterations","SVDSetTolerances",svd->max_it,&i,PETSC_NULL);CHKERRQ(ierr);
|
ierr = PetscOptionsInt("-svd_max_it","Maximum number of iterations","SVDSetTolerances",svd->max_it,&i,PETSC_NULL);CHKERRQ(ierr);
|
ierr = PetscOptionsReal("-svd_tol","Tolerance","SVDSetTolerances",svd->tol,&r,PETSC_NULL);CHKERRQ(ierr);
|
ierr = PetscOptionsReal("-svd_tol","Tolerance","SVDSetTolerances",svd->tol,&r,PETSC_NULL);CHKERRQ(ierr);
|
ierr = SVDSetTolerances(svd,r,i);CHKERRQ(ierr);
|
ierr = SVDSetTolerances(svd,r,i);CHKERRQ(ierr);
|
|
|
i = j = k = PETSC_IGNORE;
|
i = j = k = PETSC_IGNORE;
|
ierr = PetscOptionsInt("-svd_nsv","Number of singular values to compute","SVDSetDimensions",svd->nsv,&i,PETSC_NULL);CHKERRQ(ierr);
|
ierr = PetscOptionsInt("-svd_nsv","Number of singular values to compute","SVDSetDimensions",svd->nsv,&i,PETSC_NULL);CHKERRQ(ierr);
|
ierr = PetscOptionsInt("-svd_ncv","Number of basis vectors","SVDSetDimensions",svd->ncv,&j,PETSC_NULL);CHKERRQ(ierr);
|
ierr = PetscOptionsInt("-svd_ncv","Number of basis vectors","SVDSetDimensions",svd->ncv,&j,PETSC_NULL);CHKERRQ(ierr);
|
ierr = PetscOptionsInt("-svd_mpd","Maximum dimension of projected problem","SVDSetDimensions",svd->mpd,&k,PETSC_NULL);CHKERRQ(ierr);
|
ierr = PetscOptionsInt("-svd_mpd","Maximum dimension of projected problem","SVDSetDimensions",svd->mpd,&k,PETSC_NULL);CHKERRQ(ierr);
|
ierr = SVDSetDimensions(svd,i,j,k);CHKERRQ(ierr);
|
ierr = SVDSetDimensions(svd,i,j,k);CHKERRQ(ierr);
|
|
|
ierr = PetscOptionsTruthGroupBegin("-svd_largest","compute largest singular values","SVDSetWhichSingularTriplets",&flg);CHKERRQ(ierr);
|
ierr = PetscOptionsTruthGroupBegin("-svd_largest","compute largest singular values","SVDSetWhichSingularTriplets",&flg);CHKERRQ(ierr);
|
if (flg) { ierr = SVDSetWhichSingularTriplets(svd,SVD_LARGEST);CHKERRQ(ierr); }
|
if (flg) { ierr = SVDSetWhichSingularTriplets(svd,SVD_LARGEST);CHKERRQ(ierr); }
|
ierr = PetscOptionsTruthGroupEnd("-svd_smallest","compute smallest singular values","SVDSetWhichSingularTriplets",&flg);CHKERRQ(ierr);
|
ierr = PetscOptionsTruthGroupEnd("-svd_smallest","compute smallest singular values","SVDSetWhichSingularTriplets",&flg);CHKERRQ(ierr);
|
if (flg) { ierr = SVDSetWhichSingularTriplets(svd,SVD_SMALLEST);CHKERRQ(ierr); }
|
if (flg) { ierr = SVDSetWhichSingularTriplets(svd,SVD_SMALLEST);CHKERRQ(ierr); }
|
|
|
flg = PETSC_FALSE;
|
flg = PETSC_FALSE;
|
ierr = PetscOptionsTruth("-svd_monitor_cancel","Remove any hardwired monitor routines","SVDMonitorCancel",flg,&flg,PETSC_NULL);CHKERRQ(ierr);
|
ierr = PetscOptionsTruth("-svd_monitor_cancel","Remove any hardwired monitor routines","SVDMonitorCancel",flg,&flg,PETSC_NULL);CHKERRQ(ierr);
|
if (flg) {
|
if (flg) {
|
ierr = SVDMonitorCancel(svd);CHKERRQ(ierr);
|
ierr = SVDMonitorCancel(svd);CHKERRQ(ierr);
|
}
|
}
|
|
|
ierr = PetscOptionsString("-svd_monitor_all","Monitor approximate singular values and error estimates","SVDMonitorSet","stdout",monfilename,PETSC_MAX_PATH_LEN,&flg);CHKERRQ(ierr);
|
ierr = PetscOptionsString("-svd_monitor_all","Monitor approximate singular values and error estimates","SVDMonitorSet","stdout",monfilename,PETSC_MAX_PATH_LEN,&flg);CHKERRQ(ierr);
|
if (flg) {
|
if (flg) {
|
ierr = PetscViewerASCIIMonitorCreate(((PetscObject)svd)->comm,monfilename,((PetscObject)svd)->tablevel,&monviewer);CHKERRQ(ierr);
|
ierr = PetscViewerASCIIMonitorCreate(((PetscObject)svd)->comm,monfilename,((PetscObject)svd)->tablevel,&monviewer);CHKERRQ(ierr);
|
ierr = SVDMonitorSet(svd,SVDMonitorAll,monviewer,(PetscErrorCode (*)(void*))PetscViewerASCIIMonitorDestroy);CHKERRQ(ierr);
|
ierr = SVDMonitorSet(svd,SVDMonitorAll,monviewer,(PetscErrorCode (*)(void*))PetscViewerASCIIMonitorDestroy);CHKERRQ(ierr);
|
ierr = SVDSetTrackAll(svd,PETSC_TRUE);CHKERRQ(ierr);
|
ierr = SVDSetTrackAll(svd,PETSC_TRUE);CHKERRQ(ierr);
|
}
|
}
|
ierr = PetscOptionsString("-svd_monitor_conv","Monitor approximate singular values and error estimates as they converge","SVDMonitorSet","stdout",monfilename,PETSC_MAX_PATH_LEN,&flg);CHKERRQ(ierr);
|
ierr = PetscOptionsString("-svd_monitor_conv","Monitor approximate singular values and error estimates as they converge","SVDMonitorSet","stdout",monfilename,PETSC_MAX_PATH_LEN,&flg);CHKERRQ(ierr);
|
if (flg) {
|
if (flg) {
|
ierr = PetscNew(SVDMONITOR_CONV,&ctx);CHKERRQ(ierr);
|
ierr = PetscNew(SVDMONITOR_CONV,&ctx);CHKERRQ(ierr);
|
ierr = PetscViewerASCIIMonitorCreate(((PetscObject)svd)->comm,monfilename,((PetscObject)svd)->tablevel,&ctx->viewer);CHKERRQ(ierr);
|
ierr = PetscViewerASCIIMonitorCreate(((PetscObject)svd)->comm,monfilename,((PetscObject)svd)->tablevel,&ctx->viewer);CHKERRQ(ierr);
|
ierr = SVDMonitorSet(svd,SVDMonitorConverged,ctx,(PetscErrorCode (*)(void*))SVDMonitorDestroy_Converged);CHKERRQ(ierr);
|
ierr = SVDMonitorSet(svd,SVDMonitorConverged,ctx,(PetscErrorCode (*)(void*))SVDMonitorDestroy_Converged);CHKERRQ(ierr);
|
}
|
}
|
ierr = PetscOptionsString("-svd_monitor","Monitor first unconverged approximate singular value and error estimate","SVDMonitorSet","stdout",monfilename,PETSC_MAX_PATH_LEN,&flg);CHKERRQ(ierr);
|
ierr = PetscOptionsString("-svd_monitor","Monitor first unconverged approximate singular value and error estimate","SVDMonitorSet","stdout",monfilename,PETSC_MAX_PATH_LEN,&flg);CHKERRQ(ierr);
|
if (flg) {
|
if (flg) {
|
ierr = PetscViewerASCIIMonitorCreate(((PetscObject)svd)->comm,monfilename,((PetscObject)svd)->tablevel,&monviewer);CHKERRQ(ierr);
|
ierr = PetscViewerASCIIMonitorCreate(((PetscObject)svd)->comm,monfilename,((PetscObject)svd)->tablevel,&monviewer);CHKERRQ(ierr);
|
ierr = SVDMonitorSet(svd,SVDMonitorFirst,monviewer,(PetscErrorCode (*)(void*))PetscViewerASCIIMonitorDestroy);CHKERRQ(ierr);
|
ierr = SVDMonitorSet(svd,SVDMonitorFirst,monviewer,(PetscErrorCode (*)(void*))PetscViewerASCIIMonitorDestroy);CHKERRQ(ierr);
|
}
|
}
|
flg = PETSC_FALSE;
|
flg = PETSC_FALSE;
|
ierr = PetscOptionsTruth("-svd_monitor_draw","Monitor first unconverged approximate singular value and error estimate graphically","SVDMonitorSet",flg,&flg,PETSC_NULL);CHKERRQ(ierr);
|
ierr = PetscOptionsTruth("-svd_monitor_draw","Monitor first unconverged approximate singular value and error estimate graphically","SVDMonitorSet",flg,&flg,PETSC_NULL);CHKERRQ(ierr);
|
if (flg) {
|
if (flg) {
|
ierr = SVDMonitorSet(svd,SVDMonitorLG,PETSC_NULL,PETSC_NULL);CHKERRQ(ierr);
|
ierr = SVDMonitorSet(svd,SVDMonitorLG,PETSC_NULL,PETSC_NULL);CHKERRQ(ierr);
|
}
|
}
|
flg = PETSC_FALSE;
|
flg = PETSC_FALSE;
|
ierr = PetscOptionsTruth("-svd_monitor_draw_all","Monitor error estimates graphically","SVDMonitorSet",flg,&flg,PETSC_NULL);CHKERRQ(ierr);
|
ierr = PetscOptionsTruth("-svd_monitor_draw_all","Monitor error estimates graphically","SVDMonitorSet",flg,&flg,PETSC_NULL);CHKERRQ(ierr);
|
if (flg) {
|
if (flg) {
|
ierr = SVDMonitorSet(svd,SVDMonitorLGAll,PETSC_NULL,PETSC_NULL);CHKERRQ(ierr);
|
ierr = SVDMonitorSet(svd,SVDMonitorLGAll,PETSC_NULL,PETSC_NULL);CHKERRQ(ierr);
|
ierr = SVDSetTrackAll(svd,PETSC_TRUE);CHKERRQ(ierr);
|
ierr = SVDSetTrackAll(svd,PETSC_TRUE);CHKERRQ(ierr);
|
}
|
}
|
|
|
ierr = PetscOptionsEnd();CHKERRQ(ierr);
|
ierr = PetscOptionsEnd();CHKERRQ(ierr);
|
ierr = IPSetFromOptions(svd->ip);CHKERRQ(ierr);
|
ierr = IPSetFromOptions(svd->ip);CHKERRQ(ierr);
|
if (svd->ops->setfromoptions) {
|
if (svd->ops->setfromoptions) {
|
ierr = (*svd->ops->setfromoptions)(svd);CHKERRQ(ierr);
|
ierr = (*svd->ops->setfromoptions)(svd);CHKERRQ(ierr);
|
}
|
}
|
PetscFunctionReturn(0);
|
PetscFunctionReturn(0);
|
}
|
}
|
|
|
#undef __FUNCT__
|
#undef __FUNCT__
|
#define __FUNCT__ "SVDSetTrackAll"
|
#define __FUNCT__ "SVDSetTrackAll"
|
/*@
|
/*@
|
SVDSetTrackAll - Specifies if the solver must compute the residual norm of all
|
SVDSetTrackAll - Specifies if the solver must compute the residual norm of all
|
approximate singular value or not.
|
approximate singular value or not.
|
|
|
Collective on SVD
|
Collective on SVD
|
|
|
Input Parameters:
|
Input Parameters:
|
+ svd - the singular value solver context
|
+ svd - the singular value solver context
|
- trackall - whether to compute all residuals or not
|
- trackall - whether to compute all residuals or not
|
|
|
Notes:
|
Notes:
|
If the user sets trackall=PETSC_TRUE then the solver computes (or estimates)
|
If the user sets trackall=PETSC_TRUE then the solver computes (or estimates)
|
the residual norm for each singular value approximation. Computing the residual is
|
the residual norm for each singular value approximation. Computing the residual is
|
usually an expensive operation and solvers commonly compute only the residual
|
usually an expensive operation and solvers commonly compute only the residual
|
associated to the first unconverged singular value.
|
associated to the first unconverged singular value.
|
|
|
The options '-svd_monitor_all' and '-svd_monitor_draw_all' automatically
|
The options '-svd_monitor_all' and '-svd_monitor_draw_all' automatically
|
activate this option.
|
activate this option.
|
|
|
Level: intermediate
|
Level: intermediate
|
|
|
.seealso: SVDGetTrackAll()
|
.seealso: SVDGetTrackAll()
|
@*/
|
@*/
|
PetscErrorCode SVDSetTrackAll(SVD svd,PetscTruth trackall)
|
PetscErrorCode SVDSetTrackAll(SVD svd,PetscTruth trackall)
|
{
|
{
|
PetscFunctionBegin;
|
PetscFunctionBegin;
|
PetscValidHeaderSpecific(svd,SVD_COOKIE,1);
|
PetscValidHeaderSpecific(svd,SVD_COOKIE,1);
|
svd->trackall = trackall;
|
svd->trackall = trackall;
|
PetscFunctionReturn(0);
|
PetscFunctionReturn(0);
|
}
|
}
|
|
|
#undef __FUNCT__
|
#undef __FUNCT__
|
#define __FUNCT__ "SVDGetTrackAll"
|
#define __FUNCT__ "SVDGetTrackAll"
|
/*@
|
/*@
|
SVDGetTrackAll - Returns the flag indicating whether all residual norms must
|
SVDGetTrackAll - Returns the flag indicating whether all residual norms must
|
be computed or not.
|
be computed or not.
|
|
|
Not Collective
|
Not Collective
|
|
|
Input Parameter:
|
Input Parameter:
|
. svd - the singular value solver context
|
. svd - the singular value solver context
|
|
|
Output Parameter:
|
Output Parameter:
|
. trackall - the returned flag
|
. trackall - the returned flag
|
|
|
Level: intermediate
|
Level: intermediate
|
|
|
.seealso: SVDSetTrackAll()
|
.seealso: SVDSetTrackAll()
|
@*/
|
@*/
|
PetscErrorCode SVDGetTrackAll(SVD svd,PetscTruth *trackall)
|
PetscErrorCode SVDGetTrackAll(SVD svd,PetscTruth *trackall)
|
{
|
{
|
PetscFunctionBegin;
|
PetscFunctionBegin;
|
PetscValidHeaderSpecific(svd,SVD_COOKIE,1);
|
PetscValidHeaderSpecific(svd,SVD_COOKIE,1);
|
PetscValidPointer(trackall,2);
|
PetscValidPointer(trackall,2);
|
*trackall = svd->trackall;
|
*trackall = svd->trackall;
|
PetscFunctionReturn(0);
|
PetscFunctionReturn(0);
|
}
|
}
|
|
|
|
|
#undef __FUNCT__
|
#undef __FUNCT__
|
#define __FUNCT__ "SVDSetOptionsPrefix"
|
#define __FUNCT__ "SVDSetOptionsPrefix"
|
/*@C
|
/*@C
|
SVDSetOptionsPrefix - Sets the prefix used for searching for all
|
SVDSetOptionsPrefix - Sets the prefix used for searching for all
|
SVD options in the database.
|
SVD options in the database.
|
|
|
Collective on SVD
|
Collective on SVD
|
|
|
Input Parameters:
|
Input Parameters:
|
+ svd - the singular value solver context
|
+ svd - the singular value solver context
|
- prefix - the prefix string to prepend to all SVD option requests
|
- prefix - the prefix string to prepend to all SVD option requests
|
|
|
Notes:
|
Notes:
|
A hyphen (-) must NOT be given at the beginning of the prefix name.
|
A hyphen (-) must NOT be given at the beginning of the prefix name.
|
The first character of all runtime options is AUTOMATICALLY the
|
The first character of all runtime options is AUTOMATICALLY the
|
hyphen.
|
hyphen.
|
|
|
For example, to distinguish between the runtime options for two
|
For example, to distinguish between the runtime options for two
|
different SVD contexts, one could call
|
different SVD contexts, one could call
|
.vb
|
.vb
|
SVDSetOptionsPrefix(svd1,"svd1_")
|
SVDSetOptionsPrefix(svd1,"svd1_")
|
SVDSetOptionsPrefix(svd2,"svd2_")
|
SVDSetOptionsPrefix(svd2,"svd2_")
|
.ve
|
.ve
|
|
|
Level: advanced
|
Level: advanced
|
|
|
.seealso: SVDAppendOptionsPrefix(), SVDGetOptionsPrefix()
|
.seealso: SVDAppendOptionsPrefix(), SVDGetOptionsPrefix()
|
@*/
|
@*/
|
PetscErrorCode SVDSetOptionsPrefix(SVD svd,const char *prefix)
|
PetscErrorCode SVDSetOptionsPrefix(SVD svd,const char *prefix)
|
{
|
{
|
PetscErrorCode ierr;
|
PetscErrorCode ierr;
|
PetscTruth flg1,flg2;
|
PetscTruth flg1,flg2;
|
EPS eps;
|
EPS eps;
|
|
|
PetscFunctionBegin;
|
PetscFunctionBegin;
|
PetscValidHeaderSpecific(svd,SVD_COOKIE,1);
|
PetscValidHeaderSpecific(svd,SVD_COOKIE,1);
|
ierr = PetscObjectSetOptionsPrefix((PetscObject)svd,prefix);CHKERRQ(ierr);
|
ierr = PetscObjectSetOptionsPrefix((PetscObject)svd,prefix);CHKERRQ(ierr);
|
ierr = PetscTypeCompare((PetscObject)svd,SVDCROSS,&flg1);CHKERRQ(ierr);
|
ierr = PetscTypeCompare((PetscObject)svd,SVDCROSS,&flg1);CHKERRQ(ierr);
|
ierr = PetscTypeCompare((PetscObject)svd,SVDCYCLIC,&flg2);CHKERRQ(ierr);
|
ierr = PetscTypeCompare((PetscObject)svd,SVDCYCLIC,&flg2);CHKERRQ(ierr);
|
if (flg1) {
|
if (flg1) {
|
ierr = SVDCrossGetEPS(svd,&eps);CHKERRQ(ierr);
|
ierr = SVDCrossGetEPS(svd,&eps);CHKERRQ(ierr);
|
} else if (flg2) {
|
} else if (flg2) {
|
ierr = SVDCyclicGetEPS(svd,&eps);CHKERRQ(ierr);
|
ierr = SVDCyclicGetEPS(svd,&eps);CHKERRQ(ierr);
|
}
|
}
|
if (flg1 || flg2) {
|
if (flg1 || flg2) {
|
ierr = EPSSetOptionsPrefix(eps,prefix);CHKERRQ(ierr);
|
ierr = EPSSetOptionsPrefix(eps,prefix);CHKERRQ(ierr);
|
ierr = EPSAppendOptionsPrefix(eps,"svd_");CHKERRQ(ierr);
|
ierr = EPSAppendOptionsPrefix(eps,"svd_");CHKERRQ(ierr);
|
}
|
}
|
ierr = IPSetOptionsPrefix(svd->ip,prefix);CHKERRQ(ierr);
|
ierr = IPSetOptionsPrefix(svd->ip,prefix);CHKERRQ(ierr);
|
ierr = IPAppendOptionsPrefix(svd->ip,"svd_");CHKERRQ(ierr);
|
ierr = IPAppendOptionsPrefix(svd->ip,"svd_");CHKERRQ(ierr);
|
PetscFunctionReturn(0);
|
PetscFunctionReturn(0);
|
}
|
}
|
|
|
#undef __FUNCT__
|
#undef __FUNCT__
|
#define __FUNCT__ "SVDAppendOptionsPrefix"
|
#define __FUNCT__ "SVDAppendOptionsPrefix"
|
/*@C
|
/*@C
|
SVDAppendOptionsPrefix - Appends to the prefix used for searching for all
|
SVDAppendOptionsPrefix - Appends to the prefix used for searching for all
|
SVD options in the database.
|
SVD options in the database.
|
|
|
Collective on SVD
|
Collective on SVD
|
|
|
Input Parameters:
|
Input Parameters:
|
+ svd - the singular value solver context
|
+ svd - the singular value solver context
|
- prefix - the prefix string to prepend to all SVD option requests
|
- prefix - the prefix string to prepend to all SVD option requests
|
|
|
Notes:
|
Notes:
|
A hyphen (-) must NOT be given at the beginning of the prefix name.
|
A hyphen (-) must NOT be given at the beginning of the prefix name.
|
The first character of all runtime options is AUTOMATICALLY the hyphen.
|
The first character of all runtime options is AUTOMATICALLY the hyphen.
|
|
|
Level: advanced
|
Level: advanced
|
|
|
.seealso: SVDSetOptionsPrefix(), SVDGetOptionsPrefix()
|
.seealso: SVDSetOptionsPrefix(), SVDGetOptionsPrefix()
|
@*/
|
@*/
|
PetscErrorCode SVDAppendOptionsPrefix(SVD svd,const char *prefix)
|
PetscErrorCode SVDAppendOptionsPrefix(SVD svd,const char *prefix)
|
{
|
{
|
PetscErrorCode ierr;
|
PetscErrorCode ierr;
|
PetscTruth flg1,flg2;
|
PetscTruth flg1,flg2;
|
EPS eps;
|
EPS eps;
|
|
|
PetscFunctionBegin;
|
PetscFunctionBegin;
|
PetscValidHeaderSpecific(svd,SVD_COOKIE,1);
|
PetscValidHeaderSpecific(svd,SVD_COOKIE,1);
|
ierr = PetscObjectAppendOptionsPrefix((PetscObject)svd,prefix);CHKERRQ(ierr);
|
ierr = PetscObjectAppendOptionsPrefix((PetscObject)svd,prefix);CHKERRQ(ierr);
|
ierr = PetscTypeCompare((PetscObject)svd,SVDCROSS,&flg1);CHKERRQ(ierr);
|
ierr = PetscTypeCompare((PetscObject)svd,SVDCROSS,&flg1);CHKERRQ(ierr);
|
ierr = PetscTypeCompare((PetscObject)svd,SVDCYCLIC,&flg2);CHKERRQ(ierr);
|
ierr = PetscTypeCompare((PetscObject)svd,SVDCYCLIC,&flg2);CHKERRQ(ierr);
|
if (flg1) {
|
if (flg1) {
|
ierr = SVDCrossGetEPS(svd,&eps);CHKERRQ(ierr);
|
ierr = SVDCrossGetEPS(svd,&eps);CHKERRQ(ierr);
|
} else if (flg2) {
|
} else if (flg2) {
|
ierr = SVDCyclicGetEPS(svd,&eps);CHKERRQ(ierr);
|
ierr = SVDCyclicGetEPS(svd,&eps);CHKERRQ(ierr);
|
}
|
}
|
if (flg1 || flg2) {
|
if (flg1 || flg2) {
|
ierr = EPSSetOptionsPrefix(eps,((PetscObject)svd)->prefix);CHKERRQ(ierr);
|
ierr = EPSSetOptionsPrefix(eps,((PetscObject)svd)->prefix);CHKERRQ(ierr);
|
ierr = EPSAppendOptionsPrefix(eps,"svd_");CHKERRQ(ierr);
|
ierr = EPSAppendOptionsPrefix(eps,"svd_");CHKERRQ(ierr);
|
}
|
}
|
ierr = IPSetOptionsPrefix(svd->ip,((PetscObject)svd)->prefix);CHKERRQ(ierr);
|
ierr = IPSetOptionsPrefix(svd->ip,((PetscObject)svd)->prefix);CHKERRQ(ierr);
|
ierr = IPAppendOptionsPrefix(svd->ip,"svd_");CHKERRQ(ierr);
|
ierr = IPAppendOptionsPrefix(svd->ip,"svd_");CHKERRQ(ierr);
|
PetscFunctionReturn(0);
|
PetscFunctionReturn(0);
|
}
|
}
|
|
|
#undef __FUNCT__
|
#undef __FUNCT__
|
#define __FUNCT__ "SVDGetOptionsPrefix"
|
#define __FUNCT__ "SVDGetOptionsPrefix"
|
/*@C
|
/*@C
|
SVDGetOptionsPrefix - Gets the prefix used for searching for all
|
SVDGetOptionsPrefix - Gets the prefix used for searching for all
|
SVD options in the database.
|
SVD options in the database.
|
|
|
Not Collective
|
Not Collective
|
|
|
Input Parameters:
|
Input Parameters:
|
. svd - the singular value solver context
|
. svd - the singular value solver context
|
|
|
Output Parameters:
|
Output Parameters:
|
. prefix - pointer to the prefix string used is returned
|
. prefix - pointer to the prefix string used is returned
|
|
|
Notes: On the fortran side, the user should pass in a string 'prefix' of
|
Notes: On the fortran side, the user should pass in a string 'prefix' of
|
sufficient length to hold the prefix.
|
sufficient length to hold the prefix.
|
|
|
Level: advanced
|
Level: advanced
|
|
|
.seealso: SVDSetOptionsPrefix(), SVDAppendOptionsPrefix()
|
.seealso: SVDSetOptionsPrefix(), SVDAppendOptionsPrefix()
|
@*/
|
@*/
|
PetscErrorCode SVDGetOptionsPrefix(SVD svd,const char *prefix[])
|
PetscErrorCode SVDGetOptionsPrefix(SVD svd,const char *prefix[])
|
{
|
{
|
PetscErrorCode ierr;
|
PetscErrorCode ierr;
|
PetscFunctionBegin;
|
PetscFunctionBegin;
|
PetscValidHeaderSpecific(svd,SVD_COOKIE,1);
|
PetscValidHeaderSpecific(svd,SVD_COOKIE,1);
|
PetscValidPointer(prefix,2);
|
PetscValidPointer(prefix,2);
|
ierr = PetscObjectGetOptionsPrefix((PetscObject)svd,prefix);CHKERRQ(ierr);
|
ierr = PetscObjectGetOptionsPrefix((PetscObject)svd,prefix);CHKERRQ(ierr);
|
PetscFunctionReturn(0);
|
PetscFunctionReturn(0);
|
}
|
}
|
|
|
|
|
|
|