#*********************************************************************** # # use with test problem Square # # demonstrate multigraph solver for a variety of pdes # # use three different uniform meshes and linear elements # #*********************************************************************** # S dtol=1.e-3,method=1 # # first pass # t iadapt=4, irefn=25 s ifirst=1, mxnwtt=1, mxcg=200, ispd=1 g igrsw=1, t="-Delta u = 1" f ifun=0, nx=1, ny=1, lines=1 # u bx=1000 s ifirst=1, mxnwtt=1, mxcg=200, ispd=0 g igrsw=1, t="-Delta u + 1000 u_x = 1" f ifun=0, nx=1, ny=1, lines=1 # u by=1000 s ifirst=1, mxnwtt=1, mxcg=200, ispd=0 g igrsw=1, t="-Delta u + 1000 u_x + 1000 u_y = 1" f ifun=0, nx=1, ny=1, lines=1 # u bx=0, by=0, cu1=1000 s ifirst=1, mxnwtt=1, mxcg=200, ispd=1 g igrsw=1, t="-Delta u - 1000 u = 1" f ifun=0, nx=1, ny=1, lines=1 # u cu1=-1000 s ifirst=1, mxnwtt=1, mxcg=200, ispd=1 g igrsw=1, t="-Delta u + 1000 u = 1" f ifun=0, nx=1, ny=1, lines=1 # u cu1=0, x1=1.e-3 s ifirst=1, mxnwtt=1, mxcg=200, ispd=1 g igrsw=1, t="-.001 u_xx - u_yy = 1" f ifun=0, nx=1, ny=1, lines=1 # u x1=1, cir=1000 s ifirst=1, mxnwtt=1, mxcg=200, ispd=0 g igrsw=1, t="-Delta u - 1000 ( (y-1/2) u_x - (x-1/2) u_y ) = 1" f ifun=0, nx=1, ny=1, lines=1 # u cir=0, left=0, right=0, f0=1 t iadapt=4, irefn=25 s ifirst=1, mxnwtt=1, mxcg=200, ispd=1 g igrsw=1, t="-Delta u = 1 + x - y (periodic bc)" f ifun=0, nx=1, ny=1, lines=1 # # second pass # u left=2, right=2, f0=0 t iadapt=4, irefn=50 s ifirst=1, mxnwtt=1, mxcg=200, ispd=1 g igrsw=1, t="-Delta u = 1" f ifun=0, nx=1, ny=1, lines=1 # u bx=1000 s ifirst=1, mxnwtt=1, mxcg=200, ispd=0 g igrsw=1, t="-Delta u + 1000 u_x = 1" f ifun=0, nx=1, ny=1, lines=1 # u by=1000 s ifirst=1, mxnwtt=1, mxcg=200, ispd=0 g igrsw=1, t="-Delta u + 1000 u_x + 1000 u_y = 1" f ifun=0, nx=1, ny=1, lines=1 # u bx=0, by=0, cu1=1000 s ifirst=1, mxnwtt=1, mxcg=200, ispd=1 g igrsw=1, t="-Delta u - 1000 u = 1" f ifun=0, nx=1, ny=1, lines=1 # u cu1=-1000 s ifirst=1, mxnwtt=1, mxcg=200, ispd=1 g igrsw=1, t="-Delta u + 1000 u = 1" f ifun=0, nx=1, ny=1, lines=1 # u cu1=0, x1=1.e-3 s ifirst=1, mxnwtt=1, mxcg=200, ispd=1 g igrsw=1, t="-.001 u_xx - u_yy = 1" f ifun=0, nx=1, ny=1, lines=1 # u x1=1, cir=1000 s ifirst=1, mxnwtt=1, mxcg=200, ispd=0 g igrsw=1, t="-Delta u - 1000 ( (y-1/2) u_x - (x-1/2) u_y ) = 1" f ifun=0, nx=1, ny=1, lines=1 # u cir=0, left=0, right=0, f0=1 t iadapt=4, irefn=50 s ifirst=1, mxnwtt=1, mxcg=200, ispd=1 g igrsw=1, t="-Delta u = 1 + x - y (periodic bc)" f ifun=0, nx=1, ny=1, lines=1 # # third pass # u left=2, right=2, f0=0 t iadapt=4, irefn=100 s ifirst=1, mxnwtt=1, mxcg=200, ispd=1 g igrsw=1, t="-Delta u = 1" f ifun=0, nx=1, ny=1, lines=1 # u bx=1000 s ifirst=1, mxnwtt=1, mxcg=200, ispd=0 g igrsw=1, t="-Delta u + 1000 u_x = 1" f ifun=0, nx=1, ny=1, lines=1 # u by=1000 s ifirst=1, mxnwtt=1, mxcg=200, ispd=0 g igrsw=1, t="-Delta u + 1000 u_x + 1000 u_y = 1" f ifun=0, nx=1, ny=1, lines=1 # u bx=0, by=0, cu1=1000 s ifirst=1, mxnwtt=1, mxcg=200, ispd=1 g igrsw=1, t="-Delta u - 1000 u = 1" f ifun=0, nx=1, ny=1, lines=1 # u cu1=-1000 s ifirst=1, mxnwtt=1, mxcg=200, ispd=1 g igrsw=1, t="-Delta u + 1000 u = 1" f ifun=0, nx=1, ny=1, lines=1 # u cu1=0, x1=1.e-3 s ifirst=1, mxnwtt=1, mxcg=200, ispd=1 g igrsw=1, t="-.001 u_xx - u_yy = 1" f ifun=0, nx=1, ny=1, lines=1 # u x1=1, cir=1000 s ifirst=1, mxnwtt=1, mxcg=200, ispd=0 g igrsw=1, t="-Delta u - 1000 ( (y-1/2) u_x - (x-1/2) u_y ) = 1" f ifun=0, nx=1, ny=1, lines=1 # u cir=0, left=0, right=0, f0=1 t iadapt=4, irefn=100 s ifirst=1, mxnwtt=1, mxcg=200, ispd=1 g igrsw=1, t="-Delta u = 1 + x - y (periodic bc)" f ifun=0, nx=1, ny=1, lines=1 # q