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;
|
