1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
|
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML3.2 EN">
<HTML>
<HEAD> <link rel="canonical" href="http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/KSP/KSPFischerGuessCreate.html" />
<META NAME="GENERATOR" CONTENT="DOCTEXT">
<TITLE>KSPFischerGuessCreate</TITLE>
</HEAD>
<BODY BGCOLOR="FFFFFF">
<div id="version" align=right><b>petsc-3.4.2 2013-07-02</b></div>
<A NAME="KSPFischerGuessCreate"><H1>KSPFischerGuessCreate</H1></A>
Implements Paul Fischer's initial guess algorithm Method 1 and 2 for situations where a linear system is solved repeatedly
<H3><FONT COLOR="#CC3333">Synopsis</FONT></H3>
<PRE>
PetscErrorCode KSPFischerGuessCreate(KSP ksp,PetscInt method,PetscInt maxl,KSPFischerGuess *itg)
</PRE>
<H3><FONT COLOR="#CC3333">References</FONT></H3>
http://ntrs.nasa.gov/archive/nasa/casi.ntrs.nasa.gov/19940020363_1994020363.pdf
<P>
Notes: the algorithm is different from the paper because we do not CHANGE the right hand side of the new
problem and solve the problem with an initial guess of zero, rather we solve the original new problem
with a nonzero initial guess (this is done so that the linear solver convergence tests are based on
the original RHS.) But we use the xtilde = x - xguess as the new direction so that it is not
mostly orthogonal to the previous solutions.
<P>
These are not intended to be used directly, they are called by <A HREF="../KSP/KSP.html#KSP">KSP</A> automatically when the
<A HREF="../KSP/KSP.html#KSP">KSP</A> option <A HREF="../KSP/KSPSetFischerGuess.html#KSPSetFischerGuess">KSPSetFischerGuess</A>(<A HREF="../KSP/KSP.html#KSP">KSP</A>,<A HREF="../Sys/PetscInt.html#PetscInt">PetscInt</A>,<A HREF="../Sys/PetscInt.html#PetscInt">PetscInt</A>) or -ksp_guess_fischer <int,int>
<P>
Method 2 is only for positive definite matrices, since it uses the A norm.
<P>
This is not currently programmed as a PETSc class because there are only two methods; if more methods
are introduced it should be changed. For example the Knoll guess should be included
<P>
<P>
<P><B><P><B><FONT COLOR="#CC3333">Level:</FONT></B>advanced
<BR><FONT COLOR="#CC3333">Location:</FONT></B><A HREF="../../../src/ksp/ksp/interface/iguess.c.html#KSPFischerGuessCreate">src/ksp/ksp/interface/iguess.c</A>
<BR><A HREF="./index.html">Index of all KSP routines</A>
<BR><A HREF="../../index.html">Table of Contents for all manual pages</A>
<BR><A HREF="../singleindex.html">Index of all manual pages</A>
</BODY></HTML>
|
