File: demo2.auto

package info (click to toggle)
auto-07p 0.9.1%2Bdfsg-7
  • links: PTS, VCS
  • area: main
  • in suites: buster
  • size: 16,200 kB
  • sloc: fortran: 22,644; f90: 19,340; python: 19,045; ansic: 11,116; sh: 1,079; makefile: 618; perl: 339
file content (20 lines) | stat: -rwxr-xr-x 434 bytes parent folder | download | duplicates (5) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
 
demo('bvp')

bvp = run('bvp')
branchpoints = bvp("BP")
for solution in branchpoints:
    bp = load(solution, ISW=-1, NTST=50)
    # Compute forwards
    print "Solution label", bp["LAB"], "forwards"
    fw = run(bp)
    # Compute backwards
    print "Solution label", bp["LAB"], "backwards"
    bw = run(bp,DS='-')
    both = fw + bw
    merged = merge(both)
    bvp = bvp + merged

bvp=relabel(bvp)
save(bvp, 'bvp')
plot(bvp)
wait()