manual/s_ecco/lsopt_flow_1.tex

{\scriptsize
\begin{verbatim}
  lsopt_top
      |
      |---- check arguments
      |---- CALL INSTORE
      |       |
      |       |---- determine whether OPWARMI available:
      |                * if no:  cold start: create OPWARMI
      |                * if yes: warm start: read from OPWARMI
      |             create or open OPWARMD
      |
      |---- check consistency between OPWARMI and model parameters
      | 
      |---- >>> if COLD start: <<<
      |      |  first simulation with f.g. xx_0; output: first ff_0, gg_0
      |      |  set first preconditioner value xdiff_0 to 1
      |      |  store xx(0), gg(0), xdiff(0) to OPWARMD (first 3 entries)
      |      |
      |     >>> else: WARM start: <<<
      |         read xx(i), gg(i) from OPWARMD (first 2 entries)
      |         for first warm start after cold start, i=0
      |
      |
      |
      |---- /// if ITMAX > 0: perform optimization (increment loop index i)
      |      (
      |      )---- save current values of gg(i-1) -> gold(i-1), ff -> fold(i-1)
      |      (---- CALL LSUPDXX
      |      )       |
      |      (       |---- >>> if jmax=0 <<<
      |      )       |      |  first optimization after cold start:
      |      (       |      |  preconditioner estimated via ff_0 - ff_(first guess)
      |      )       |      |  dd(i-1) = -gg(i-1)*preco
      |      (       |      |  
      |      )       |     >>> if jmax > 0 <<<
      |      (       |         dd(i-1) = -gg(i-1)
      |      )       |         CALL HESSUPD
      |      (       |           |
      |      )       |           |---- dd(i-1) modified via Hessian approx.
      |      (       |
      |      )       |---- >>> if <dd,gg> >= 0 <<<
      |      (       |         ifail = 4
      |      )       |
      |      (       |---- compute step size: tact(i-1)
      |      )       |---- compute update: xdiff(i) = xx(i-1) + tact(i-1)*dd(i-1)
      |      (
      |      )---- >>> if ifail = 4 <<<
      |      (         goto 1000
      |      )
      |      (---- CALL OPTLINE / LSLINE
      |      )       |
     ...    ...     ...
\end{verbatim}
}
1	heimbach	1.1	{\scriptsize
2			\begin{verbatim}
3			lsopt_top
4			\|
5			\|---- check arguments
6			\|---- CALL INSTORE
7			\| \|
8			\| \|---- determine whether OPWARMI available:
9			\| * if no: cold start: create OPWARMI
10			\| * if yes: warm start: read from OPWARMI
11			\| create or open OPWARMD
12			\|
13			\|---- check consistency between OPWARMI and model parameters
14			\|
15			\|---- >>> if COLD start: <<<
16			\| \| first simulation with f.g. xx_0; output: first ff_0, gg_0
17			\| \| set first preconditioner value xdiff_0 to 1
18			\| \| store xx(0), gg(0), xdiff(0) to OPWARMD (first 3 entries)
19			\| \|
20			\| >>> else: WARM start: <<<
21			\| read xx(i), gg(i) from OPWARMD (first 2 entries)
22			\| for first warm start after cold start, i=0
23			\|
24			\|
25			\|
26			\|---- /// if ITMAX > 0: perform optimization (increment loop index i)
27			\| (
28			\| )---- save current values of gg(i-1) -> gold(i-1), ff -> fold(i-1)
29			\| (---- CALL LSUPDXX
30			\| ) \|
31			\| ( \|---- >>> if jmax=0 <<<
32			\| ) \| \| first optimization after cold start:
33			\| ( \| \| preconditioner estimated via ff_0 - ff_(first guess)
34			\| ) \| \| dd(i-1) = -gg(i-1)*preco
35			\| ( \| \|
36			\| ) \| >>> if jmax > 0 <<<
37			\| ( \| dd(i-1) = -gg(i-1)
38			\| ) \| CALL HESSUPD
39			\| ( \| \|
40			\| ) \| \|---- dd(i-1) modified via Hessian approx.
41			\| ( \|
42			\| ) \|---- >>> if <dd,gg> >= 0 <<<
43			\| ( \| ifail = 4
44			\| ) \|
45			\| ( \|---- compute step size: tact(i-1)
46			\| ) \|---- compute update: xdiff(i) = xx(i-1) + tact(i-1)*dd(i-1)
47			\| (
48			\| )---- >>> if ifail = 4 <<<
49			\| ( goto 1000
50			\| )
51			\| (---- CALL OPTLINE / LSLINE
52			\| ) \|
53			... ... ...
54			\end{verbatim}
55			}