#*********************************************************************** # # use this journal file with test problem Square # # nonlinear eigenvalue problem # #*********************************************************************** # # initialize coefficient functions and create mesh # Smethod=1,dtol=0,mxnwtt=10 u c0=0, icont=15 t iadapt=4, irefn=8 # # follow trivial branch, locate first eigenvalue # s rltrgt=0, ispd=1, mxcg=30 s rltrgt=10 s itask=1, rltrgt=22 w # # follow bifurcating branch for first eigenvalue # s itask=2 s rtrgt=.1 s rltrgt=25 s rltrgt=50 s rltrgt=100 s rltrgt=150 s rltrgt=300 s rltrgt=500 # # refine the mesh # t iadapt=1, ndtrgt=5000 s itask=7,mxnwtt=2 t s t s i inplsw=6, lines=0 f nx=1, ny=-1, ifun=0, lines=2 # # follow trivial branch, locate second eigenvalue # r s rltrgt=30,mxnwtt=10 s rltrgt=50 s itask=1, rltrgt=60 w # # follow bifurcating branch for second eigenvalue # s itask=2 s rtrgt=.1 s rltrgt=60 s rltrgt=110 s rltrgt=150 s rltrgt=300 s rltrgt=500 # # refine the mesh # t iadapt=1, ndtrgt=5000 s itask=7,mxnwtt=2 t s t s i inplsw=6, lines=0 f nx=1, ny=-1, ifun=0, lines=2 # # follow trivial branch, locate third eigenvalue # r s rltrgt=70,mxnwtt=10 s rltrgt=80 s itask=1, rltrgt=85 w # # follow bifurcating branch for third eigenvalue # s itask=2 s rtrgt=.1 s rltrgt=100 s rltrgt=150 s rltrgt=200 s rltrgt=300 s rltrgt=400 s rltrgt=500 # # refine the mesh # t iadapt=1, ndtrgt=5000 s itask=7,mxnwtt=2 t s t s i inplsw=6, lines=0 f nx=1, ny=-1, ifun=0, lines=2 # # follow trivial branch, locate fourth eigenvalue # r s rltrgt=100,mxnwtt=10 s rltrgt=101 s itask=1, rltrgt=103 w # # follow bifurcating branch for fourth eigenvalue # s itask=2 s rtrgt=.1 s rltrgt=110 s rltrgt=125 s rltrgt=150 s rltrgt=200 s rltrgt=300 s rltrgt=400 s rltrgt=500 # # refine the mesh # t iadapt=1, ndtrgt=5000 s itask=7,mxnwtt=2 t s t s i inplsw=6, lines=0 f nx=1, ny=-1, ifun=0, lines=2 # g igrsw=3 q