#*********************************************************************** # # use with test problem Square # # symmetry breaking bifurcation problem # #*********************************************************************** # # initialize coefficient functions and create mesh # u a1x=.1, a2y=.1, cu1=-1, cu0=0, icont=13, left=1, top=1, right=1, bottom=1 t iadapt=4, irefn=4 # # follow main branch to first bifurcation point # s ispd=1, rltrgt=0, mxcg=50, dtol=0 s rltrgt=.3 s itask=1, rtrgt=2 s rtrgt=2 s rtrgt=4 s rtrgt=9 s rtrgt=10, itask=1 w rwfile=file1.rw # # follow main branch past bifurcation point # s rtrgt=12 s rtrgt=20 # # return to bifurcation point, follow secondary branch # r rwfile=file1.rw s itask=2 s rtrgt=8.9 s rtrgt=8.8 s rtrgt=8.5 f lines=2 s rtrgt=8.1 s rtrgt=8.05 # # locate the second bifurcation point # s rtrgt=8, itask=1 w rwfile=file2.rw # # follow tertiary branch # s itask=2 s rtrgt=8.4 f s rtrgt=9 s rtrgt=10 s rtrgt=11 s rltrgt=.09 s rltrgt=.085 s rltrgt=.08 s rtrgt=10 s rtrgt=9 s rltrgt=.035 s rltrgt=.03 s rtrgt=10 s rtrgt=12 s rtrgt=15 s rtrgt=20 f # # return to second bifurcation point, follow secondary branch # r rwfile=file2.rw s rtrgt=8 s rtrgt=7.5 s rtrgt=7, itask=1 s rtrgt=6 s rtrgt=5 s rltrgt=.1 s rltrgt=.08 s rltrgt=.06 s rltrgt=.04 s rtrgt=5 s rtrgt=6 s rtrgt=8 s rtrgt=12 s rtrgt=20 f g igrsw=3 #q