basic automatic benchmark

.gitlab-ci.yml

@@ -71,9 +71,7 @@ mseatest :
   stage : test
   script:
     - pip3 install --user dist/*manylinux1_x86_64*
-    - cd demo
-    - curl -O ftp://braque.mema.ucl.ac.be/msea/data-test-1.tar.gz
-    - tar xf data-test-1.tar.gz
+    - cd tests
     - python3 testshp.py 
 
 mfdeploy-test :

demo/testshp.py → tests/testshp.py

+import os
+import urllib.request
+import tarfile
+if not os.path.isdir("data") :
+    urllib.request.urlretrieve("ftp://braque.mema.ucl.ac.be/msea/data-test-1.tar.gz", "data-test-1.tar.gz")
+    f = tarfile.open("data-test-1.tar.gz","r:*")
+    f.extractall()
+
 import msea
 from msea.geometry import CurveType
+import msea.geometry
 import numpy as np
 from osgeo import osr 
 
@@ -14,12 +23,18 @@ domain_srs.ImportFromEPSG(32631)
 #domain_srs.ImportFromProj4("+proj=utm +ellps=WGS84 +zone=31")
 
 domain = msea.geometry.Domain(domain_srs)
-domain.add_boundary_curves_shp("test/data_no_duplicate.shp","physical",CurveType.POLYLINE)
-#domain.add_interior_curves_shp("test/interior.shp",None,CurveType.STRICTPOLYLINE)
-#domain.add_interior_points_shp("test/interior.shp")
-bath_field = msea.field.Raster("medit.tiff")
+domain.add_boundary_curves_shp("data/data_no_duplicate.shp","physical",CurveType.POLYLINE)
+#domain.add_interior_curves_shp("data/interior.shp",None,CurveType.STRICTPOLYLINE)
+#domain.add_interior_points_shp("data/interior.shp")
+bath_field = msea.field.Raster("data/medit.tiff")
 dist_coast_field = msea.field.Distance(domain,100,["coast","island"])
 dist_porquerolles_field = msea.field.Distance(domain,20,["porquerolles"])
 #coarse = msea.geometry.coarsen_boundaries(domain,(8e5,4.68e6),domain_srs,mesh_size,20)
-msea.gmsh.mesh(domain,"test.msh",mesh_size,intermediate_file_name="log")
+msea.gmsh.mesh(domain,"data.msh",mesh_size,intermediate_file_name="log")
 msea.gmsh.convert_to_gis("test.msh",domain_srs,"test.gpkg")
+
+msea.gmsh.gmsh.model.add("test")
+msea.gmsh.gmsh.open("test.msh")
+tag,nodes = msea.gmsh.gmsh.model.mesh.getElementsByType(2)
+ntri = len(tag)
+assert(ntri>40000 and ntri<41000)