File: test_continuation.param

package info (click to toggle)
getfem 5.4.4%2Bdfsg1-5
  • links: PTS, VCS
  • area: main
  • in suites: sid, trixie
  • size: 31,640 kB
  • sloc: cpp: 126,151; ansic: 24,798; python: 9,244; sh: 3,648; perl: 1,829; makefile: 1,367
file content (56 lines) | stat: -rw-r--r-- 2,737 bytes parent folder | download | duplicates (2)
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
% Copyright (C) 2017-2020 Yves Renard.
%
% This file is a part of GetFEM++
%
% GetFEM++  is  free software;  you  can  redistribute  it  and/or modify it
% under  the  terms  of the  GNU  Lesser General Public License as published
% by  the  Free Software Foundation;  either version 3 of the License,  or
% (at your option) any later version along with the GCC Runtime Library
% Exception either version 3.1 or (at your option) any later version.
% This program  is  distributed  in  the  hope  that it will be useful,  but
% WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
% or  FITNESS  FOR  A PARTICULAR PURPOSE.  See the GNU Lesser General Public
% License and GCC Runtime Library Exception for more details.
% You  should  have received a copy of the GNU Lesser General Public License
% along  with  this program;  if not, write to the Free Software Foundation,
% Inc., 51 Franklin St, Fifth Floor, Boston, MA  02110-1301, USA.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% parameters for program test_continuation                                %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%   discretisation parameters:	                     	      %%%%%
NX = 10;            	       % space step
FEM_TYPE = 'FEM_PK(1,1)';      % P1 for segments
INTEGRATION = 'IM_GAUSS1D(3)'; % Gauss-Legendre integration on the
                               % segment of order 3

%%%%%   continuation data:		                     	      %%%%%
DATAPATH = 'data/';
%BP_ROOTFILENAME = 'continuation_step_62_bp';	% root of the file with 
					% bifurcation data to be loaded
IND_BRANCH = 2;		% index of the branch to be continued when starting 
			% from a bifurcation point
DIRECTION = 1.;		% direction of the initial tangent
LAMBDA0 = 0.;		% initial value of parameter
NBSTEP = 80;		% number of continuation steps

SINGULARITIES = 2;	% 0: ignoring singular points
			% 1: automatic detection of limit points
			% 2: even detection and treatment of bifurcation pts 
H_INIT = 2E-2;		% initial step length
H_MAX = 2E-1;		% maximal step length
H_MIN = 2E-5;		% minimal step length
H_INC = 1.3;		% scale factor for increasing the step length
H_DEC = 0.5;		% scale factor for decreasing the step length

MAXITER = 5;		% maximum iterations of the Newton correction
THR_ITER = 4;		% threshold for the number of iterations for
			% increasing the step length
RESIDUAL = 1E-6;     	% residual
DIFFERENCE = 1E-6;	% difference of two forthcoming iteratives
COS = 0.997; 		% cosine of the angle between tangents at
			% two forthcoming points
RESIDUAL_SOLVE = 1E-8;	% initial residual

%%%%%   output parameters     		                             %%%%%
NOISY = 1;