diff options
Diffstat (limited to 'CarpetExtra/SpaceTimeToy/par/spacetimetoy_periodic.gnu')
| -rw-r--r-- | CarpetExtra/SpaceTimeToy/par/spacetimetoy_periodic.gnu | 30 |
1 files changed, 30 insertions, 0 deletions
diff --git a/CarpetExtra/SpaceTimeToy/par/spacetimetoy_periodic.gnu b/CarpetExtra/SpaceTimeToy/par/spacetimetoy_periodic.gnu new file mode 100644 index 000000000..b78bdf751 --- /dev/null +++ b/CarpetExtra/SpaceTimeToy/par/spacetimetoy_periodic.gnu @@ -0,0 +1,30 @@ +# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/SpaceTimeToy/par/spacetimetoy_periodic.gnu,v 1.1 2003/06/18 18:24:29 schnetter Exp $ + +set grid +set style data linespoints + +kx=0.222222222222222 +ky=0.151515151515152 +kz=0.222222222222222 +omega=sqrt(kx**2+ky**2+kz**2) +dt=0.3 +phi(x)=cos(pi*(kx*x+omega*dt*t)) +psi(x)=-pi*omega*sin(pi*(kx*x+omega*dt*t)) + +t=0 + + + +# unigrid function +p [-4.5:4.5] "spacetimetoy_periodic_coarse/phi.xl" i t u 10:13, "spacetimetoy_periodic/phi.xl" i 2*t u 10:13, "spacetimetoy_periodic_fine/phi.xl" i 4*t u 10:13, phi(x) + +# unigrid error +p [-4.5:4.5] "spacetimetoy_periodic_coarse/phi.xl" i t u 10:($13-phi($10)), "spacetimetoy_periodic/phi.xl" i 2*t u 10:(4*($13-phi($10))), "spacetimetoy_periodic_fine/phi.xl" i 4*t u 10:(16*($13-phi($10))) + + + +# refinement function +p [-4.5:4.5] "spacetimetoy_periodic_coarse_rl2/phi.xl" i 3*t u 10:13, "spacetimetoy_periodic_rl2/phi.xl" i 6*t u 10:13, "spacetimetoy_periodic_fine_rl2/phi.xl" i 12*t u 10:13, phi(x) + +# refinement error +p [-4.5:4.5] "spacetimetoy_periodic_coarse_rl2/phi.xl" i 3*t u 10:($13-phi($10)), "spacetimetoy_periodic_rl2/phi.xl" i 6*t u 10:(4*($13-phi($10))), "spacetimetoy_periodic_fine_rl2/phi.xl" i 12*t u 10:(16*($13-phi($10))) |