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
|
/*
* GENERIC.c
*/
#include <model_headers.h>
/* ------------------------------------------------------------------------
function used to define the vector field or map
------------------------------------------------------------------------ */
int user_ds_func(f,x,p)
double *f,*x,*p;
{
}
/* ------------------------------------------------------------------------
function used to define the Jacobian
------------------------------------------------------------------------ */
/*
int user_jac(m,x,p)
double **m, *x, *p;
{
}
*/
/* ------------------------------------------------------------------------
function used to define inverse or approximate inverse
------------------------------------------------------------------------ */
/*
int user_inv(y,x,p)
double *y,*x,*p;
{
}
*/
/* ------------------------------------------------------------------------
function used to define aux functions of the varbs, time, or params
------------------------------------------------------------------------ */
/*
int user_aux_func(f,x,p)
double *f,*x,*p;
{
}
*/
/* ------------------------------------------------------------------------
Procedure to define default data for the dynamical system. NOTE: You
may change the entries for each variable but DO NOT change the list of
items. If a variable is unused, NULL or zero the entry, as appropriate.
------------------------------------------------------------------------ */
int user_init()
{
/* ------------ define the dynamical system in this segment --------------- */
int n_varb=2; /* dim of phase space */
static char *variable_names[]={"x","y"}; /* list of phase varb names */
static double variables[]={0.,0.}; /* default varb initial values */
static double variable_min[]={0.,0.}; /* default varb min for display */
static double variable_max[]={1.,1.}; /* default varb max for display */
static char *indep_varb_name="time"; /* name of indep variable */
static double indep_varb_min=0.; /* default indep varb min for display */
static double indep_varb_max=10000.; /* default indep varb max for display */
int n_param=2; /* dim of parameter space */
static char *parameter_names[]={"a","b"}; /* list of param names */
static double parameters[]={0.,0.}; /* initial parameter values */
static double parameter_min[]={0.,0.}; /* default param min for display */
static double parameter_max[]={1.,1.}; /* default param max for display */
int n_funct=2; /* number of user-defined functions */
static char *funct_names[]={"f1", "f2"}; /* list of funct names; {""} if none */
static double funct_min[]={0.,0.}; /* default funct min for display */
static double funct_max[]={1.,1.}; /* default funct max for display */
int manifold_type=EUCLIDEAN; /* PERIODIC (a periodic varb) or EUCLIDEAN */
static int periodic_varb[]={FALSE, FALSE}; /* if PERIODIC, which varbs are periodic? */
static double period_start[]={0.,0.}; /* if PERIODIC, begin fundamental domain */
static double period_end[]={1., 1.}; /* if PERIODIC, end of fundamental domain */
int mapping_toggle=TRUE; /* this is a map? TRUE or FALSE */
int inverse_toggle=FALSE; /* if so, is inverse FALSE, APPROX_INV, */
/* or EXPLICIT_INV? FALSE for vec field */
/* In this section, input NULL or the name of the function which contains... */
int (*def_name)()=NULL; /* the eqns of motion */
int (*jac_name)()=NULL; /* the jacobian (deriv w.r.t. space) */
int (*aux_func_name)()=NULL; /* the auxiliary functions */
int (*inv_name)()=NULL; /* the inverse or approx inverse */
int (*dfdt_name)()=NULL; /* the deriv w.r.t time */
int (*dfdparam_name)()=NULL; /* the derivs w.r.t. parameters */
c_filename = __FILE__; /* display this file for help with this system */
/* ------------------ end of dynamical system definition ------------------ */
#include <ds_define.c>
return 0;
}
|
