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
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
|
<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN" "http://www.w3.org/TR/html401/loose.dtd">
<html>
<!-- Copyright (C) 2004, 2005, 2006, 2007 Alain Lahellec
Copyright (C) 2004, 2005, 2006, 2007 Patrice Dumas
Copyright (C) 2004, Ste'phane Hallegatte
Permission is granted to copy, distribute and/or modify this document
under the terms of the GNU Free Documentation License, Version 1.1 or
any later version published by the Free Software Foundation; with no
Invariant Sections, with no Front-Cover text and with no Back-Cover Text.
A copy of the license is included in the section entitled "GNU Free
Documentation License."
-->
<!-- Created on a sunny day by texi2html
texi2html was written by:
Lionel Cons <Lionel.Cons@cern.ch> (original author)
Karl Berry <karl@freefriends.org>
Olaf Bachmann <obachman@mathematik.uni-kl.de>
and many others.
Maintained by: Many creative people.
Send bugs and suggestions to <texi2html-bug@nongnu.org>
-->
<head>
<title>Miniker 102 manual: 3.9 Entering model size explicitely</title>
<meta name="description" content="Miniker 102 manual: 3.9 Entering model size explicitely">
<meta name="keywords" content="Miniker 102 manual: 3.9 Entering model size explicitely">
<meta name="resource-type" content="document">
<meta name="distribution" content="global">
<meta name="Generator" content="texi2html">
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<style type="text/css">
<!--
@import "mini_ker_tex4ht_math.css";
@import "mini_ker_tex4ht_tex.css";
a.summary-letter {text-decoration: none}
blockquote.smallquotation {font-size: smaller}
pre.display {font-family: serif}
pre.format {font-family: serif}
pre.menu-comment {font-family: serif}
pre.menu-preformatted {font-family: serif}
pre.smalldisplay {font-family: serif; font-size: smaller}
pre.smallexample {font-size: smaller}
pre.smallformat {font-family: serif; font-size: smaller}
pre.smalllisp {font-size: smaller}
span.roman {font-family:serif; font-weight:normal;}
span.sansserif {font-family:sans-serif; font-weight:normal;}
ul.toc {list-style: none}
-->
</style>
</head>
<body lang="en" bgcolor="#FFFFFF" text="#000000" link="#0000FF" vlink="#800080" alink="#FF0000">
<a name="Explicit-model-size"></a>
<ul class="toc"><li><a href="mini_ker.html#Top">Miniker 102 manual</a> </li>
<li><ul class="toc"><li><a href="Advanced-programming.html#Advanced-programming">3. Advanced Miniker programming</a> </li>
</ul></li>
</ul>
<table cellpadding="1" cellspacing="1" border="0">
<tr><td valign="middle" align="left">[<a href="Observations-and-data.html#Data" title="Previous section in reading order"> < </a>]</td>
<td valign="middle" align="left">[<a href="#Size-sequence" title="Next section in reading order"> > </a>]</td>
<td valign="middle" align="left"> </td>
<td valign="middle" align="left">[<a href="mini_ker.html#SEC_Contents" title="Table of contents">Contents</a>]</td>
<td valign="middle" align="left">[<a href="Concepts-index.html#Concepts-index" title="Index">Index</a>]</td>
<td valign="middle" align="left">[<a href="mini_ker_abt.html#SEC_About" title="About (help)"> ? </a>]</td>
</tr></table>
<a name="Entering-model-size-explicitely"></a>
<h2 class="section"> 3.9 Entering model size explicitely </h2>
<ul class="toc">
<li> <a href="#Size-sequence">3.9.1 The explicit size sequence</a> </li>
<li> <a href="#Model-with-explicit-size">3.9.2 Entering the model equations, with explicit sizes</a> </li>
</ul>
<p>It is possible to enter the model dimensions explicitely, instead of
generating them automatically, as it was done previously.
This feature is turned on by <code>sel dimetaphi</code>
in ‘<tt>selseq.kumac</tt>’ with cmz
and <code>dimetaphi</code> added to the <code>SEL</code> variable in
the ‘<tt>Makefile</tt>’ with make.
</p>
<hr size="2">
<a name="Size-sequence"></a>
<table cellpadding="1" cellspacing="1" border="0">
<tr><td valign="middle" align="left">[<a href="#Explicit-model-size" title="Previous section in reading order"> < </a>]</td>
<td valign="middle" align="left">[<a href="#Model-with-explicit-size" title="Next section in reading order"> > </a>]</td>
<td valign="middle" align="left"> </td>
<td valign="middle" align="left">[<a href="mini_ker.html#SEC_Contents" title="Table of contents">Contents</a>]</td>
<td valign="middle" align="left">[<a href="Concepts-index.html#Concepts-index" title="Index">Index</a>]</td>
<td valign="middle" align="left">[<a href="mini_ker_abt.html#SEC_About" title="About (help)"> ? </a>]</td>
</tr></table>
<a name="The-explicit-size-sequence"></a>
<h3 class="subsection"> 3.9.1 The explicit size sequence </h3>
<a name="index-dimetaphi"></a>
<a name="index-model-size"></a>
<a name="index-np-1"></a>
<a name="index-mp-1"></a>
<a name="index-maxstep"></a>
<a name="index-dimetaphi-1"></a>
<p>The dimension of the model is entered in the sequence ‘<tt>dimetaphi</tt>’,
using the fortran <code>parameter np</code> for <code>eta(.)</code> and
<code>mp</code> for <code>ff(.)</code>.
For the Lotka-Volterra model, we have two cell components and only one transfer.
</p>
<table><tr><td> </td><td><pre class="example">parameter (np=2,mp=1);
</pre></td></tr></table>
<p>You should not change the layout of the parameter statement as the
mortran preprocessor matches the line.
</p>
<p>You also have to provide other parameters even if you don’t have any
use for them. If you don’t it will trigger fortran errors.
It includes the <code>maxstep</code> parameter that can have any value but 0,
<code>lp</code> and <code>mobs</code> that should be 0 in the example, and <code>nxp</code>,
<code>nyp</code> and <code>nzp</code> that should also be 0.
The layout is the following:
</p>
<table><tr><td> </td><td><pre class="example">parameter (np=2,mp=1);
parameter (mobs=0);
parameter (nxp=0,nyp=0,nzp=0);
parameter (lp=0);
parameter (maxstep=1);
</pre></td></tr></table>
<p>If there are observations, (see <a href="Observations-and-data.html#Observations">Observations</a>), the
size of the observation vector is set in the ‘<tt>dimetaphi</tt>’ sequence
by the <code>mobs</code> parameter. For example if there is one observation:
</p>
<table><tr><td> </td><td><pre class="example">parameter (mobs=1);
</pre></td></tr></table>
<p>To specify parameters (see <a href="Parameters.html#Parameters">Parameters</a>), the number of such parameters
has to be declared in ‘<tt>dimetaphi</tt>’ with the parameter <code>lp</code>.
Then, if there are two parameters, they are first declared with
</p>
<table><tr><td> </td><td><pre class="example">parameter (lp=2);
</pre></td></tr></table>
<hr size="2">
<a name="Model-with-explicit-size"></a>
<table cellpadding="1" cellspacing="1" border="0">
<tr><td valign="middle" align="left">[<a href="#Size-sequence" title="Previous section in reading order"> < </a>]</td>
<td valign="middle" align="left">[<a href="Programming-with-cmz-directives.html#Programming-with-cmz-directives" title="Next section in reading order"> > </a>]</td>
<td valign="middle" align="left"> </td>
<td valign="middle" align="left">[<a href="mini_ker.html#SEC_Contents" title="Table of contents">Contents</a>]</td>
<td valign="middle" align="left">[<a href="Concepts-index.html#Concepts-index" title="Index">Index</a>]</td>
<td valign="middle" align="left">[<a href="mini_ker_abt.html#SEC_About" title="About (help)"> ? </a>]</td>
</tr></table>
<a name="Entering-the-model-equations_002c-with-explicit-sizes"></a>
<h3 class="subsection"> 3.9.2 Entering the model equations, with explicit sizes </h3>
<a name="index-model-equations"></a>
<a name="index-Phi_005ftef_0028_002e_0029"></a>
<a name="index-deta_005ftef_0028_002e_0029"></a>
<a name="index-eta_0028_002e_0029_002c-explicit-sizes"></a>
<a name="index-ff_0028_002e_0029_002c-explicit-sizes"></a>
<p>When sizes are explicit, another possibility exists for entering
the model equations. The use of symbolic names, as described in
<a href="A-model-description.html#Model-equations">Model equations</a> is still possible, and it also becomes possible to
set directly the equations associated with the <code>eta(.)</code>
and <code>ff(.)</code> vectors.
</p>
<p>In case the symbolic names are not used,
the model equations for cells and transfers are entered using a mortran macro,
<code>f_set</code><a name="DOCF4" href="mini_ker_fot.html#FOOT4">(4)</a>, setting the <code>eta(.)</code> evolution with
<code>deta_tef(.)</code>
and the transfer definitions <code>ff(.)</code> with <code>Phi_tef(.)</code>.
</p>
<dl>
<dt><a name="index-f_005fset"></a><u>Macro:</u> <b>f_set</b><i> Phi_tef(<var>i</var>) = f(eta(.),ff(.))</i></dt>
<dd><p>This macro defines the transfer <var>i</var> static equation.
<code>f</code> is a fortran
expression which may be function of cell state variables,
‘<samp>eta(1)</samp>’…‘<samp>eta(np)</samp>’ and transfers
‘<samp>ff(1)</samp>’…‘<samp>ff(mp)</samp>’.
</p></dd></dl>
<p>In the case of the predator-prey model, the transfer definition for
φmeet is:
</p><table><tr><td> </td><td><pre class="example">f_set Phi_tef(1) = eta(1)*eta(2);
</pre></td></tr></table>
<dl>
<dt><a name="index-f_005fset-1"></a><u>Macro:</u> <b>f_set</b><i> deta_tef(<var>i</var>) = g(eta(<var>i</var>),ff(.))</i></dt>
<dd><p>This macro defines the cell state component <var>i</var> time evolution model.
<code>g</code> is a expression which may be function of cell state variables,
‘<samp>eta(1)</samp>’…‘<samp>eta(np)</samp>’ and transfers
‘<samp>ff(1)</samp>’…‘<samp>ff(mp)</samp>’.
</p></dd></dl>
<p>The two cell equations of the predator-prey model are, with index 1 for the
prey ( ηprey) and index 2 for the predator ( ηpred):
</p>
<table><tr><td> </td><td><pre class="example">f_set deta_tef(1) = apar*eta(1)-apar*ff(1);
f_set deta_tef(2) = - cpar*eta(2) + cpar*ff(1);
</pre></td></tr></table>
<p>The whole model is:
</p>
<table><tr><td> </td><td><pre class="example">!%%%%%%%%%%%%%%%%%%%%%%
! Transfer definition
!%%%%%%%%%%%%%%%%%%%%%%
! rencontres (meeting)
f_set Phi_tef(1) = eta(1)*eta(2);
!%%%%%%%%%%%%%%%%%%%%%%
! Cell definition
!%%%%%%%%%%%%%%%%%%%%%%
! eta(1) : prey
! eta(2) : predator
f_set deta_tef(1) = apar*eta(1)-apar*ff(1);
f_set deta_tef(2) = - cpar*eta(2) + cpar*ff(1);
</pre></td></tr></table>
<p>The starting points for cells are entered like:
</p><table><tr><td> </td><td><pre class="example">! initial state
! -------------
eta(1) = 1.;
eta(2) = 1.;
</pre></td></tr></table>
<p>If there are observations, they are entered as special transferts with
index above <code>mp</code>, for example:
</p>
<table><tr><td> </td><td><pre class="example">f_set Phi_tef(mp+1) = ff(1) ;
</pre></td></tr></table>
<hr size="2">
<table cellpadding="1" cellspacing="1" border="0">
<tr><td valign="middle" align="left">[<a href="mini_ker.html#SEC_Contents" title="Table of contents">Contents</a>]</td>
<td valign="middle" align="left">[<a href="Concepts-index.html#Concepts-index" title="Index">Index</a>]</td>
<td valign="middle" align="left">[<a href="mini_ker_abt.html#SEC_About" title="About (help)"> ? </a>]</td>
</tr></table>
<p>
<font size="-1">
This document was generated by <em>a tester</em> on <em>a sunny day</em> using <a href="http://www.nongnu.org/texi2html/"><em>texi2html</em></a>.
</font>
<br>
</p>
</body>
</html>
|
