EXTERN_C_END
#undef __FUNCT__
#define __FUNCT__ "TaoAppSetFiniteDifferencesOptions"
/* @
TaoAppSetFiniteDifferencesOptions - Sets various TAO parameters from user options
Collective on TAO_APPLICATION
Input Parameters:
+ taoapp - the TAO Application (optional)
Level: beginner
.keywords: options, finite differences
.seealso: TaoSolveApplication();
@ */
int TaoAppSetFiniteDifferencesOptions(TAO_APPLICATION taoapp){
int info;
PetscTruth flg;
PetscFunctionBegin;
PetscValidHeaderSpecific(taoapp,TAO_APP_COOKIE,1);
flg=PETSC_FALSE;
info = PetscOptionsName("-tao_fd","use finite differences for Hessian","TaoAppDefaultComputeHessian",&flg);CHKERRQ(info);
if (flg) {
info = TaoAppSetHessianRoutine(taoapp,TaoAppDefaultComputeHessian,(void*)taoapp);CHKERRQ(info);
info = PetscInfo(taoapp,"Setting default finite difference Hessian matrix\n"); CHKERRQ(info);
}
flg=PETSC_FALSE;
info = PetscOptionsName("-tao_fdgrad","use finite differences for gradient","TaoAppDefaultComputeGradient",&flg);CHKERRQ(info);
if (flg) {
info = TaoAppSetGradientRoutine(taoapp,TaoAppDefaultComputeGradient,(void*)taoapp);CHKERRQ(info);
info = PetscInfo(taoapp,"Setting default finite difference gradient routine\n"); CHKERRQ(info);
}
flg=PETSC_FALSE;
info = PetscOptionsName("-tao_fd_coloring","use finite differences with coloring to compute Hessian","TaoAppDefaultComputeHessianColor",&flg);CHKERRQ(info);
if (flg) {
info = TaoAppSetHessianRoutine(taoapp,TaoAppDefaultComputeHessianColor,(void*)taoapp);CHKERRQ(info);
info = PetscInfo(taoapp,"Use finite differencing with coloring to compute Hessian \n"); CHKERRQ(info);
}
PetscFunctionReturn(0);
}
static PetscErrorCode PCSetFromOptions_ILU(PetscOptionItems *PetscOptionsObject,PC pc)
{
PetscErrorCode ierr;
PetscInt itmp;
PetscBool flg,set;
PC_ILU *ilu = (PC_ILU*)pc->data;
PetscReal tol;
PetscFunctionBegin;
ierr = PetscOptionsHead(PetscOptionsObject,"ILU Options");CHKERRQ(ierr);
ierr = PCSetFromOptions_Factor(PetscOptionsObject,pc);CHKERRQ(ierr);
ierr = PetscOptionsInt("-pc_factor_levels","levels of fill","PCFactorSetLevels",(PetscInt)((PC_Factor*)ilu)->info.levels,&itmp,&flg);CHKERRQ(ierr);
if (flg) ((PC_Factor*)ilu)->info.levels = itmp;
ierr = PetscOptionsBool("-pc_factor_diagonal_fill","Allow fill into empty diagonal entry","PCFactorSetAllowDiagonalFill",((PC_Factor*)ilu)->info.diagonal_fill ? PETSC_TRUE : PETSC_FALSE,&flg,&set);CHKERRQ(ierr);
if (set) ((PC_Factor*)ilu)->info.diagonal_fill = (PetscReal) flg;
ierr = PetscOptionsName("-pc_factor_nonzeros_along_diagonal","Reorder to remove zeros from diagonal","PCFactorReorderForNonzeroDiagonal",&flg);CHKERRQ(ierr);
if (flg) {
tol = PETSC_DECIDE;
ierr = PetscOptionsReal("-pc_factor_nonzeros_along_diagonal","Reorder to remove zeros from diagonal","PCFactorReorderForNonzeroDiagonal",ilu->nonzerosalongdiagonaltol,&tol,0);CHKERRQ(ierr);
ierr = PCFactorReorderForNonzeroDiagonal(pc,tol);CHKERRQ(ierr);
}
ierr = PetscOptionsTail();CHKERRQ(ierr);
PetscFunctionReturn(0);
}
/*@C
PetscDrawSetFromOptions - Sets the graphics type from the options database.
Defaults to a PETSc X windows graphics.
Collective on PetscDraw
Input Parameter:
. draw - the graphics context
Options Database Keys:
+ -nox - do not use X graphics (ignore graphics calls, but run program correctly)
- -nox_warning - when X windows support is not installed this prevents the warning message
from being printed
Level: intermediate
Notes:
Must be called after PetscDrawCreate() before the PetscDrawtor is used.
Concepts: drawing^setting options
Concepts: graphics^setting options
.seealso: PetscDrawCreate(), PetscDrawSetType()
@*/
PetscErrorCode PetscDrawSetFromOptions(PetscDraw draw)
{
PetscErrorCode ierr;
PetscBool flg,nox;
char vtype[256];
const char *def;
PetscBool save;
#if !defined(PETSC_USE_WINDOWS_GRAPHICS) && !defined(PETSC_HAVE_X)
PetscBool warn;
#endif
char filename[PETSC_MAX_PATH_LEN];
PetscBool movie = PETSC_FALSE;
PetscFunctionBegin;
PetscValidHeaderSpecific(draw,PETSC_DRAW_CLASSID,1);
if (!PetscDrawList) {
ierr = PetscDrawRegisterAll(PETSC_NULL);CHKERRQ(ierr);
}
if (((PetscObject)draw)->type_name) {
def = ((PetscObject)draw)->type_name;
} else {
ierr = PetscOptionsHasName(PETSC_NULL,"-nox",&nox);CHKERRQ(ierr);
def = PETSC_DRAW_NULL;
#if defined(PETSC_USE_WINDOWS_GRAPHICS) && !defined(PETSC_HAVE_X)
if (!nox) def = PETSC_DRAW_WIN32;
#elif defined(PETSC_HAVE_X)
if (!nox) def = PETSC_DRAW_X;
#else
ierr = PetscOptionsHasName(PETSC_NULL,"-nox_warning",&warn);CHKERRQ(ierr);
if (!nox && !warn) {
(*PetscErrorPrintf)("PETSc installed without X windows on this machine\nproceeding without graphics\n");
}
#endif
}
ierr = PetscObjectOptionsBegin((PetscObject)draw);CHKERRQ(ierr);
ierr = PetscOptionsList("-draw_type","Type of graphical output","PetscDrawSetType",PetscDrawList,def,vtype,256,&flg);CHKERRQ(ierr);
if (flg) {
ierr = PetscDrawSetType(draw,vtype);CHKERRQ(ierr);
} else if (!((PetscObject)draw)->type_name) {
ierr = PetscDrawSetType(draw,def);CHKERRQ(ierr);
}
ierr = PetscOptionsName("-nox","Run without graphics","None",&nox);CHKERRQ(ierr);
ierr = PetscOptionsBool("-draw_save_movie","Make a movie from the images saved","PetscDrawSetSave",movie,&movie,PETSC_NULL);CHKERRQ(ierr);
ierr = PetscOptionsString("-draw_save","Save graphics to file","PetscDrawSetSave",filename,filename,PETSC_MAX_PATH_LEN,&save);CHKERRQ(ierr);
if (save) {
ierr = PetscDrawSetSave(draw,filename,movie);CHKERRQ(ierr);
}
/* process any options handlers added with PetscObjectAddOptionsHandler() */
ierr = PetscObjectProcessOptionsHandlers((PetscObject)draw);CHKERRQ(ierr);
ierr = PetscOptionsEnd();CHKERRQ(ierr);
PetscFunctionReturn(0);
}
PetscErrorCode KSPSetFromOptions_LSQR(KSP ksp)
{
PetscErrorCode ierr;
KSP_LSQR *lsqr = (KSP_LSQR*)ksp->data;
PetscFunctionBegin;
ierr = PetscOptionsHead("KSP LSQR Options");CHKERRQ(ierr);
ierr = PetscOptionsName("-ksp_LSQR_set_standard_error","Set Standard Error Estimates of Solution","KSPLSQRSetStandardErrorVec",&lsqr->se_flg);CHKERRQ(ierr);
ierr = PetscOptionsTail();CHKERRQ(ierr);
PetscFunctionReturn(0);
}
/*@
TaoSetOptions - Sets various TAO parameters from user options
Collective on TAO_APPLICATION
Input Parameters:
+ taoapp - the TAO Application (optional)
- tao - the TAO optimization solver (optional)
Level: beginner
Note:
This routine should be called after TaoSetupApplicationSolver()
Note:
This routine must be called if there are multiple processors involved and
the MPI Communicator is different than MPI_COMM_WORLD.
.keywords: options
.seealso: TaoSolveApplication()
@*/
int TaoSetOptions(TAO_APPLICATION taoapp, TAO_SOLVER tao){
int info;
const char *prefix=0;
PetscTruth flg;
MPI_Comm comm=MPI_COMM_WORLD;
PetscFunctionBegin;
if (tao){
PetscValidHeaderSpecific(tao,TAO_COOKIE,2);
info = PetscObjectGetOptionsPrefix((PetscObject)tao,&prefix); CHKERRQ(info);
info = PetscObjectGetComm((PetscObject)tao,&comm);CHKERRQ(info);
info = PetscOptionsBegin(comm,prefix,"TAO PETSC APPLICATIONS ","solver");CHKERRQ(info);
info = TaoSetFromOptions(tao); CHKERRQ(info);
flg=PETSC_FALSE;
info = PetscOptionsName("-tao_xmonitor","Use graphics convergence","TaoPetscXMonitor",&flg);CHKERRQ(info);
if (flg){
info = TaoSetPetscXMonitor(tao); CHKERRQ(info);
}
info = PetscOptionsEnd();CHKERRQ(info);
}
if (taoapp){
info = TaoAppSetFromOptions(taoapp); CHKERRQ(info);
}
if (tao && taoapp){
info = TaoSetupApplicationSolver(taoapp,tao);CHKERRQ(info);
info = PetscOptionsName("-tao_lmvmh","User supplies approximate hessian for LMVM solvers","TaoLMVMSetH0",&flg);
if (flg){
info=TaoBLMVMSetH0(tao,TAO_TRUE);CHKERRQ(info);
info=TaoLMVMSetH0(tao,TAO_TRUE);CHKERRQ(info);
}
}
PetscFunctionReturn(0);
}
/*@
MatSetFromOptions - Creates a matrix where the type is determined
from the options database. Generates a parallel MPI matrix if the
communicator has more than one processor. The default matrix type is
AIJ, using the routines MatCreateSeqAIJ() and MatCreateAIJ() if
you do not select a type in the options database.
Collective on Mat
Input Parameter:
. A - the matrix
Options Database Keys:
+ -mat_type seqaij - AIJ type, uses MatCreateSeqAIJ()
. -mat_type mpiaij - AIJ type, uses MatCreateAIJ()
. -mat_type seqdense - dense type, uses MatCreateSeqDense()
. -mat_type mpidense - dense type, uses MatCreateDense()
. -mat_type seqbaij - block AIJ type, uses MatCreateSeqBAIJ()
- -mat_type mpibaij - block AIJ type, uses MatCreateBAIJ()
Even More Options Database Keys:
See the manpages for particular formats (e.g., MatCreateSeqAIJ())
for additional format-specific options.
Level: beginner
.keywords: matrix, create
.seealso: MatCreateSeqAIJ((), MatCreateAIJ(),
MatCreateSeqDense(), MatCreateDense(),
MatCreateSeqBAIJ(), MatCreateBAIJ(),
MatCreateSeqSBAIJ(), MatCreateSBAIJ(),
MatConvert()
@*/
PetscErrorCode MatSetFromOptions(Mat B)
{
PetscErrorCode ierr;
const char *deft = MATAIJ;
char type[256];
PetscBool flg,set;
PetscFunctionBegin;
PetscValidHeaderSpecific(B,MAT_CLASSID,1);
ierr = PetscObjectOptionsBegin((PetscObject)B);CHKERRQ(ierr);
if (B->rmap->bs < 0) {
PetscInt newbs = -1;
ierr = PetscOptionsInt("-mat_block_size","Set the blocksize used to store the matrix","MatSetBlockSize",newbs,&newbs,&flg);CHKERRQ(ierr);
if (flg) {
ierr = PetscLayoutSetBlockSize(B->rmap,newbs);CHKERRQ(ierr);
ierr = PetscLayoutSetBlockSize(B->cmap,newbs);CHKERRQ(ierr);
}
}
ierr = PetscOptionsFList("-mat_type","Matrix type","MatSetType",MatList,deft,type,256,&flg);CHKERRQ(ierr);
if (flg) {
ierr = MatSetType(B,type);CHKERRQ(ierr);
} else if (!((PetscObject)B)->type_name) {
ierr = MatSetType(B,deft);CHKERRQ(ierr);
}
ierr = PetscOptionsName("-mat_is_symmetric","Checks if mat is symmetric on MatAssemblyEnd()","MatIsSymmetric",&B->checksymmetryonassembly);CHKERRQ(ierr);
ierr = PetscOptionsReal("-mat_is_symmetric","Checks if mat is symmetric on MatAssemblyEnd()","MatIsSymmetric",0.0,&B->checksymmetrytol,NULL);CHKERRQ(ierr);
ierr = PetscOptionsBool("-mat_null_space_test","Checks if provided null space is correct in MatAssemblyEnd()","MatSetNullSpaceTest",PETSC_FALSE,&B->checknullspaceonassembly,NULL);CHKERRQ(ierr);
if (B->ops->setfromoptions) {
ierr = (*B->ops->setfromoptions)(B);CHKERRQ(ierr);
}
flg = PETSC_FALSE;
ierr = PetscOptionsBool("-mat_new_nonzero_location_err","Generate an error if new nonzeros are created in the matrix structure (useful to test preallocation)","MatSetOption",flg,&flg,&set);CHKERRQ(ierr);
if (set) {ierr = MatSetOption(B,MAT_NEW_NONZERO_LOCATION_ERR,flg);CHKERRQ(ierr);}
flg = PETSC_FALSE;
ierr = PetscOptionsBool("-mat_new_nonzero_allocation_err","Generate an error if new nonzeros are allocated in the matrix structure (useful to test preallocation)","MatSetOption",flg,&flg,&set);CHKERRQ(ierr);
if (set) {ierr = MatSetOption(B,MAT_NEW_NONZERO_ALLOCATION_ERR,flg);CHKERRQ(ierr);}
/* process any options handlers added with PetscObjectAddOptionsHandler() */
ierr = PetscObjectProcessOptionsHandlers((PetscObject)B);CHKERRQ(ierr);
ierr = PetscOptionsEnd();CHKERRQ(ierr);
PetscFunctionReturn(0);
}
PetscErrorCode KSPSetFromOptions_LSQR(KSP ksp)
{
PetscErrorCode ierr;
KSP_LSQR *lsqr = (KSP_LSQR*)ksp->data;
char monfilename[PETSC_MAX_PATH_LEN];
PetscViewer monviewer;
PetscBool flg;
PetscFunctionBegin;
ierr = PetscOptionsHead("KSP LSQR Options");CHKERRQ(ierr);
ierr = PetscOptionsName("-ksp_lsqr_set_standard_error","Set Standard Error Estimates of Solution","KSPLSQRSetStandardErrorVec",&lsqr->se_flg);CHKERRQ(ierr);
ierr = PetscOptionsString("-ksp_lsqr_monitor","Monitor residual norm and norm of residual of normal equations","KSPMonitorSet","stdout",monfilename,PETSC_MAX_PATH_LEN,&flg);CHKERRQ(ierr);
if (flg) {
ierr = PetscViewerASCIIOpen(PetscObjectComm((PetscObject)ksp),monfilename,&monviewer);CHKERRQ(ierr);
ierr = KSPMonitorSet(ksp,KSPLSQRMonitorDefault,monviewer,(PetscErrorCode (*)(void**))PetscViewerDestroy);CHKERRQ(ierr);
}
ierr = PetscOptionsTail();CHKERRQ(ierr);
PetscFunctionReturn(0);
}
EXTERN_C_END
#undef __FUNCT__
#define __FUNCT__ "PCSetFromOptions_ILU"
static PetscErrorCode PCSetFromOptions_ILU(PC pc)
{
PetscErrorCode ierr;
PetscInt itmp;
PetscBool flg;
/* PetscReal dt[3]; */
PC_ILU *ilu = (PC_ILU*)pc->data;
PetscReal tol;
PetscFunctionBegin;
ierr = PetscOptionsHead("ILU Options");CHKERRQ(ierr);
ierr = PCSetFromOptions_Factor(pc);CHKERRQ(ierr);
ierr = PetscOptionsInt("-pc_factor_levels","levels of fill","PCFactorSetLevels",(PetscInt)((PC_Factor*)ilu)->info.levels,&itmp,&flg);CHKERRQ(ierr);
if (flg) ((PC_Factor*)ilu)->info.levels = itmp;
flg = PETSC_FALSE;
ierr = PetscOptionsBool("-pc_factor_diagonal_fill","Allow fill into empty diagonal entry","PCFactorSetAllowDiagonalFill",flg,&flg,PETSC_NULL);CHKERRQ(ierr);
((PC_Factor*)ilu)->info.diagonal_fill = (double) flg;
/*
dt[0] = ((PC_Factor*)ilu)->info.dt;
dt[1] = ((PC_Factor*)ilu)->info.dtcol;
dt[2] = ((PC_Factor*)ilu)->info.dtcount;
PetscInt dtmax = 3;
ierr = PetscOptionsRealArray("-pc_factor_drop_tolerance,","<dt,dtcol,maxrowcount>","PCFactorSetDropTolerance",dt,&dtmax,&flg);CHKERRQ(ierr);
if (flg) {
ierr = PCFactorSetDropTolerance(pc,dt[0],dt[1],(PetscInt)dt[2]);CHKERRQ(ierr);
}
*/
ierr = PetscOptionsName("-pc_factor_nonzeros_along_diagonal","Reorder to remove zeros from diagonal","PCFactorReorderForNonzeroDiagonal",&flg);CHKERRQ(ierr);
if (flg) {
tol = PETSC_DECIDE;
ierr = PetscOptionsReal("-pc_factor_nonzeros_along_diagonal","Reorder to remove zeros from diagonal","PCFactorReorderForNonzeroDiagonal",ilu->nonzerosalongdiagonaltol,&tol,0);CHKERRQ(ierr);
ierr = PCFactorReorderForNonzeroDiagonal(pc,tol);CHKERRQ(ierr);
}
ierr = PetscOptionsTail();CHKERRQ(ierr);
PetscFunctionReturn(0);
}