blob: 13c569350a394b23c04b23409740d5512e23c19c (
plain)
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
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
|
\subsection{The \texttt{FPSPH:} module}
The {\tt FPSPH:} module performs a single SPH iteration. This module is intended to be used in a {\tt REPEAT UNTIL}
loop, implemented in CLE-2000 macro-language, where the macro calculation is explicitely called, as depicted in
Fig.~\ref{fig:fig_fpsph}.
\vskip 0.2cm
\begin{figure}[h!]
\begin{center}
\includegraphics[scale=0.85]{fpsph.eps}
\caption{Fixed point SPH iterations.}\label{fig:fig_fpsph}
\end{center}
\end{figure}
The calling specifications are:
\begin{DataStructure}{Structure \moc{FPSPH:}}
\dusa{OPTIM} \moc{:=} \moc{FPSPH:} $[$ \dusa{OPTIM} $]$ \dusa{MACROLIB} \dusa{MACROREF} \moc{::} \dstr{fpsph\_data}
\end{DataStructure}
\noindent where
\begin{ListeDeDescription}{mmmmmmmm}
\item[\dusa{OPTIM}] \texttt{character*12} name of the \dds{optimize} object ({\tt L\_OPTIMIZE} signature) containing the
SPH factors. At the first call, object \dusa{OPTIM} must appear on LHS to receive its initial values. On subsequent calls, object
\dusa{OPTIM} must appear on both LHS and RHS to be able to update the previous values.
\item[\dusa{MACROLIB}] \texttt{character*12} name of the read-only extended \dds{macrolib} object ({\tt L\_MACROLIB} signature) containing the
macroscopic cross sections used by the macro-calculation and fluxes produced by the macro-calculation.
\item[\dusa{MACROREF}] \texttt{character*12} name of the read-only extended \dds{macrolib} object ({\tt L\_MACROLIB} signature) containing the
reference macroscopic cross sections and fluxes.
\item[\dstr{fpsph\_data}] structure containing the data to the module \texttt{FPSPH:} (see Sect.~\ref{sect:lnsr_data}).
\end{ListeDeDescription}
\clearpage
\subsubsection{Data input for module \texttt{FPSPH:}}\label{sect:lnsr_data}
\begin{DataStructure}{Structure \moc{fpsph\_data}}
$[$ \moc{EDIT} \dusa{iprint} $]$ \\
$[~$\moc{SPH} $\{$ \moc{PN} $|$ \moc{SN} $\}~]$ \\
$[$ \moc{GRPMIN} \dusa{ngr1} $]~[$ \moc{GRPMAX} \dusa{ngr2} $]$\\
$[$ \moc{OUT-STEP-EPS} \dusa{$\epsilon_{ext}$} \\
$[$ \moc{VAR-VAL-MIN} \dusa{varmin} $]~[$ \moc{VAR-VAL-MAX} \dusa{varmax} $]$ \\
$[$ \moc{OUT-CONV-TST} {\tt >>} \dusa{$l_{conv}$} {\tt <<} $[$ {\tt >>} \dusa{$rms_{conv}$} {\tt <<} $]~]$ \\
;
\end{DataStructure}
\noindent where
\begin{ListeDeDescription}{mmmmmmmm}
\item[\moc{EDIT}] keyword used to set \dusa{iprint}.
\item[\dusa{iprint}] index used to control the printing in module.
\item[\moc{PN}] keyword to activate a calculation of heterogeneous SPH factors of diffusion, PN or SPN type.
\item[\moc{SN}] keyword to activate a calculation of heterogeneous SPH factors of PIJ, IC, SN or MOC type.
This is the default option.
\item[\moc{GRPMIN}] keyword used to set the first energy group where SPH correction is applied. By default,
the first energy group index is used.
\item[\dusa{ngr1}] minimum energy group index where SPH correction is applied.
\item[\moc{GRPMAX}] keyword used to set the last energy group where SPH correction is applied. By default,
the total number of energy groups in \dusa{MACROLIB} is used.
\item[\dusa{ngr2}] maximum energy group index where SPH correction is applied.
\item[\moc{OUT-STEP-EPS}] keyword used to set the tolerance of SPH iteration convergence inside module {\tt FPSPH:} (default value
is $1.0 \times 10^{-4}$).
\item[\dusa{$\epsilon_{ext}$}] tolerance value (real or double precision).
\item[\moc{VAR-VAL-MIN}] keyword to specify the minimum values of the SPH factors. These values can also be set in a previous call
to module {\tt GRAD:}.
\item[\dusa{varmin}] single real value used for all SPH factors.
\item[\moc{VAR-VAL-MAX}] keyword to specify the maximum values of the SPH factors. These values can also be set in a previous call
to module {\tt GRAD:}.
\item[\dusa{varmax}] single real value used for all SPH factors.
\item[\moc{OUT-CONV-TST}] keyword used to determine if the SPH convergence has been reached.
\item[\dusa{$l_{conv}$}] $=1$ means that SPH convergence has been reached; $=0$ otherwise.
\item[\dusa{$rms_{conv}$}] RMS error of the converged SPH factors.
\end{ListeDeDescription}
\clearpage
|
