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
|
#========================================
# AUTO demo sspg
# BPCONT for non-generic non-periodic BVP
#========================================
pgm="sspg"
print "\n***"+pgm+": run 1: trivial solution***"
r1=run(pgm,sv=pgm+'.1')
bp=r1('BP1')
print "\n***"+pgm+": run 1: branch switch***"
run(bp,ISW=-1,NMX=5,NPR=1,sv=pgm+'.1.sw')
print "\n***"+pgm+": run 2: switch to BP continuation***"
sw=run(bp,ICP=[1,2],ISW=2,NMX=50,NPR=10,sv=pgm+'.2.sw')
print "\n***"+pgm+": run 2: BP continuation (fw)***"
run(sw,sv=pgm+'.2.fw')
print "\n***"+pgm+": run 2: BP continuation (bw)***"
run(sw,DS='-',sv=pgm+'.2.bw')
print "\n***Clean the directory***"
cl()
|
