Browse Source

moving plots to thesis

master
Stephen McQuay 9 years ago
parent
commit
0eb9c4102d
83 changed files with 0 additions and 1542 deletions
  1. +0
    -14
      plots/exact_2D.plt
  2. +0
    -28
      plots/gen_plots.zsh
  3. +0
    -8
      plots/main.plt
  4. +0
    -25
      plots/qhull_vs_delny/bench.py
  5. +0
    -19
      plots/qhull_vs_delny/output.bash.cgd
  6. +0
    -255
      plots/qhull_vs_delny/output.python.cgd
  7. +0
    -2
      plots/qhull_vs_delny/plot.gp
  8. +0
    -23
      plots/qhull_vs_delny/plot.py
  9. +0
    -99
      plots/qhull_vs_delny/results.cqd
  10. +0
    -10
      plots/resolution/resolution.2D.64.out
  11. +0
    -19
      plots/resolution/resolution.2D.regular.plt
  12. +0
    -6
      plots/resolution/resolution.3D.gmsh.out
  13. +0
    -19
      plots/resolution/resolution.3D.gmsh.plt
  14. +0
    -33
      plots/resolution/resolution_3D.py
  15. +0
    -11
      plots/rms/2D/rms.1-4.2.plt
  16. +0
    -11
      plots/rms/2D/rms.1-4.3.plt
  17. +0
    -11
      plots/rms/2D/rms.1-4.4.plt
  18. +0
    -11
      plots/rms/2D/rms.1-4.5.plt
  19. +0
    -12
      plots/rms/2D/rms.1.2-5.plt
  20. +0
    -11
      plots/rms/2D/rms.2.2-5.plt
  21. +0
    -11
      plots/rms/2D/rms.3.2-5.plt
  22. +0
    -11
      plots/rms/2D/rms.4.2-5.plt
  23. +0
    -27
      plots/rms/2D/rms.onering.plt
  24. +0
    -14
      plots/rms/2D/rms.out
  25. +0
    -10
      plots/rms/3D/rms.1-3.2.plt
  26. +0
    -10
      plots/rms/3D/rms.1-3.3.plt
  27. +0
    -10
      plots/rms/3D/rms.1-3.4.plt
  28. +0
    -10
      plots/rms/3D/rms.1-3.5.plt
  29. +0
    -13
      plots/rms/3D/rms.1.2-5.plt
  30. +0
    -13
      plots/rms/3D/rms.2.2-5.plt
  31. +0
    -11
      plots/rms/3D/rms.3.2-5.plt
  32. +0
    -22
      plots/rms/3D/rms.onering.plt
  33. +0
    -14
      plots/rms/3D/rms.out
  34. +0
    -6
      plots/rms/all.plt
  35. +0
    -12
      plots/scalability/efficiency.plt
  36. +0
    -22
      plots/scalability/hist.py
  37. +0
    -35
      plots/scalability/scalability_histograms.py
  38. +0
    -26
      plots/scalability/scalability_plot_gen.py
  39. +0
    -19
      plots/scalability/scale.out
  40. +0
    -12
      plots/scalability/server.out
  41. +0
    -12
      plots/scalability/server.plt
  42. +0
    -10
      plots/scalability/speedup.plt
  43. +0
    -27
      plots/timing/2D/onering.plt
  44. +0
    -20
      plots/timing/2D/time.1-4.2.plt
  45. +0
    -11
      plots/timing/2D/time.1-4.3.plt
  46. +0
    -11
      plots/timing/2D/time.1-4.4.plt
  47. +0
    -21
      plots/timing/2D/time.1-4.5.plt
  48. +0
    -11
      plots/timing/2D/time.1.2-5.plt
  49. +0
    -11
      plots/timing/2D/time.2.2-5.plt
  50. +0
    -11
      plots/timing/2D/time.3.2-5.plt
  51. +0
    -11
      plots/timing/2D/time.4.2-5.plt
  52. +0
    -14
      plots/timing/2D/time.out
  53. +0
    -22
      plots/timing/3D/onering.plt
  54. +0
    -16
      plots/timing/3D/time.1-3.2.plt
  55. +0
    -10
      plots/timing/3D/time.1-3.3.plt
  56. +0
    -10
      plots/timing/3D/time.1-3.4.plt
  57. +0
    -18
      plots/timing/3D/time.1-3.5.plt
  58. +0
    -11
      plots/timing/3D/time.1.2-5.plt
  59. +0
    -11
      plots/timing/3D/time.2.2-5.plt
  60. +0
    -11
      plots/timing/3D/time.3.2-5.plt
  61. +0
    -14
      plots/timing/3D/time.out
  62. +0
    -6
      plots/timing/all.plt
  63. +0
    -42
      plots/timing/timing.aggregate.plt
  64. +0
    -27
      plots/truthiness/2D/onering.plt
  65. +0
    -12
      plots/truthiness/2D/truthy.1-4.2.plt
  66. +0
    -12
      plots/truthiness/2D/truthy.1-4.3.plt
  67. +0
    -12
      plots/truthiness/2D/truthy.1-4.4.plt
  68. +0
    -12
      plots/truthiness/2D/truthy.1-4.5.plt
  69. +0
    -13
      plots/truthiness/2D/truthy.1.2-5.plt
  70. +0
    -12
      plots/truthiness/2D/truthy.2.2-5.plt
  71. +0
    -12
      plots/truthiness/2D/truthy.3.2-5.plt
  72. +0
    -12
      plots/truthiness/2D/truthy.4.2-5.plt
  73. +0
    -14
      plots/truthiness/2D/truthy.out
  74. +0
    -22
      plots/truthiness/3D/onering.plt
  75. +0
    -11
      plots/truthiness/3D/truthy.1-3.2.plt
  76. +0
    -11
      plots/truthiness/3D/truthy.1-3.3.plt
  77. +0
    -11
      plots/truthiness/3D/truthy.1-3.4.plt
  78. +0
    -11
      plots/truthiness/3D/truthy.1-3.5.plt
  79. +0
    -14
      plots/truthiness/3D/truthy.1.2-5.plt
  80. +0
    -14
      plots/truthiness/3D/truthy.2.2-5.plt
  81. +0
    -12
      plots/truthiness/3D/truthy.3.2-5.plt
  82. +0
    -14
      plots/truthiness/3D/truthy.out
  83. +0
    -8
      plots/truthiness/all.plt

+ 0
- 14
plots/exact_2D.plt View File

@@ -1,14 +0,0 @@
load 'main.plt'
set output 'exact.eps'

set yrange [0:1]
set xrange [0:1]

set pm3d at b
set ticslevel 0.08

set nokey

set isosample 40,40

splot (sin(x*pi) * cos(y*pi))**2

+ 0
- 28
plots/gen_plots.zsh View File

@@ -1,28 +0,0 @@
#!/bin/zsh

ROOTDIR=$HOME/src/research/plots
FIGDIR=$HOME/src/thesis/figures/results/

DIRS=(\
resolution \
scalability \
timing \
timing/2D \
timing/3D \
rms/2D \
rms/3D \
truthiness/2D \
truthiness/3D \
)

for dir in $DIRS
do
cd $dir
for pltfile in *.plt
do
gnuplot $pltfile
done
cd $ROOTDIR
done

mv $(find | grep eps$) $FIGDIR

+ 0
- 8
plots/main.plt View File

@@ -1,8 +0,0 @@
set terminal postscript enhanced eps color
#set terminal postscript enhanced eps

set grid

set key box

set size 0.8,0.8

+ 0
- 25
plots/qhull_vs_delny/bench.py View File

@@ -1,25 +0,0 @@
import sys
import time
import interp.bootstrap
import numpy as np
from interp.grid.delaunay import dgrid, get_qdelaunay_dump_str

from delaunay import Triangulation


# for count in xrange(100000,1000000, 100000):
count = int(sys.argv[1])
v = np.random.random((count, 3))
g = dgrid(v)

g_start = time.time()
g.construct_connectivity()
g_end = time.time()

dt_start = time.time()
dt = Triangulation(g.verts)
dt_end = time.time()

with open("output.cgd", 'a') as f:
s = "%d %f %f\n" % (count, g_end - g_start, dt_end - dt_start)
f.write(s)

+ 0
- 19
plots/qhull_vs_delny/output.bash.cgd View File

@@ -1,19 +0,0 @@
10000 4.358884 4.071420
20000 12.284323 13.088460
30000 23.463748 26.769070
40000 36.918351 46.537505
50000 54.396764 70.307496
60000 74.990658 98.235962
70000 98.728850 131.401528
80000 123.678278 171.596682
90000 154.050214 213.684310
100000 186.437281 259.919163
110000 221.873001 310.210802
120000 257.232807 369.233422
130000 298.404708 430.279021
140000 344.442933 499.202024
150000 391.153554 563.993615
160000 436.987274 643.019125
170000 491.391332 732.660288
180000 547.777618 819.785378
190000 611.478075 906.047083

+ 0
- 255
plots/qhull_vs_delny/output.python.cgd View File

@@ -1,255 +0,0 @@
100 0.047999 0.009422
200 0.070899 0.019659
300 0.096712 0.063634
400 0.129279 0.047459
500 0.219131 0.055515
600 0.186475 0.128657
700 0.289947 0.082739
800 0.340396 0.102652
900 0.294601 0.214491
1000 0.421002 0.131719
1100 0.461535 0.143235
1200 0.384233 0.281429
1300 0.547769 0.175322
1400 0.598529 0.291340
1500 0.471884 0.363554
1600 0.666989 0.332696
1700 0.737200 0.355773
1800 0.752624 0.372474
1900 0.813278 0.395763
2000 0.827568 0.555087
2100 0.907530 0.447374
2200 0.950769 0.468655
2300 1.122449 0.495623
2400 1.135867 0.514808
2500 1.089476 0.734838
2600 1.295460 0.570625
2700 1.182967 0.791656
2800 1.355481 0.826170
2900 1.438569 0.863660
3000 1.325295 0.885664
3100 1.581403 0.695968
3200 1.575659 0.952604
3300 1.655784 1.001539
3400 1.727168 1.045440
3500 1.778233 1.072406
3600 1.814249 1.088649
3700 1.832326 1.409561
3800 1.945704 1.165096
3900 2.008672 1.202718
4000 2.057761 1.242398
4100 2.122421 1.273546
4200 2.153430 1.650734
4300 2.516681 1.359595
4400 2.253733 1.731281
4500 2.336863 1.788446
4600 2.378671 1.834615
4700 2.437512 1.862757
4800 2.830098 1.957137
4900 2.564581 1.960565
5000 2.912312 2.004686
5100 2.983080 2.057447
5200 3.058453 2.107918
5300 3.114177 2.628001
5400 3.275035 2.242694
5500 3.197436 2.723701
5600 3.326968 2.821540
5700 3.404702 2.868773
5800 3.476629 2.910723
5900 3.625217 3.040586
6000 3.590140 3.044023
6100 3.747166 3.144620
6200 3.834644 3.203244
6300 3.886890 3.280042
6400 3.926654 3.291087
6500 4.464715 3.407888
6600 4.075471 3.413190
6700 4.536596 3.458579
6800 4.620090 3.545290
6900 4.718980 4.267050
7000 4.971424 3.725282
7100 4.956617 3.737774
7200 5.093261 3.857758
7300 5.632831 3.920162
7400 5.306636 3.984243
7500 5.336264 4.821789
7600 5.489230 5.016427
7700 5.589738 5.016706
7800 6.313402 4.288711
7900 6.415683 5.199262
8000 6.413859 5.307123
8100 6.426843 5.265939
8200 6.766083 5.445775
8300 6.638535 5.499324
8400 6.815593 5.613655
8500 6.827982 5.604333
8600 6.849010 6.532516
8700 7.184065 5.816340
8800 7.128844 6.762998
8900 7.311674 5.932302
9000 7.439522 6.090411
9100 7.473794 7.192710
9200 7.619600 7.167273
9300 8.466331 6.272211
9400 7.759528 7.293581
9500 7.960313 7.492307
9600 8.095775 7.585801
9700 8.167192 7.705659
9800 9.166040 7.856572
9900 8.457956 7.883608
10000 9.344719 8.001859
10100 9.642055 8.174025
10200 9.656403 8.260292
10300 9.656492 9.470269
10400 10.071541 8.538257
10500 9.950824 9.702932
10600 10.366404 11.537015
10700 11.578142 10.147674
10800 10.441344 10.156902
10900 11.561428 9.065712
11000 10.775458 10.983379
11100 12.049373 10.639839
11200 12.215302 11.300136
11300 12.798212 11.227852
11400 11.552312 11.166603
11500 11.395557 12.491557
11600 12.537739 11.278265
11700 12.732686 11.404658
11800 12.956403 11.619160
11900 13.111756 11.724366
12000 13.162634 13.249459
12100 13.596588 12.062135
12200 13.440376 13.535533
12300 13.629667 13.687059
12400 13.795565 13.854947
12500 15.302078 12.586318
12600 14.293713 14.468024
12700 16.968780 14.285985
12800 15.106672 15.076618
12900 14.856244 14.931597
13000 14.852538 14.938293
13100 16.584244 15.192819
13200 15.235971 15.228923
13300 16.764129 15.761469
13400 17.170068 15.750559
13500 17.573425 16.539195
13600 18.092185 16.338719
13700 17.309370 17.754095
13800 17.733436 16.266628
13900 17.968198 19.838720
14000 17.992370 18.541615
14100 18.260844 18.651974
14200 18.674546 19.330625
14300 18.681486 19.258071
14400 19.239808 19.558980
14500 20.891203 19.528348
14600 20.950823 19.768551
14700 21.176991 19.861108
14800 19.723195 19.898361
14900 21.421963 20.195574
15000 20.091257 22.622511
15100 20.168624 22.587613
15200 20.297649 22.911602
15300 20.481518 23.028209
15400 22.757843 23.436407
15500 22.908474 21.451511
15600 22.637640 23.905791
15700 23.122585 24.022241
15800 23.268845 24.239530
15900 23.842870 24.648094
16000 23.759917 24.554258
16100 24.703206 25.605368
16200 24.812305 25.546528
16300 24.860067 25.724616
16400 24.933327 25.721597
16500 27.223346 26.214275
16600 25.375054 28.650909
16700 27.556047 26.422677
16800 27.949738 26.733193
16900 27.938289 29.459037
17000 26.528751 29.850431
17100 28.714271 30.112784
17200 26.782098 30.487300
17300 29.206300 30.714398
17400 29.023970 30.660169
17500 29.728589 31.432677
17600 30.057684 31.561094
17700 32.722450 31.911768
17800 30.467706 31.943814
17900 30.730185 32.262856
18000 30.900582 35.203025
18100 33.730091 32.949535
18200 31.755282 33.193677
18300 34.427319 33.467022
18400 34.668552 33.760462
18500 35.077295 34.096138
18600 35.137330 34.233175
18700 35.200482 37.364114
18800 35.930076 35.000927
18900 33.740058 38.149483
19000 36.319830 38.407051
19100 36.448262 38.725132
19200 37.231378 39.259957
19300 37.449582 39.501841
19400 37.994335 40.048953
19500 37.708301 43.070084
19600 38.041702 43.529877
19700 41.349886 40.863797
19800 39.063304 41.045847
19900 39.019994 44.486624
20000 39.571041 45.018980
20100 40.127330 45.532098
20200 40.319021 45.691957
20300 40.498957 46.059254
20400 43.880946 46.627267
20500 41.162964 46.717860
20600 44.192133 46.899932
20700 44.627189 47.372401
20800 44.736632 47.556623
20900 45.658531 48.250996
21000 45.708075 48.446644
21100 46.495488 49.188960
21200 49.544554 49.416731
21300 47.084084 49.949468
21400 47.506258 50.251095
21500 47.476536 54.007076
21600 48.218477 54.706843
21700 51.572074 51.358229
21800 52.018751 51.748216
21900 52.172109 55.929243
22000 49.714049 56.321718
22100 53.399162 56.990122
22200 50.353217 57.105554
22300 53.054872 56.610281
22300 15.050173 15.603131
23300 23.574880 22.674618
24300 26.044030 26.281024
25300 28.112445 29.941780
26300 31.267650 31.683991
27300 32.385395 34.929903
28300 34.724820 38.467582
29300 38.499965 42.432000
30300 40.724985 44.726637
31300 43.037760 48.240965
32300 47.141275 51.070441
33300 49.420085 55.375932
34300 54.008640 57.483027
35300 54.893517 61.485955
36300 59.278087 66.502684
37300 61.418032 72.472749
38300 65.521988 74.227871
39300 68.432421 79.789643
40300 72.580188 82.202839
41300 75.431865 88.730005
42300 82.076166 89.980316
43300 84.260960 96.160103
44300 89.615419 103.184661
45000 45.710591 58.462027
55000 97.647944 117.987362
65000 140.092906 163.141031
75000 187.001122 228.359633
85000 253.616062 294.059444
95000 299.802204 378.583525
100000 206.991352 254.952767
150000 398.023458 597.453935
200000 969.711886 1400.735031

+ 0
- 2
plots/qhull_vs_delny/plot.gp View File

@@ -1,2 +0,0 @@
plot "./output.cgd-1" using 1:2 with lines, "./output.cgd-1" using 1:3 with lines, "./output.cgd" using 1:2 with lines, "./output.cgd" using 1:3 with lines
pause -1

+ 0
- 23
plots/qhull_vs_delny/plot.py View File

@@ -1,23 +0,0 @@
import numpy as np
import matplotlib.pyplot as plt

bash = np.loadtxt('output.bash.cgd', unpack = True)
py = np.loadtxt('output.python.cgd', unpack = True)


plt.figure(1)
ax = plt.subplot(211)
ax.set_title("up")
plt.plot(bash[0], bash[1], 'b--')
plt.ylabel('time (sec)')

ax = plt.subplot(212)
ax.set_title("down")
plt.plot(py[0], py[1])


plt.xlabel('total points')
plt.ylabel('time (sec)')


plt.savefig('asdf.png')

+ 0
- 99
plots/qhull_vs_delny/results.cqd View File

@@ -1,99 +0,0 @@
100 0.042380 0.009482
200 0.070966 0.019559
300 0.099712 0.063360
400 0.129388 0.046986
500 0.211595 0.056179
600 0.246149 0.069776
700 0.287796 0.085229
800 0.338569 0.099770
900 0.280211 0.207107
1000 0.408283 0.130189
1100 0.458021 0.144868
1200 0.504123 0.159955
1300 0.548773 0.176861
1400 0.581178 0.286109
1500 0.471314 0.360158
1600 0.660644 0.323710
1700 0.714677 0.349696
1800 0.761600 0.373077
1900 0.795049 0.393191
2000 0.980112 0.426283
2100 0.884566 0.440110
2200 0.944961 0.466153
2300 0.976980 0.491693
2400 1.041750 0.516152
2500 1.078512 0.729253
2600 1.284825 0.571990
2700 1.162411 0.783929
2800 1.206563 0.816419
2900 1.399989 0.854620
3000 1.298645 0.887150
3100 1.570723 0.691291
3200 1.546390 0.941965
3300 1.672172 1.003558
3400 1.705987 1.039525
3500 1.755851 1.068202
3600 1.798685 1.083653
3700 1.807331 1.405700
3800 1.905298 1.151212
3900 1.994835 1.198307
4000 2.033566 1.235601
4100 2.094191 1.266220
4200 2.138027 1.632787
4300 2.455475 1.332195
4400 2.241637 1.698046
4500 2.315390 1.760717
4600 2.355638 1.813777
4700 2.381101 1.829883
4800 2.763281 1.918376
4900 2.534302 1.922099
5000 2.866890 1.978336
5100 2.960812 2.030787
5200 3.014747 2.082091
5300 3.016663 2.566534
5400 3.214297 2.192325
5500 3.610524 2.280557
5600 3.271708 2.774567
5700 3.331651 2.836686
5800 3.855931 2.419305
5900 3.532042 2.981868
6000 3.544821 2.999845
6100 3.700257 3.086741
6200 3.775651 3.168054
6300 3.838548 3.222814
6400 3.867673 3.287251
6500 4.549294 3.474946
6600 4.156695 4.154661
6700 4.828235 3.651414
6800 4.721340 3.573548
6900 4.614265 4.210497
7000 4.886040 3.722642
7100 4.888223 4.426519
7200 5.030514 3.814733
7300 5.159213 3.890596
7400 5.277981 3.967803
7500 5.265495 4.738340
7600 5.282582 4.789858
7700 5.469494 4.928022
7800 6.236024 5.043690
7900 6.223296 5.113952
8000 6.335423 5.201895
8100 6.340897 5.213795
8200 6.555179 5.390177
8300 6.524318 5.334415
8400 6.676610 5.439902
8500 6.773503 5.518722
8600 6.840255 6.501196
8700 7.076299 5.790069
8800 7.073919 6.743595
8900 7.309526 5.959033
9000 7.391387 5.981054
9100 7.507321 7.072105
9200 7.562241 7.145818
9300 8.576819 6.296824
9400 7.720963 7.295802
9500 7.884473 7.435793
9600 7.996540 7.552492
9700 8.096361 7.623015
9800 9.084937 7.820823
9900 8.382479 7.873313

+ 0
- 10
plots/resolution/resolution.2D.64.out View File

@@ -1,10 +0,0 @@
0.100000 1.704717e-02 1.485405e-02 1.003847e-02 5.869152e-03 3.516177e-03 1.783731e-03 1.795719e-01 5.735480e-02
0.050000 3.979204e-03 1.400396e-03 6.258123e-04 1.888244e-04 6.720196e-05 6.408175e-05 2.290350e-02 1.744233e-01
0.033333 1.627815e-03 3.174392e-04 1.210472e-04 2.024364e-05 5.079075e-06 1.720423e-05 2.406938e-01 1.619885e-02
0.025000 9.313592e-04 1.198217e-04 4.024488e-05 4.332817e-06 8.988869e-07 2.407076e-06 5.822483e-03 5.539409e-03
0.020000 5.959409e-04 5.529551e-05 1.680567e-05 1.277112e-06 2.221020e-07 9.613279e-07 1.526722e-03 1.000030e-03
0.016667 4.018160e-04 2.913706e-05 7.708216e-06 4.800731e-07 6.942728e-08 7.095225e-07 1.321824e-03 1.359354e-03
0.014286 2.884502e-04 1.655307e-05 4.189791e-06 1.830632e-07 2.721073e-08 3.431936e-07 2.271145e-04 3.885443e-04
0.012500 2.245358e-04 1.080741e-05 2.410492e-06 1.030748e-07 1.097364e-08 2.725363e-07 4.457260e-04 5.130550e-04
0.011111 1.741243e-04 7.925485e-06 1.504109e-06 5.204435e-08 5.477144e-09 1.355392e-07 2.035055e-04 3.372173e-04
0.010000 1.433374e-04 5.381351e-06 9.908304e-07 2.703056e-08 2.756173e-09 1.113693e-07 2.099056e-04 1.644788e-03

+ 0
- 19
plots/resolution/resolution.2D.regular.plt View File

@@ -1,19 +0,0 @@
load '../main.plt'

set output 'resolution.2D.eps'

# set title '2-D Meshes'

set log xy
set key right bottom

set xrange [0.008:0.18]

# set xlabel 'Mesh Spacing'
set ylabel 'RMS of error'

plot './resolution.2D.64.out' \
u 1:3 t '{/Symbol n} = 2' w l,\
'' u 1:4 t '{/Symbol n} = 3' w l,\
'' u 1:5 t '{/Symbol n} = 4' w l,\
'' u 1:6 t '{/Symbol n} = 5' w l;

+ 0
- 6
plots/resolution/resolution.3D.gmsh.out View File

@@ -1,6 +0,0 @@
# generated by running resolution_3D.py, and then by calculating x by dividing
# the length of the domain by the cube root of number of points.

0.0454730588061 0.000270108596188 8.53500404066e-05 1.41791618738e-05 4.9979436523e-06
0.00997060697355 2.0540042226e-06 1.69775361716e-07 1.01412678674e-08 8.0860277895e-10
0.00461278449044 2.18695319382e-07 1.52441703341e-08 3.24157996636e-10 2.18935788738e-11

+ 0
- 19
plots/resolution/resolution.3D.gmsh.plt View File

@@ -1,19 +0,0 @@
load '../main.plt'

set output 'resolution.3D.eps'

# set title '3-D Meshes'

set log xy
set key right bottom

set xrange [0.004:0.055]

set xlabel 'Mesh Spacing'
set ylabel 'RMS of error'

plot './resolution.3D.gmsh.out' \
u 1:2 t '{/Symbol n} = 2' w l,\
'' u 1:3 t '{/Symbol n} = 3' w l,\
'' u 1:4 t '{/Symbol n} = 4' w l,\
'' u 1:5 t '{/Symbol n} = 5' w l;

+ 0
- 33
plots/resolution/resolution_3D.py View File

@@ -1,33 +0,0 @@
"""
This script was run to generate the resolution study. it requres the generated
sqlite3 db of results from the gmsh folder.
"""


import sys
import sqlite3
import numpy as np

from interp.tools import rms

con = sqlite3.connect(sys.argv[1])
cur = con.cursor()

orders = np.array(cur.execute('select distinct ord from results').fetchall())[:,0]
eps = np.array(cur.execute('select distinct ep from results').fetchall())[:,0]
res = np.array(cur.execute('select distinct res from results').fetchall())[:,0]

params = (
(2,64),
(3,64),
(4,64),
(5,64),
)
for r in res:
for order, ep in params:
# print np.array([np.abs(i[FINAL] - i[EXACT]) for i in data if i[ORDER] == order and i[EP] == ep])),
cur_results = cur.execute('select abs(exact - final) from results where res = ? and ep = ? and ord = ?',
(int(r), int(ep),int(order)))
cur_ = np.array(cur_results.fetchall())
print rms(np.array(cur_)),
print

+ 0
- 11
plots/rms/2D/rms.1-4.2.plt View File

@@ -1,11 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'rms.2.1-4.2.eps'

set title '2-D Interpolation Error Data, Order 2'

plot 'rms.out' u 1:2 t 'Mesh 1' w lp, \
'' u 1:6 t 'Mesh 2' w lp, \
'' u 1:10 t 'Mesh 3' w lp, \
'' u 1:14 t 'Mesh 4' w lp;

+ 0
- 11
plots/rms/2D/rms.1-4.3.plt View File

@@ -1,11 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'rms.2.1-4.3.eps'

set title '2-D Interpolation Error Data, Order 3'

plot 'rms.out' u 1:3 t 'Mesh 1' w lp, \
'' u 1:7 t 'Mesh 2' w lp, \
'' u 1:11 t 'Mesh 3' w lp, \
'' u 1:15 t 'Mesh 4' w lp;

+ 0
- 11
plots/rms/2D/rms.1-4.4.plt View File

@@ -1,11 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'rms.2.1-4.4.eps'

set title '2-D Interpolation Error Data, Order 4'

plot 'rms.out' u 1:4 t 'Mesh 1' w lp, \
'' u 1:8 t 'Mesh 2' w lp, \
'' u 1:12 t 'Mesh 3' w lp, \
'' u 1:16 t 'Mesh 4' w lp;

+ 0
- 11
plots/rms/2D/rms.1-4.5.plt View File

@@ -1,11 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'rms.2.1-4.5.eps'

set title '2-D Interpolation Error Data, Order 5'

plot 'rms.out' u 1:5 t 'Mesh 1' w lp, \
'' u 1:9 t 'Mesh 2' w lp, \
'' u 1:13 t 'Mesh 3' w lp, \
'' u 1:17 t 'Mesh 4' w lp;

+ 0
- 12
plots/rms/2D/rms.1.2-5.plt View File

@@ -1,12 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'rms.2.1.2-5.eps'

# set title 'Mesh 1'
set yrange [1e-10:]

plot 'rms.out' u 1:2 t '{/Symbol n} = 2' w p, \
'' u 1:3 t '{/Symbol n} = 3' w p, \
'' u 1:4 t '{/Symbol n} = 4' w p, \
'' u 1:5 t '{/Symbol n} = 5' w p;

+ 0
- 11
plots/rms/2D/rms.2.2-5.plt View File

@@ -1,11 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'rms.2.2.2-5.eps'

# set title 'Mesh 2'

plot 'rms.out' u 1:6 t '{/Symbol n} = 2' w p, \
'' u 1:7 t '{/Symbol n} = 3' w p, \
'' u 1:8 t '{/Symbol n} = 4' w p, \
'' u 1:9 t '{/Symbol n} = 5' w p;

+ 0
- 11
plots/rms/2D/rms.3.2-5.plt View File

@@ -1,11 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'rms.2.3.2-5.eps'

# set title 'Mesh 3'

plot 'rms.out' u 1:10 t '{/Symbol n} = 2' w p, \
'' u 1:11 t '{/Symbol n} = 3' w p, \
'' u 1:12 t '{/Symbol n} = 4' w p, \
'' u 1:13 t '{/Symbol n} = 5' w p;

+ 0
- 11
plots/rms/2D/rms.4.2-5.plt View File

@@ -1,11 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'rms.2.4.2-5.eps'

# set title 'Mesh 4'

plot 'rms.out' u 1:14 t '{/Symbol n} = 2' w p, \
'' u 1:15 t '{/Symbol n} = 3' w p, \
'' u 1:16 t '{/Symbol n} = 4' w p, \
'' u 1:17 t '{/Symbol n} = 5' w p;

+ 0
- 27
plots/rms/2D/rms.onering.plt View File

@@ -1,27 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'rms.2.aggregate.eps'

set title "2-D Interpolation Error Data, Aggregate"

plot 'rms.out' \
u 1:2 t 'Order 2, Mesh 1' w lp, \
'' u 1:3 t 'Order 3, Mesh 1' w lp, \
'' u 1:4 t 'Order 4, Mesh 1' w lp, \
'' u 1:5 t 'Order 5, Mesh 1' w lp, \
\
'' u 1:6 t 'Order 2, Mesh 2' w lp, \
'' u 1:7 t 'Order 3, Mesh 2' w lp, \
'' u 1:8 t 'Order 4, Mesh 2' w lp, \
'' u 1:9 t 'Order 5, Mesh 2' w lp, \
\
'' u 1:10 t 'Order 2, Mesh 3' w lp, \
'' u 1:11 t 'Order 3, Mesh 3' w lp, \
'' u 1:12 t 'Order 4, Mesh 3' w lp, \
'' u 1:13 t 'Order 5, Mesh 3' w lp, \
\
'' u 1:14 t 'Order 2, Mesh 3' w lp, \
'' u 1:15 t 'Order 3, Mesh 3' w lp, \
'' u 1:16 t 'Order 4, Mesh 3' w lp, \
'' u 1:17 t 'Order 5, Mesh 3' w lp;

+ 0
- 14
plots/rms/2D/rms.out View File

@@ -1,14 +0,0 @@
# i found these resolutions [1 2 3 4], orders: [2 3 4 5], and extra points: [ 4 6 8 12 16 20 32 48 64 96 128 192 256]
4 9.653693e-04 1.153226e-05 9.389895e-06 8.797822e-07 2.224437e-02 2.712190e-03 2.284719e-04 4.496180e-05 5.758518e-02 1.330732e-03 3.052339e-04 6.705222e-05 8.837961e-03 1.386682e-03 1.181360e-04 3.211895e-05
6 2.915436e-04 8.236172e-06 2.725166e-07 2.271404e-08 2.474664e-02 1.735101e-03 1.648392e-04 1.023013e-05 1.026672e-02 8.868542e-04 1.280747e-04 8.152238e-06 2.691455e-02 5.416025e-02 2.688884e-04 5.519579e-05
8 2.478579e-04 6.410394e-06 2.122306e-07 1.786786e-08 6.698900e-04 1.211739e-06 4.507145e-08 1.259001e-08 1.447646e-02 5.915672e-03 3.033947e-04 2.782823e-05 3.045717e-02 8.220426e-04 3.751131e-03 7.505716e-05
12 2.817965e-04 7.564166e-06 2.263891e-07 1.959486e-08 8.511874e-05 7.473324e-07 2.123070e-08 1.549617e-09 3.191349e-03 1.048907e-06 1.109421e-07 3.179551e-09 6.015045e-02 3.469820e-03 2.390970e-04 2.496696e-04
16 3.495202e-04 8.973882e-06 2.969632e-07 2.269972e-08 9.662410e-05 9.263365e-07 2.283407e-08 6.352336e-10 3.723101e-05 6.305445e-08 3.227009e-09 2.589288e-12 1.805396e-01 8.124292e-03 7.249283e-04 2.512040e-05
20 3.647239e-04 8.520724e-06 2.756078e-07 2.129599e-08 1.188420e-04 1.173823e-06 2.423284e-08 6.362691e-10 1.994467e-05 4.945311e-08 3.163340e-09 3.350373e-14 1.563370e-05 1.482477e-08 1.304885e-10 2.846050e-09
32 4.675779e-04 8.985031e-06 2.815930e-07 2.133162e-08 1.765312e-04 1.751708e-06 2.972239e-08 6.385047e-10 3.040437e-05 6.741910e-08 4.495782e-09 3.178679e-14 6.740579e-06 1.318622e-08 3.963718e-12 3.162278e-16
48 6.323092e-04 9.805009e-06 2.779090e-07 2.152389e-08 2.441426e-04 2.596996e-06 3.853129e-08 7.174750e-10 5.556858e-05 9.779271e-08 4.496267e-09 7.154649e-14 1.415039e-05 2.082282e-08 6.025949e-12 3.166228e-14
64 7.972738e-04 1.148514e-05 2.973366e-07 2.016172e-08 3.084407e-04 3.388851e-06 4.476765e-08 3.259496e-09 9.047942e-05 1.593592e-07 5.526528e-09 4.683235e-13 2.618652e-05 3.044288e-08 9.470078e-12 3.167965e-14
96 1.089351e-03 1.418822e-05 3.199458e-07 2.209300e-08 4.307070e-04 4.939473e-06 6.009096e-08 3.279036e-09 1.553180e-04 2.554138e-07 7.771500e-09 5.930467e-13 5.904557e-05 5.144691e-08 1.771492e-11 7.094012e-14
128 1.346511e-03 1.734560e-05 3.387161e-07 2.069123e-08 5.625393e-04 6.439954e-06 7.303612e-08 5.584036e-09 2.286338e-04 4.060982e-07 6.417359e-09 9.077802e-13 1.063582e-04 8.215802e-08 3.217094e-10 3.210894e-13
192 1.731959e-03 2.360744e-05 4.093571e-07 2.259225e-08 8.429591e-04 9.366251e-06 1.053834e-07 6.462134e-09 3.663338e-04 8.219567e-07 7.202510e-09 1.674440e-12 2.310849e-04 1.772614e-07 7.274127e-10 1.032741e-12
256 2.002899e-03 3.041290e-05 4.779848e-07 2.420134e-08 1.109571e-03 1.215176e-05 1.358001e-07 6.494317e-09 5.018140e-04 1.334592e-06 9.142651e-09 4.203941e-12 3.767699e-04 3.204174e-07 3.267461e-09 3.159821e-12

+ 0
- 10
plots/rms/3D/rms.1-3.2.plt View File

@@ -1,10 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'rms.3.1-3.2.eps'

set title "3-D Interpolation Error Data, Order 2"

plot 'rms.out' u 1:2 t 'Mesh 1' w lp, \
'' u 1:6 t 'Mesh 2' w lp, \
'' u 1:10 t 'Mesh 3' w lp;

+ 0
- 10
plots/rms/3D/rms.1-3.3.plt View File

@@ -1,10 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'rms.3.1-3.3.eps'

set title "3-D Interpolation Error Data, Order 3"

plot 'rms.out' u 1:3 t 'Mesh 1' w lp, \
'' u 1:7 t 'Mesh 2' w lp, \
'' u 1:11 t 'Mesh 3' w lp;

+ 0
- 10
plots/rms/3D/rms.1-3.4.plt View File

@@ -1,10 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'rms.3.1-3.4.eps'

set title "3-D Interpolation Error Data, Order 4"

plot 'rms.out' u 1:4 t 'Mesh 1' w lp, \
'' u 1:8 t 'Mesh 2' w lp, \
'' u 1:12 t 'Mesh 3' w lp;

+ 0
- 10
plots/rms/3D/rms.1-3.5.plt View File

@@ -1,10 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'rms.3.1-3.5.eps'

set title "3-D Interpolation Error Data, Order 5"

plot 'rms.out' u 1:5 t 'Mesh 1' w lp, \
'' u 1:9 t 'Mesh 2' w lp, \
'' u 1:13 t 'Mesh 3' w lp;

+ 0
- 13
plots/rms/3D/rms.1.2-5.plt View File

@@ -1,13 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'rms.3.1.2-5.eps'

set title "3-D Interpolation Error Data, Mesh 1"

unset xlabel

plot 'rms.out' u 1:2 t 'Order 2' w lp, \
'' u 1:3 t 'Order 3' w lp, \
'' u 1:4 t 'Order 4' w lp, \
'' u 1:5 t 'Order 5' w lp;

+ 0
- 13
plots/rms/3D/rms.2.2-5.plt View File

@@ -1,13 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'rms.3.2.2-5.eps'

set title "3-D Interpolation Error Data, Mesh 2"

unset xlabel

plot 'rms.out' u 1:6 t 'Order 2' w lp, \
'' u 1:7 t 'Order 3' w lp, \
'' u 1:8 t 'Order 4' w lp, \
'' u 1:9 t 'Order 5' w lp;

+ 0
- 11
plots/rms/3D/rms.3.2-5.plt View File

@@ -1,11 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'rms.3.3.2-5.eps'

set title "3-D Interpolation Error Data, Mesh 3"

plot 'rms.out' u 1:10 t 'Order 2' w lp, \
'' u 1:11 t 'Order 3' w lp, \
'' u 1:12 t 'Order 4' w lp, \
'' u 1:13 t 'Order 5' w lp;

+ 0
- 22
plots/rms/3D/rms.onering.plt View File

@@ -1,22 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'rms.3.aggregate.eps'

set title "3-D Interpolation Error Data, Aggregate"

plot 'rms.out' \
u 1:2 t 'Order 2, Mesh 1' w lp, \
'' u 1:3 t 'Order 3, Mesh 1' w lp, \
'' u 1:4 t 'Order 4, Mesh 1' w lp, \
'' u 1:5 t 'Order 5, Mesh 1' w lp, \
\
'' u 1:6 t 'Order 2, Mesh 2' w lp, \
'' u 1:7 t 'Order 3, Mesh 2' w lp, \
'' u 1:8 t 'Order 4, Mesh 2' w lp, \
'' u 1:9 t 'Order 5, Mesh 2' w lp, \
\
'' u 1:10 t 'Order 2, Mesh 3' w lp, \
'' u 1:11 t 'Order 3, Mesh 3' w lp, \
'' u 1:12 t 'Order 4, Mesh 3' w lp, \
'' u 1:13 t 'Order 5, Mesh 3' w lp;

+ 0
- 14
plots/rms/3D/rms.out View File

@@ -1,14 +0,0 @@
# i found these resolutions [1 2 3], orders: [2 3 4 5], and extra points: [ 4 6 8 12 16 20 32 48 64 96 128 192 256]
4 ?5.803865e+13 ?6.468366e+10 ?9.446612e+08 4.142812e-02 1.288835e-02 1.677748e-03 4.217798e-02 3.252782e-03 6.342185e-02 5.850643e+00 4.395137e-03 5.464553e-01
6 ?1.506279e+13 ?6.329637e+08 1.458046e-05 2.701640e-02 1.234589e-03 1.083698e-03 2.591163e-02 5.251461e-03 1.506222e-04 7.098341e-01 1.219063e-03 1.341399e-04
8 ?3.503680e+13 3.586842e-06 6.304341e-07 2.210574e-02 7.848561e-04 1.601867e-04 3.757592e-02 6.845328e-04 3.306416e-04 1.632203e-02 8.592921e-04 9.923875e-05
12 2.567755e-04 2.290057e-06 2.099488e-07 1.423990e+00 1.147174e-02 6.093459e-04 5.862785e-03 9.020114e-04 1.016040e-04 6.652001e-03 5.977288e-04 3.539456e-05
16 2.044888e-04 2.089695e-06 1.780128e-07 1.231168e-01 1.934193e-05 1.869495e-06 9.785679e-03 3.104862e-04 7.676741e-05 5.124481e-03 3.262511e-04 1.833029e-04
20 2.004471e-04 1.917685e-06 1.755959e-07 1.495925e-04 2.151168e-06 1.490905e-07 1.994276e-02 3.910980e-04 1.498353e-04 1.582137e-02 1.854517e-04 5.026134e-05
32 2.147791e-04 1.778664e-06 1.762228e-07 6.238698e-05 1.143586e-07 1.237475e-08 1.062038e-02 4.924117e-05 5.850621e-06 1.480912e-01 3.178775e-04 8.517430e-05
48 2.361992e-04 1.952376e-06 1.980072e-07 7.389347e-05 1.428927e-07 1.439005e-08 1.262571e-05 9.681331e-09 7.871841e-11 2.349378e-02 1.348102e-03 2.129725e-04
64 2.701086e-04 2.054004e-06 2.186953e-07 8.535004e-05 1.697754e-07 1.524417e-08 1.417916e-05 1.014127e-08 3.241580e-10 4.997944e-06 8.086028e-10 2.189358e-11
96 3.074006e-04 2.155472e-06 2.179027e-07 1.045950e-04 2.182197e-07 1.855175e-08 1.994606e-05 1.396961e-08 1.087190e-10 3.864147e-06 6.646445e-10 9.905357e-12
128 3.450126e-04 2.256681e-06 2.324303e-07 1.162485e-04 2.607905e-07 2.045799e-08 2.731520e-05 1.420313e-08 4.607101e-10 5.198090e-06 8.295068e-10 3.162330e-10
192 4.247757e-04 2.498312e-06 2.537302e-07 1.468682e-04 3.346440e-07 2.317954e-08 4.134303e-05 1.869158e-08 4.681613e-10 8.013911e-06 1.097134e-09 3.171439e-10
256 5.054097e-04 2.641297e-06 2.694648e-07 1.701789e-04 4.017138e-07 2.577219e-08 5.140416e-05 2.207550e-08 6.498217e-10 1.124241e-05 1.217000e-09 3.262365e-10

+ 0
- 6
plots/rms/all.plt View File

@@ -1,6 +0,0 @@
set xlabel 'Number of Extra Points, S_k'
set ylabel 'RMS of total error'

set log xy

set key bottom left

+ 0
- 12
plots/scalability/efficiency.plt View File

@@ -1,12 +0,0 @@
load '../main.plt'

set output 'efficiency.eps'

set nokey

set ylabel "Efficiency (E_p)"
set xlabel "Number of Participating Minions"

set yrange [0:1]

plot 'scale.out' u 1:3 w p

+ 0
- 22
plots/scalability/hist.py View File

@@ -1,22 +0,0 @@
import sys
import shelve
import pickle

import numpy as np

if len(sys.argv) != 2:
print "usage: %s <interp.shelve>" % sys.argv[0]
sys.exit(1)


s = shelve.open(sys.argv[1])
d = dict(s)
s.close()

for line in (i[1] for i in sorted(d.iteritems(), key = lambda x: x[1]['stats']['participants'])):
run = line['stats']
print "#", run['participants']
x = np.array(run['tasks'].values())
for i in x:
print i
print

+ 0
- 35
plots/scalability/scalability_histograms.py View File

@@ -1,35 +0,0 @@
import sys
import shelve
import pickle

import numpy as np
import matplotlib
matplotlib.use('Agg')
import matplotlib.pyplot as plt

if len(sys.argv) != 2:
print "usage: %s <interp.shelve>" % sys.argv[0]
sys.exit(1)


s = shelve.open(sys.argv[1])
d = dict(s)
s.close()

c = 0
for line in (i[1] for i in sorted(d.iteritems(), key = lambda x: x[1]['stats']['participants'])):
c += 1
run = line['stats']
x = np.array(run['tasks'].values())

# the histogram of the data
plt.figure(c)

plt.title("Histogram of Interpolation Performance")
plt.xlabel("Number of Interpolations Performed")
plt.ylabel("Number of minions")

n, bins, patches = plt.hist(x) # , 50, normed=1, facecolor='green', alpha=0.5)
plt.grid(True)
output_file_name = 'asdf.%0.3d.%d.png' % (run['participants'], str(run['count']).count("0"))
plt.savefig(output_file_name)

+ 0
- 26
plots/scalability/scalability_plot_gen.py View File

@@ -1,26 +0,0 @@
import sys
import shelve


s = shelve.open(sys.argv[1])
d = dict(s)
s.close()

# {'count': 10000, 'extra': 256, 'receive': 1148.0290439128876, 'submit': 4.298105001449585, 'participants': 1, 'order': 5}
# 'tasks': defaultdict(<type 'int'>, {'m5-2-1.local-31671': 10000}),

for line in (i[1] for i in sorted(d.iteritems(), key = lambda x: x[1]['stats']['participants'])):
run = line['stats']

thruput = run['count'] / run['receive']

p = run['participants']

# speedup (http://en.wikipedia.org/wiki/Speedup):
T_1 = 0.114802904391
S_p = T_1 / (run['receive']/run['count'])

E_p = S_p / p

print p, S_p, E_p


+ 0
- 19
plots/scalability/scale.out View File

@@ -1,19 +0,0 @@
1 0.999999999997 0.999999999997
2 1.94216782645 0.971083913223
4 3.8842715987 0.971067899674
5 4.74931333527 0.949862667054
8 8.1352967993 ?1.01691209991
11 10.9213356195 0.992848692681
16 16.1043018639 ?1.00651886649
22 21.6914424408 0.985974656402
32 31.6479069695 0.988997092798
45 44.4369078566 0.987486841258
64 63.198130042 0.987470781906
90 88.930924868 0.988121387422
128 126.013638476 0.984481550592
181 163.291098498 0.902160765184
196 172.310663615 0.879136038854
224 135.769780068 0.606115089587
256 150.571313747 0.588169194325
384 139.374902817 0.362955476086
512 156.074759811 0.304833515255

+ 0
- 12
plots/scalability/server.out View File

@@ -1,12 +0,0 @@
8 2
16 7
45 30
64 40
90 65
128 100
181 160
196 175
224 185
256 195
384 200
512 200

+ 0
- 12
plots/scalability/server.plt View File

@@ -1,12 +0,0 @@
load '../main.plt'

set output 'server.eps'

set nokey

set xlabel "Number of Participating Minions"
set ylabel "Approximate CPU utilization (% CPU)"

set yrange [0:210]

plot 'server.out' u 1:2 w p t 'Speedup'

+ 0
- 10
plots/scalability/speedup.plt View File

@@ -1,10 +0,0 @@
load '../main.plt'

set output 'speedup.eps'

set nokey

set ylabel "Speedup (S_p)"
set xlabel "Number of Participating Minions"

plot 'scale.out' u 1:2 w p t 'Speedup'

+ 0
- 27
plots/timing/2D/onering.plt View File

@@ -1,27 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'timing.2.aggregate.eps'

unset xlabel
# set title '2-D Timing'

plot 'time.out' u 1:2 t 'Order 2, res 1' w lp, \
'' u 1:3 t 'Order 3, res 1' w lp, \
'' u 1:4 t 'Order 4, res 1' w lp, \
'' u 1:5 t 'Order 5, res 1' w lp, \
\
'' u 1:6 t 'Order 2, res 2' w lp, \
'' u 1:7 t 'Order 3, res 2' w lp, \
'' u 1:8 t 'Order 4, res 2' w lp, \
'' u 1:9 t 'Order 5, res 2' w lp, \
\
'' u 1:10 t 'Order 2, res 3' w lp, \
'' u 1:11 t 'Order 3, res 3' w lp, \
'' u 1:12 t 'Order 4, res 3' w lp, \
'' u 1:13 t 'Order 5, res 3' w lp, \
\
'' u 1:14 t 'Order 2, res 4' w lp, \
'' u 1:15 t 'Order 3, res 4' w lp, \
'' u 1:16 t 'Order 4, res 4' w lp, \
'' u 1:17 t 'Order 5, res 4' w lp;

+ 0
- 20
plots/timing/2D/time.1-4.2.plt View File

@@ -1,20 +0,0 @@
load '../../main.plt'
load '../all.plt'


set title '2-D Order 2 Timing'

set output 'time.2.1-4.2.eps'
plot 'time.out' u 1:2 t 'Mesh 1' w lp, \
'' u 1:6 t 'Mesh 2' w lp, \
'' u 1:10 t 'Mesh 3' w lp, \
'' u 1:14 t 'Mesh 4' w lp;

set output 'time.1.eps'

unset xlabel
plot 'time.out' u 1:2 t 'Mesh 1' w lp, \
'' u 1:6 t 'Mesh 2' w lp, \
'' u 1:10 t 'Mesh 3' w lp, \
'' u 1:14 t 'Mesh 4' w lp;


+ 0
- 11
plots/timing/2D/time.1-4.3.plt View File

@@ -1,11 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'time.2.1-4.3.eps'

set title '2-D Order 3 Timing'

plot 'time.out' u 1:3 t 'Mesh 1' w lp, \
'' u 1:7 t 'Mesh 2' w lp, \
'' u 1:11 t 'Mesh 3' w lp, \
'' u 1:15 t 'Mesh 4' w lp;

+ 0
- 11
plots/timing/2D/time.1-4.4.plt View File

@@ -1,11 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'time.2.1-4.4.eps'

set title '2-D Order 4 Timing'

plot 'time.out' u 1:4 t 'Mesh 1' w lp, \
'' u 1:8 t 'Mesh 2' w lp, \
'' u 1:12 t 'Mesh 3' w lp, \
'' u 1:16 t 'Mesh 4' w lp;

+ 0
- 21
plots/timing/2D/time.1-4.5.plt View File

@@ -1,21 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'time.2.1-4.5.eps'

set title '2-D Order 5 Timing'

plot 'time.out' u 1:5 t 'Mesh 1' w lp, \
'' u 1:9 t 'Mesh 2' w lp, \
'' u 1:13 t 'Mesh 3' w lp, \
'' u 1:17 t 'Mesh 4' w lp;

set output 'time.2.eps'

unset xlabel
unset ylabel

plot 'time.out' u 1:5 t 'Mesh 1' w lp, \
'' u 1:9 t 'Mesh 2' w lp, \
'' u 1:13 t 'Mesh 3' w lp, \
'' u 1:17 t 'Mesh 4' w lp;

+ 0
- 11
plots/timing/2D/time.1.2-5.plt View File

@@ -1,11 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'time.2.1.2-5.eps'

set title '2-D Mesh 1 Timing'

plot 'time.out' u 1:2 t '{/Symbol n} = 2' w p, \
'' u 1:3 t '{/Symbol n} = 3' w p, \
'' u 1:4 t '{/Symbol n} = 4' w p, \
'' u 1:5 t '{/Symbol n} = 5' w p;

+ 0
- 11
plots/timing/2D/time.2.2-5.plt View File

@@ -1,11 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'time.2.2.2-5.eps'

set title '2-D Mesh 2 Timing'

plot 'time.out' u 1:6 t '{/Symbol n} = 2' w p, \
'' u 1:7 t '{/Symbol n} = 3' w p, \
'' u 1:8 t '{/Symbol n} = 4' w p, \
'' u 1:9 t '{/Symbol n} = 5' w p;

+ 0
- 11
plots/timing/2D/time.3.2-5.plt View File

@@ -1,11 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'time.2.3.2-5.eps'

set title '2-D Mesh 3 Timing'

plot 'time.out' u 1:10 t '{/Symbol n} = 2' w p, \
'' u 1:11 t '{/Symbol n} = 3' w p, \
'' u 1:12 t '{/Symbol n} = 4' w p, \
'' u 1:13 t '{/Symbol n} = 5' w p;

+ 0
- 11
plots/timing/2D/time.4.2-5.plt View File

@@ -1,11 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'time.2.4.2-5.eps'

set title '2-D Mesh 4 Timing'

plot 'time.out' u 1:14 t '{/Symbol n} = 2' w p, \
'' u 1:15 t '{/Symbol n} = 3' w p, \
'' u 1:16 t '{/Symbol n} = 4' w p, \
'' u 1:17 t '{/Symbol n} = 5' w p;

+ 0
- 14
plots/timing/2D/time.out View File

@@ -1,14 +0,0 @@
# i found these resolutions [1 2 3 4], orders: [2 3 4 5], and extra points: [ 4 6 8 12 16 20 32 48 64 96 128 192 256]
4 ?8.418652e-03 3.655416e-03 3.898457e-03 ?5.967275e-03 3.366499e-03 3.630556e-03 3.967421e-03 4.110508e-03 3.460547e-03 3.771054e-03 4.046290e-03 4.253433e-03 3.589308e-03 3.961004e-03 4.245772e-03 4.426731e-03
6 3.635591e-03 4.052618e-03 4.306668e-03 4.506541e-03 3.741066e-03 4.042297e-03 4.311079e-03 4.517118e-03 3.821324e-03 4.200951e-03 4.493928e-03 4.678068e-03 4.033009e-03 4.415525e-03 4.722003e-03 4.867687e-03
8 4.015927e-03 4.457202e-03 4.726348e-03 4.927077e-03 4.070683e-03 4.436711e-03 4.744807e-03 4.950891e-03 4.262951e-03 4.669412e-03 4.971558e-03 5.169330e-03 4.518866e-03 4.932510e-03 5.243011e-03 5.405859e-03
12 4.818652e-03 5.337744e-03 5.631116e-03 5.856908e-03 4.927565e-03 5.341905e-03 5.677934e-03 5.919848e-03 5.183425e-03 5.659151e-03 5.993829e-03 6.202919e-03 5.562590e-03 6.028122e-03 6.395465e-03 6.563699e-03
16 5.574414e-03 6.136970e-03 6.433900e-03 6.698951e-03 5.730611e-03 6.182127e-03 6.540695e-03 6.819537e-03 6.062893e-03 6.576594e-03 6.952685e-03 7.177316e-03 6.566743e-03 7.060981e-03 7.447590e-03 7.661908e-03
20 6.333457e-03 6.961534e-03 7.269407e-03 7.533588e-03 6.530682e-03 7.038809e-03 7.463932e-03 7.779570e-03 6.996469e-03 7.520437e-03 7.941494e-03 8.192018e-03 7.621880e-03 8.114247e-03 8.561446e-03 8.789118e-03
32 8.602427e-03 9.298019e-03 9.799798e-03 1.011445e-02 8.989583e-03 9.586569e-03 1.017834e-02 1.052187e-02 9.762936e-03 1.031987e-02 1.088784e-02 1.117662e-02 1.073170e-02 1.123353e-02 1.183430e-02 1.213934e-02
48 1.172114e-02 1.244392e-02 1.321292e-02 1.348411e-02 1.240673e-02 1.298877e-02 1.381157e-02 1.420465e-02 1.349452e-02 1.403827e-02 1.496650e-02 1.527638e-02 1.488403e-02 1.529200e-02 1.631513e-02 1.668931e-02
64 1.477927e-02 1.543926e-02 1.651096e-02 1.677096e-02 1.571161e-02 1.624126e-02 1.727634e-02 1.778081e-02 1.705369e-02 1.757175e-02 1.883169e-02 1.918908e-02 1.890653e-02 1.933516e-02 2.068021e-02 2.116732e-02
96 2.091696e-02 2.158498e-02 2.320393e-02 2.359935e-02 2.224743e-02 2.284590e-02 2.447949e-02 2.497847e-02 2.423504e-02 2.487589e-02 2.685860e-02 2.721632e-02 2.707030e-02 2.770393e-02 2.962969e-02 3.030726e-02
128 2.692815e-02 2.769574e-02 2.986872e-02 3.043787e-02 2.875882e-02 2.940888e-02 3.169157e-02 3.237488e-02 3.136918e-02 3.213348e-02 3.480177e-02 3.542475e-02 3.517311e-02 3.672961e-02 3.858470e-02 3.941328e-02
192 3.898889e-02 4.004129e-02 4.325339e-02 4.400221e-02 4.187063e-02 4.269522e-02 4.625397e-02 4.710145e-02 4.580470e-02 4.784677e-02 5.074559e-02 5.154375e-02 5.141480e-02 5.434784e-02 5.628456e-02 5.773756e-02
256 5.087446e-02 5.210588e-02 5.631130e-02 5.736148e-02 5.475548e-02 5.651388e-02 6.040707e-02 6.139642e-02 5.998892e-02 6.304337e-02 6.609333e-02 6.764806e-02 6.752912e-02 7.154435e-02 7.400029e-02 7.559165e-02

+ 0
- 22
plots/timing/3D/onering.plt View File

@@ -1,22 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'timing.3.aggregate.eps'


# set title '3-D Timing'

plot 'time.out' u 1:2 t 'Order 2, Mesh 1' w lp, \
'' u 1:3 t 'Order 3, Mesh 1' w lp, \
'' u 1:4 t 'Order 4, Mesh 1' w lp, \
'' u 1:5 t 'Order 5, Mesh 1' w lp, \
\
'' u 1:6 t 'Order 2, Mesh 2' w lp, \
'' u 1:7 t 'Order 3, Mesh 2' w lp, \
'' u 1:8 t 'Order 4, Mesh 2' w lp, \
'' u 1:9 t 'Order 5, Mesh 2' w lp, \
\
'' u 1:10 t 'Order 2, Mesh 3' w lp, \
'' u 1:11 t 'Order 3, Mesh 3' w lp, \
'' u 1:12 t 'Order 4, Mesh 3' w lp, \
'' u 1:13 t 'Order 5, Mesh 3' w lp;

+ 0
- 16
plots/timing/3D/time.1-3.2.plt View File

@@ -1,16 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'time.3.1-3.2.eps'

set title '3-D Order 2 Timing'

plot 'time.out' u 1:2 t 'Mesh 1' w lp, \
'' u 1:6 t 'Mesh 2' w lp, \
'' u 1:10 t 'Mesh 3' w lp;

set output 'time.3.eps'

plot 'time.out' u 1:2 t 'Mesh 1' w lp, \
'' u 1:6 t 'Mesh 2' w lp, \
'' u 1:10 t 'Mesh 3' w lp;

+ 0
- 10
plots/timing/3D/time.1-3.3.plt View File

@@ -1,10 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'time.3.1-3.3.eps'

set title '3-D Order 3 Timing'

plot 'time.out' u 1:3 t 'Mesh 1' w lp, \
'' u 1:7 t 'Mesh 2' w lp, \
'' u 1:11 t 'Mesh 3' w lp;

+ 0
- 10
plots/timing/3D/time.1-3.4.plt View File

@@ -1,10 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'time.3.1-3.4.eps'

set title '3-D Order 4 Timing'

plot 'time.out' u 1:4 t 'Mesh 1' w lp, \
'' u 1:8 t 'Mesh 2' w lp, \
'' u 1:12 t 'Mesh 3' w lp;

+ 0
- 18
plots/timing/3D/time.1-3.5.plt View File

@@ -1,18 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'time.3.1-3.5.eps'

set title '3-D Order 5 Timing'

plot 'time.out' u 1:5 t 'Mesh 1' w lp, \
'' u 1:9 t 'Mesh 2' w lp, \
'' u 1:13 t 'Mesh 3' w lp;

set output 'time.4.eps'

unset ylabel

plot 'time.out' u 1:5 t 'Mesh 1' w lp, \
'' u 1:9 t 'Mesh 2' w lp, \
'' u 1:13 t 'Mesh 3' w lp;

+ 0
- 11
plots/timing/3D/time.1.2-5.plt View File

@@ -1,11 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'time.3.1.2-5.eps'

set title '3-D Mesh 1 Timing'

plot 'time.out' u 1:2 t 'Order 2' w lp, \
'' u 1:3 t 'Order 3' w lp, \
'' u 1:4 t 'Order 4' w lp, \
'' u 1:5 t 'Order 5' w lp;

+ 0
- 11
plots/timing/3D/time.2.2-5.plt View File

@@ -1,11 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'time.3.2.2-5.eps'

set title '3-D Mesh 2 Timing'

plot 'time.out' u 1:6 t 'Order 2' w lp, \
'' u 1:7 t 'Order 3' w lp, \
'' u 1:8 t 'Order 4' w lp, \
'' u 1:9 t 'Order 5' w lp;

+ 0
- 11
plots/timing/3D/time.3.2-5.plt View File

@@ -1,11 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'time.3.3.2-5.eps'

set title '3-D Mesh 3 Timing'

plot 'time.out' u 1:10 t 'Order 2' w lp, \
'' u 1:11 t 'Order 3' w lp, \
'' u 1:12 t 'Order 4' w lp, \
'' u 1:13 t 'Order 5' w lp;

+ 0
- 14
plots/timing/3D/time.out View File

@@ -1,14 +0,0 @@
# i found these resolutions [1 2 3], orders: [2 3 4 5], and extra points: [ 4 6 8 12 16 20 32 48 64 96 128 192 256]
4 5.336315e-03 7.098224e-03 ?1.451139e-02 5.387542e-03 6.369798e-03 6.753827e-03 5.837306e-03 6.806631e-03 6.940095e-03 6.452059e-03 7.355210e-03 7.433917e-03
6 5.963069e-03 6.748175e-03 6.965468e-03 5.834118e-03 6.859724e-03 7.203237e-03 6.302207e-03 7.442983e-03 7.636667e-03 7.018682e-03 8.144408e-03 8.274832e-03
8 6.171001e-03 7.272441e-03 7.468094e-03 6.347245e-03 7.439834e-03 7.794676e-03 6.873887e-03 8.064043e-03 8.327575e-03 7.785675e-03 8.986467e-03 9.142167e-03
12 7.136869e-03 8.348455e-03 8.588379e-03 7.416486e-03 8.666553e-03 9.031010e-03 8.147807e-03 9.415864e-03 9.784622e-03 9.313637e-03 1.070905e-02 1.096778e-02
16 8.018513e-03 9.349211e-03 9.611884e-03 8.423628e-03 9.813188e-03 1.018199e-02 9.357845e-03 1.078487e-02 1.115367e-02 1.079949e-02 1.234503e-02 1.272981e-02
20 8.900174e-03 1.036481e-02 1.065289e-02 9.431523e-03 1.096711e-02 1.137259e-02 1.057305e-02 1.215228e-02 1.261098e-02 1.251150e-02 1.411872e-02 1.452954e-02
32 1.153077e-02 1.333737e-02 1.371999e-02 1.244510e-02 1.442168e-02 1.483082e-02 1.422781e-02 1.622985e-02 1.679338e-02 1.750419e-02 1.934132e-02 1.982753e-02
48 1.505026e-02 1.736507e-02 1.781956e-02 1.643499e-02 1.893480e-02 1.951993e-02 1.931718e-02 2.181412e-02 2.241628e-02 2.421645e-02 2.644494e-02 2.694771e-02
64 1.850048e-02 2.147960e-02 2.178370e-02 2.053152e-02 2.347319e-02 2.406672e-02 2.440554e-02 2.722164e-02 2.792115e-02 3.059527e-02 3.344584e-02 3.406788e-02
96 2.536127e-02 2.941491e-02 2.976743e-02 2.863635e-02 3.257535e-02 3.322455e-02 3.423788e-02 3.818909e-02 3.905295e-02 4.357597e-02 4.763484e-02 4.835693e-02
128 3.234847e-02 3.704390e-02 3.776606e-02 3.677154e-02 4.159797e-02 4.234546e-02 4.410054e-02 4.934341e-02 5.028840e-02 5.630809e-02 6.185975e-02 6.276128e-02
192 4.620687e-02 5.245157e-02 5.339440e-02 5.261651e-02 5.932565e-02 6.048002e-02 6.357378e-02 7.131608e-02 7.240405e-02 8.184434e-02 9.045512e-02 9.157925e-02
256 6.004110e-02 6.794921e-02 6.920238e-02 6.847277e-02 7.733731e-02 7.873610e-02 8.321401e-02 9.367546e-02 9.479373e-02 1.076390e-01 1.198657e-01 1.205021e-01

+ 0
- 6
plots/timing/all.plt View File

@@ -1,6 +0,0 @@
set xlabel 'Number of Extra Points, S_k'
set ylabel 'Average Interpolation Time (s)'

set yrange [0.0:0.15]

set key top left

+ 0
- 42
plots/timing/timing.aggregate.plt View File

@@ -1,42 +0,0 @@
load '../main.plt'
load 'all.plt'

set output "timing.aggregate.eps"

set nokey

plot '3D/time.out'\
u 1:2 t 'O 2, res 1' w lp, \
'' u 1:3 t 'dim 3, O 3, res 1' w lp, \
'' u 1:4 t 'dim 3, O 4, res 1' w lp, \
'' u 1:5 t 'dim 3, O 5, res 1' w lp, \
\
'' u 1:6 t 'dim 3, O 2, res 2' w lp, \
'' u 1:7 t 'dim 3, O 3, res 2' w lp, \
'' u 1:8 t 'dim 3, O 4, res 2' w lp, \
'' u 1:9 t 'dim 3, O 5, res 2' w lp, \
\
'' u 1:10 t 'dim 3, O 2, res 3' w lp, \
'' u 1:11 t 'dim 3, O 3, res 3' w lp, \
'' u 1:12 t 'dim 3, O 4, res 3' w lp, \
'' u 1:13 t 'dim 3, O 5, res 3' w lp, \
\
'2D/time.out' u 1:2 t 'O 2, res 1' w lp, \
'' u 1:3 t 'dim 2, O 3, res 1' w lp, \
'' u 1:4 t 'dim 2, O 4, res 1' w lp, \
'' u 1:5 t 'dim 2, O 5, res 1' w lp, \
\
'' u 1:6 t 'dim 2, O 2, res 2' w lp, \
'' u 1:7 t 'dim 2, O 3, res 2' w lp, \
'' u 1:8 t 'dim 2, O 4, res 2' w lp, \
'' u 1:9 t 'dim 2, O 5, res 2' w lp, \
\
'' u 1:10 t 'dim 2, O 2, res 3' w lp, \
'' u 1:11 t 'dim 2, O 3, res 3' w lp, \
'' u 1:12 t 'dim 2, O 4, res 3' w lp, \
'' u 1:13 t 'dim 2, O 5, res 3' w lp, \
\
'' u 1:14 t 'dim 2, O 2, res 4' w lp, \
'' u 1:15 t 'dim 2, O 3, res 4' w lp, \
'' u 1:16 t 'dim 2, O 4, res 4' w lp, \
'' u 1:17 t 'dim 2, O 5, res 4' w lp;

+ 0
- 27
plots/truthiness/2D/onering.plt View File

@@ -1,27 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'truthy.2.aggregate.eps'

set title "2-D Improvement, Aggregate"

plot 'truthy.out' \
u 1:2 t 'Order 2, Mesh 1' w lp, \
'' u 1:3 t 'Order 3, Mesh 1' w lp, \
'' u 1:4 t 'Order 4, Mesh 1' w lp, \
'' u 1:5 t 'Order 5, Mesh 1' w lp, \
\
'' u 1:6 t 'Order 2, Mesh 2' w lp, \
'' u 1:7 t 'Order 3, Mesh 2' w lp, \
'' u 1:8 t 'Order 4, Mesh 2' w lp, \
'' u 1:9 t 'Order 5, Mesh 2' w lp, \
\
'' u 1:10 t 'Order 2, Mesh 3' w lp, \
'' u 1:11 t 'Order 3, Mesh 3' w lp, \
'' u 1:12 t 'Order 4, Mesh 3' w lp, \
'' u 1:13 t 'Order 5, Mesh 3' w lp, \
\
'' u 1:14 t 'Order 2, Mesh 4' w lp, \
'' u 1:15 t 'Order 3, Mesh 4' w lp, \
'' u 1:16 t 'Order 4, Mesh 4' w lp, \
'' u 1:17 t 'Order 5, Mesh 4' w lp;

+ 0
- 12
plots/truthiness/2D/truthy.1-4.2.plt View File

@@ -1,12 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'truthy.2.1-4.2.eps'

set title "2-D Improvement, Order 2"

plot 'truthy.out' \
u 1:2 t 'Mesh 1' w lp, \
'' u 1:6 t 'Mesh 2' w lp, \
'' u 1:10 t 'Mesh 3' w lp, \
'' u 1:14 t 'Mesh 4' w lp;

+ 0
- 12
plots/truthiness/2D/truthy.1-4.3.plt View File

@@ -1,12 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'truthy.2.1-4.3.eps'

set title "2-D Improvement, Order 3"

plot 'truthy.out' \
u 1:3 t 'Mesh 1' w lp, \
'' u 1:7 t 'Mesh 2' w lp, \
'' u 1:11 t 'Mesh 3' w lp, \
'' u 1:15 t 'Mesh 4' w lp;

+ 0
- 12
plots/truthiness/2D/truthy.1-4.4.plt View File

@@ -1,12 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'truthy.2.1-4.4.eps'

set title "2-D Improvement, Order 4"

plot 'truthy.out' \
u 1:4 t 'Mesh 1' w lp, \
'' u 1:8 t 'Mesh 2' w lp, \
'' u 1:12 t 'Mesh 3' w lp, \
'' u 1:16 t 'Mesh 4' w lp;

+ 0
- 12
plots/truthiness/2D/truthy.1-4.5.plt View File

@@ -1,12 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'truthy.2.1-4.5.eps'

set title "2-D Improvement, Order 5"

plot 'truthy.out' \
u 1:5 t 'Mesh 1' w lp, \
'' u 1:9 t 'Mesh 2' w lp, \
'' u 1:13 t 'Mesh 3' w lp, \
'' u 1:17 t 'Mesh 4' w lp;

+ 0
- 13
plots/truthiness/2D/truthy.1.2-5.plt View File

@@ -1,13 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'truthy.2.1.2-5.eps'
set key bottom lef

# set title 'Mesh 1'

plot 'truthy.out' \
u 1:2 t '{/Symbol n} = 2' w p, \
'' u 1:3 t '{/Symbol n} = 3' w p, \
'' u 1:4 t '{/Symbol n} = 4' w p, \
'' u 1:5 t '{/Symbol n} = 5' w p;

+ 0
- 12
plots/truthiness/2D/truthy.2.2-5.plt View File

@@ -1,12 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'truthy.2.2.2-5.eps'

# set title 'Mesh 2'

plot 'truthy.out' \
u 1:6 t '{/Symbol n} = 2' w p, \
'' u 1:7 t '{/Symbol n} = 3' w p, \
'' u 1:8 t '{/Symbol n} = 4' w p, \
'' u 1:9 t '{/Symbol n} = 5' w p;

+ 0
- 12
plots/truthiness/2D/truthy.3.2-5.plt View File

@@ -1,12 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'truthy.2.3.2-5.eps'

# set title 'Mesh 3'

plot 'truthy.out' \
u 1:10 t '{/Symbol n} = 2' w p, \
'' u 1:11 t '{/Symbol n} = 3' w p, \
'' u 1:12 t '{/Symbol n} = 4' w p, \
'' u 1:13 t '{/Symbol n} = 5' w p;

+ 0
- 12
plots/truthiness/2D/truthy.4.2-5.plt View File

@@ -1,12 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'truthy.2.4.2-5.eps'

# set title 'Mesh 4'

plot 'truthy.out' \
u 1:14 t '{/Symbol n} = 2' w p, \
'' u 1:15 t '{/Symbol n} = 3' w p, \
'' u 1:16 t '{/Symbol n} = 4' w p, \
'' u 1:17 t '{/Symbol n} = 5' w p;

+ 0
- 14
plots/truthiness/2D/truthy.out View File

@@ -1,14 +0,0 @@
# i found these resolutions [1 2 3 4], orders: [2 3 4 5], and extra points: [ 4 6 8 12 16 20 32 48 64 96 128 192 256]
4 9.020000e-01 9.700000e-01 9.900000e-01 9.870000e-01 3.840000e-01 4.050000e-01 3.840000e-01 4.080000e-01 4.360000e-01 4.520000e-01 4.230000e-01 4.280000e-01 4.240000e-01 4.570000e-01 4.520000e-01 5.050000e-01
6 9.290000e-01 9.790000e-01 9.980000e-01 9.990000e-01 5.620000e-01 6.110000e-01 5.840000e-01 5.880000e-01 4.620000e-01 4.660000e-01 4.530000e-01 4.780000e-01 4.790000e-01 4.500000e-01 4.770000e-01 4.880000e-01
8 9.410000e-01 9.830000e-01 9.990000e-01 1.000000e+00 9.710000e-01 9.970000e-01 1.000000e+00 1.000000e+00 4.340000e-01 4.750000e-01 4.520000e-01 4.790000e-01 4.630000e-01 4.880000e-01 4.960000e-01 5.000000e-01
12 9.290000e-01 9.810000e-01 9.970000e-01 9.960000e-01 9.820000e-01 9.980000e-01 1.000000e+00 1.000000e+00 9.600000e-01 9.990000e-01 1.000000e+00 1.000000e+00 4.780000e-01 4.790000e-01 4.880000e-01 5.250000e-01
16 9.210000e-01 9.740000e-01 9.980000e-01 9.960000e-01 9.820000e-01 9.980000e-01 1.000000e+00 1.000000e+00 9.940000e-01 1.000000e+00 1.000000e+00 1.000000e+00 4.960000e-01 4.850000e-01 4.740000e-01 5.150000e-01
20 9.270000e-01 9.770000e-01 9.990000e-01 9.980000e-01 9.770000e-01 9.980000e-01 1.000000e+00 1.000000e+00 9.970000e-01 1.000000e+00 1.000000e+00 1.000000e+00 9.960000e-01 1.000000e+00 1.000000e+00 1.000000e+00
32 9.090000e-01 9.760000e-01 9.980000e-01 9.960000e-01 9.660000e-01 9.980000e-01 9.990000e-01 1.000000e+00 9.910000e-01 1.000000e+00 1.000000e+00 1.000000e+00 1.000000e+00 1.000000e+00 1.000000e+00 1.000000e+00
48 8.970000e-01 9.770000e-01 9.960000e-01 9.960000e-01 9.530000e-01 9.990000e-01 9.990000e-01 1.000000e+00 9.870000e-01 1.000000e+00 1.000000e+00 1.000000e+00 9.940000e-01 1.000000e+00 1.000000e+00 1.000000e+00
64 8.760000e-01 9.710000e-01 9.960000e-01 9.960000e-01 9.360000e-01 9.940000e-01 9.990000e-01 1.000000e+00 9.840000e-01 1.000000e+00 1.000000e+00 1.000000e+00 9.890000e-01 1.000000e+00 1.000000e+00 1.000000e+00
96 8.440000e-01 9.710000e-01 9.950000e-01 9.960000e-01 9.090000e-01 9.890000e-01 9.990000e-01 1.000000e+00 9.680000e-01 1.000000e+00 1.000000e+00 1.000000e+00 9.840000e-01 1.000000e+00 1.000000e+00 1.000000e+00
128 8.130000e-01 9.680000e-01 9.950000e-01 9.960000e-01 8.820000e-01 9.860000e-01 9.990000e-01 1.000000e+00 9.580000e-01 9.990000e-01 1.000000e+00 1.000000e+00 9.710000e-01 1.000000e+00 1.000000e+00 1.000000e+00
192 7.380000e-01 9.590000e-01 9.950000e-01 9.960000e-01 8.410000e-01 9.810000e-01 9.990000e-01 9.990000e-01 9.220000e-01 9.990000e-01 1.000000e+00 1.000000e+00 9.490000e-01 9.980000e-01 1.000000e+00 1.000000e+00
256 6.770000e-01 9.520000e-01 9.930000e-01 9.960000e-01 8.040000e-01 9.750000e-01 9.970000e-01 9.990000e-01 8.680000e-01 9.980000e-01 1.000000e+00 1.000000e+00 9.240000e-01 9.980000e-01 1.000000e+00 1.000000e+00

+ 0
- 22
plots/truthiness/3D/onering.plt View File

@@ -1,22 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'truthy.3.aggregate.eps'

set title "3-D Improvement, Aggregate"

plot 'truthy.out' \
u 1:2 t 'Order 2, Mesh 1' w lp, \
'' u 1:3 t 'Order 3, Mesh 1' w lp, \
'' u 1:4 t 'Order 4, Mesh 1' w lp, \
'' u 1:5 t 'Order 5, Mesh 1' w lp, \
\
'' u 1:6 t 'Order 2, Mesh 2' w lp, \
'' u 1:7 t 'Order 3, Mesh 2' w lp, \
'' u 1:8 t 'Order 4, Mesh 2' w lp, \
'' u 1:9 t 'Order 5, Mesh 2' w lp, \
\
'' u 1:10 t 'Order 2, Mesh 3' w lp, \
'' u 1:11 t 'Order 3, Mesh 3' w lp, \
'' u 1:12 t 'Order 4, Mesh 3' w lp, \
'' u 1:13 t 'Order 5, Mesh 3' w lp;

+ 0
- 11
plots/truthiness/3D/truthy.1-3.2.plt View File

@@ -1,11 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'truthy.3.1-3.2.eps'

set title "3-D Improvement, Order 2"

plot 'truthy.out' \
u 1:2 t 'Mesh 1' w lp, \
'' u 1:6 t 'Mesh 2' w lp, \
'' u 1:10 t 'Mesh 3' w lp;

+ 0
- 11
plots/truthiness/3D/truthy.1-3.3.plt View File

@@ -1,11 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'truthy.3.1-3.3.eps'

set title "3-D Improvement, Order 3"

plot 'truthy.out' \
u 1:3 t 'Mesh 1' w lp, \
'' u 1:7 t 'Mesh 2' w lp, \
'' u 1:11 t 'Mesh 3' w lp;

+ 0
- 11
plots/truthiness/3D/truthy.1-3.4.plt View File

@@ -1,11 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'truthy.3.1-3.4.eps'

set title "3-D Improvement, Order 4"

plot 'truthy.out' \
u 1:4 t 'Mesh 1' w lp, \
'' u 1:8 t 'Mesh 2' w lp, \
'' u 1:12 t 'Mesh 3' w lp;

+ 0
- 11
plots/truthiness/3D/truthy.1-3.5.plt View File

@@ -1,11 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'truthy.3.1-3.5.eps'

set title "3-D Improvement, Order 5"

plot 'truthy.out' \
u 1:5 t 'Mesh 1' w lp, \
'' u 1:9 t 'Mesh 2' w lp, \
'' u 1:13 t 'Mesh 3' w lp;

+ 0
- 14
plots/truthiness/3D/truthy.1.2-5.plt View File

@@ -1,14 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'truthy.3.1.2-5.eps'

set title "3-D Improvement, Mesh 1"

unset xlabel

plot 'truthy.out' \
u 1:2 t 'Order 2' w lp, \
'' u 1:3 t 'Order 3' w lp, \
'' u 1:4 t 'Order 4' w lp, \
'' u 1:5 t 'Order 5' w lp;

+ 0
- 14
plots/truthiness/3D/truthy.2.2-5.plt View File

@@ -1,14 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'truthy.3.2.2-5.eps'

set title "3-D Improvement, Mesh 2"

unset xlabel

plot 'truthy.out' \
u 1:6 t 'Order 2' w lp, \
'' u 1:7 t 'Order 3' w lp, \
'' u 1:8 t 'Order 4' w lp, \
'' u 1:9 t 'Order 5' w lp;

+ 0
- 12
plots/truthiness/3D/truthy.3.2-5.plt View File

@@ -1,12 +0,0 @@
load '../../main.plt'
load '../all.plt'

set output 'truthy.3.3.2-5.eps'

set title "3-D Improvement, Mesh 3"

plot 'truthy.out' \
u 1:10 t 'Order 2' w lp, \
'' u 1:11 t 'Order 3' w lp, \
'' u 1:12 t 'Order 4' w lp, \
'' u 1:13 t 'Order 5' w lp;

+ 0
- 14
plots/truthiness/3D/truthy.out View File

@@ -1,14 +0,0 @@
# i found these resolutions [1 2 3], orders: [2 3 4 5], and extra points: [ 4 6 8 12 16 20 32 48 64 96 128 192 256]
4 3.430000e-01 3.240000e-01 3.230000e-01 3.280000e-01 3.200000e-01 3.150000e-01 3.660000e-01 3.730000e-01 3.690000e-01 4.130000e-01 4.310000e-01 4.130000e-01
6 6.740000e-01 8.790000e-01 9.270000e-01 3.390000e-01 3.400000e-01 3.230000e-01 3.740000e-01 3.900000e-01 4.040000e-01 4.430000e-01 4.460000e-01 4.120000e-01
8 8.330000e-01 9.600000e-01 9.880000e-01 3.250000e-01 3.520000e-01 3.420000e-01 3.880000e-01 4.100000e-01 3.990000e-01 4.400000e-01 4.440000e-01 4.380000e-01
12 8.930000e-01 9.810000e-01 9.910000e-01 3.710000e-01 3.920000e-01 3.920000e-01 4.230000e-01 4.080000e-01 3.960000e-01 4.230000e-01 4.830000e-01 4.390000e-01
16 9.110000e-01 9.800000e-01 9.920000e-01 8.420000e-01 9.700000e-01 9.920000e-01 4.390000e-01 4.250000e-01 4.130000e-01 4.410000e-01 4.730000e-01 4.600000e-01
20 9.140000e-01 9.820000e-01 9.920000e-01 9.390000e-01 9.970000e-01 9.980000e-01 4.110000e-01 4.220000e-01 4.090000e-01 4.600000e-01 4.780000e-01 4.530000e-01
32 9.120000e-01 9.830000e-01 9.930000e-01 9.660000e-01 9.990000e-01 9.990000e-01 9.120000e-01 9.810000e-01 9.950000e-01 4.430000e-01 4.460000e-01 4.550000e-01
48 9.070000e-01 9.820000e-01 9.910000e-01 9.520000e-01 9.970000e-01 9.980000e-01 9.850000e-01 9.990000e-01 1.000000e+00 3.920000e-01 4.290000e-01 3.890000e-01
64 9.000000e-01 9.810000e-01 9.930000e-01 9.500000e-01 9.970000e-01 9.980000e-01 9.860000e-01 9.990000e-01 1.000000e+00 9.910000e-01 1.000000e+00 1.000000e+00
96 8.830000e-01 9.830000e-01 9.930000e-01 9.430000e-01 9.970000e-01 9.990000e-01 9.840000e-01 9.990000e-01 9.990000e-01 9.890000e-01 9.990000e-01 1.000000e+00
128 8.920000e-01 9.790000e-01 9.910000e-01 9.310000e-01 9.970000e-01 9.990000e-01 9.850000e-01 9.990000e-01 9.990000e-01 9.870000e-01 9.990000e-01 1.000000e+00
192 8.640000e-01 9.800000e-01 9.940000e-01 9.230000e-01 9.940000e-01 9.990000e-01 9.730000e-01 9.990000e-01 9.990000e-01 9.870000e-01 9.990000e-01 1.000000e+00
256 8.590000e-01 9.820000e-01 9.910000e-01 9.080000e-01 9.930000e-01 9.990000e-01 9.740000e-01 9.990000e-01 9.990000e-01 9.880000e-01 9.990000e-01 1.000000e+00

+ 0
- 8
plots/truthiness/all.plt View File

@@ -1,8 +0,0 @@
set xlabel "Number of Extra Points, S_k"
set ylabel "Percent of Improved Interpolations"

set log x

# set yrange [0:1]

set key bottom right

Loading…
Cancel
Save