update benchmark

test/perf/benchmark-R-ncdf4.R

@@ -11,6 +11,11 @@ print(paste("ncdf4 version: ",packageVersion("ncdf4")))
 fname = "filename_fv.nc"
 
 process <- function(fname) {
+  # drop file caches; requires root
+  fileConn<-file("/proc/sys/vm/drop_caches",open = "wt")
+  writeLines("3", fileConn)
+  close(fileConn)
+
   nc = nc_open(fname)
 
   # how do you get the dimension from the file?

test/perf/benchmark-julia-NCDatasets.jl

@@ -14,6 +14,9 @@ function compute(v)
 end
 
 function process(fname)
+    # drop file caches; requires root
+    write("/proc/sys/vm/drop_caches","3")
+
     ds = NCDataset(fname,"r") do ds
         v = ds["v1"];
         tot = compute(v)
@@ -24,7 +27,7 @@ end
 fname = "filename_fv.nc"
 tot = process(fname)
 
-@show tot
+println("result ",tot)
 
 bm = run(@benchmarkable process(fname) samples=100 seconds=10000)
 

test/perf/benchmark-python-netCDF4.py

@@ -10,20 +10,25 @@ def compute(v):
     return tot/v.shape[0]
 
 def process(fname):
+    with open("/proc/sys/vm/drop_caches","w") as f:
+        f.write("3")
+
     with netCDF4.Dataset(fname) as ds:
         v = ds["v1"]
         tot = compute(v)
         return tot
 
-def process_example():
+
+if __name__ == "__main__":
     fname = "filename_fv.nc";
-    process(fname)
+    tot = process(fname)
 
+    print("result ",tot)
 
-setup = "from __main__ import process_example"
-print("python-netCDF4 version ",netCDF4.__version__)
+    setup = "from __main__ import process"
+    print("python-netCDF4 version ",netCDF4.__version__)
 
-benchtime = timeit.repeat("process_example()", setup=setup,number = 1, repeat = 100)
-with open("python-netCDF4.txt","w") as f:
-    for bt in benchtime:
-        print(bt,file=f)
+    benchtime = timeit.repeat(lambda: process(fname), setup=setup,number = 1, repeat = 100)
+    with open("python-netCDF4.txt","w") as f:
+        for bt in benchtime:
+            print(bt,file=f)