Docs: updating docs (cassiope bot)

Author: cassiope

docs/doc/Examples/OCC/driver4.py

@@ -4,51 +4,61 @@
 import Generator
 import Converter
 
-# Create parameter
+# Create a parameter
 epaisseur = D.Scalar(12., name='epaisseur')
 epaisseur.range = [10,15]
 
 # discrete profile
 naca = Geom.naca(12, N=51)
 bbox = Generator.bbox(naca)
 
-# Create grid
+# Create parameter grid
 grid1 = D.Grid(bbox[0:3], bbox[3:], N=(3,3,1))
-grid1.P[1][2][0].y.range = [0, 5]
+grid1.P[1][2][0].y.range = [0, 5, 0.1]
 D.Eq(epaisseur.s, grid1.P[1][2][0].x.s)
 
-# Create profile
+# Create parametric profile
 spline1 = D.Spline3( grid1, mesh=naca, name='spline1' )
 
-# Create sketch
+# Create parametric sketch
 sketch1 = D.Sketch([spline1], name='sketch1')
 
-# test
+# solve for free parameters
 D.DRIVER.solve2()
-
 #grid1.P[1][2][0].y.print()
 
+# instantiate a CAD from free parameters
+# then mesh and get sensibilities
 D.DRIVER.instantiate({'P.1.2.0.y': 0.8})
-
 sketch1.writeCAD('out.step')
-
 mesh = sketch1.mesh(0.01, 0.01, 0.01)
 D.DRIVER._diff(sketch1, mesh)
 Converter.convertArrays2File(mesh, 'out.plt')
 
-D.DRIVER.setDOE({'P.1.2.0.y': 0.1})
+# Build DOE
+#D.DRIVER.setDOE({'P.1.2.0.y': 0.1})
+D.DRIVER.setDOE()
 D.DRIVER.createDOE('doe.hdf')
 D.DRIVER.walkDOE(sketch1, 0.01, 0.01, 0.01)
-# reread from doe file
+
+# reread one snapshot from DOE file
 m = D.DRIVER.readSnaphot(0)
 Converter.convertArrays2File(m, 'reread.plt')
-full = D.DRIVER.readAllSnapshots()
-print(full.shape)
-#print(full)
-D.DRIVER.fullSvd(full)
 
-import sys; sys.exit()
+# read snapshots as matrix
+F = D.DRIVER.readAllSnapshots()
+# Create a ROM limited to K modes
+Phi,S,Vt = D.DRIVER.createROM(F, K=20)
+D.DRIVER.writeROM('rom.hdf')
+# add to file the coordinates of snapshots on POD vectors
+D.DRIVER.addAllCoefs()
+
+# reread and build a snapshot from ROM
+coords = D.DRIVER.readCoefs(0)
+m = D.DRIVER.evalROM(coords)
+Converter.convertArrays2File(m, 'reread2.plt')
 
+# instantiate CADs, mesh and display
 import CPlot, time
 for i in range(20):
     D.DRIVER.instantiate({'P.1.2.0.y': 0.3+i/50.})

docs/doc/Examples/OCC/driver5.py

@@ -1,31 +1,52 @@
-# driver: parametric circle (derivatives)
+# driver: parametric circle (with derivatives)
 import OCC.Driver as D
 import Converter
 
-# Create parameter
+# Create a parameter
 radius = D.Scalar(1., name='radius')
-radius.range = [0.1,10]
+radius.range = [0.1, 10, 0.3]
 
-# Create circle
+# Create parametric circle
 circle1 = D.Circle( (0,0,0), radius, name='circle1' )
 
-# Create sketch
+# Create parametric sketch
 sketch1 = D.Sketch([circle1], name='sketch1')
 
-# test
+# solve for free parameters
 D.DRIVER.solve2()
 
+# instantiate a CAD
 D.DRIVER.instantiate({'radius': 1.5})
-
 sketch1.writeCAD('out.step')
-
 mesh = sketch1.mesh(0.01, 0.01, 0.01)
 D.DRIVER._diff(sketch1, mesh)
 Converter.convertArrays2File(mesh, 'out.plt')
 
-#import CPlot, time
-#for i in range(50):
-#    D.DRIVER.instantiate({'radius': 10-i/10.})
-#    mesh = sketch1.mesh(0.01, 0.01, 0.01)
-#    CPlot.display(mesh)
-#    time.sleep(0.5)
+# Build DOE
+D.DRIVER.setDOE()
+D.DRIVER.createDOE('doe.hdf')
+D.DRIVER.walkDOE(sketch1, 0.01, 0.01, 0.01)
+
+# reread one snaphsot from DOE file
+m = D.DRIVER.readSnaphot(0)
+Converter.convertArrays2File(m, 'reread.plt')
+
+# read snapshots as matrix
+F = D.DRIVER.readAllSnapshots()
+D.DRIVER.createROM(F, K=-1)
+D.DRIVER.writeROM('rom.hdf')
+# add to file the coordinates of snapshots on POD vectors
+D.DRIVER.addAllCoefs()
+
+# reread and build a snapshot from ROM
+coords = D.DRIVER.readCoefs(0)
+m = D.DRIVER.evalROM(coords)
+Converter.convertArrays2File(m, 'reread2.plt')
+
+# instantiate CADs, mesh and display
+import CPlot, time
+for i in range(50):
+    D.DRIVER.instantiate({'radius': 10-i/10.})
+    mesh = sketch1.mesh(0.01, 0.01, 0.01)
+    CPlot.display(mesh)
+    time.sleep(0.5)

docs/doc/Examples/OCC/driver8.py

@@ -0,0 +1,40 @@
+# driver: parametric profile in parallel
+import OCC.Driver as D
+import Geom
+import Generator
+import Converter
+
+# Create a parameter
+epaisseur = D.Scalar(12., name='epaisseur')
+epaisseur.range = [10,15]
+
+# discrete profile
+naca = Geom.naca(12, N=51)
+bbox = Generator.bbox(naca)
+
+# Create parameter grid
+grid1 = D.Grid(bbox[0:3], bbox[3:], N=(3,3,1))
+grid1.P[1][2][0].y.range = [0, 5, 0.1]
+D.Eq(epaisseur.s, grid1.P[1][2][0].x.s)
+
+# Create parametric profile
+spline1 = D.Spline3( grid1, mesh=naca, name='spline1' )
+
+# Create parametric sketch
+sketch1 = D.Sketch([spline1], name='sketch1')
+
+# solve for free parameters
+D.DRIVER.solve2()
+
+# Build DOE
+D.DRIVER.setDOE()
+D.DRIVER.createDOE('doe.hdf')
+D.DRIVER.walkDOE(sketch1, 0.01, 0.01, 0.01)
+
+# read snapshots as matrix
+#F = D.DRIVER.readAllSnapshots()
+# Create a ROM limited to K modes
+#Phi,S,Vt = D.DRIVER.createROM(F, K=20)
+#D.DRIVER.writeROM('rom.hdf')
+# add to file the coordinates of snapshots on POD vectors
+#D.DRIVER.addAllCoefs()

docs/doc/Examples/Post/integMomentNormPT_t1.py

@@ -4,6 +4,7 @@
 import Generator.PyTree as G
 import Post.PyTree as P
 import KCore.test as test
+test.TOLERANCE = 1.e-10
 
 ni = 11; nj = 11
 def f1(x,y): return 2*x + y

docs/doc/Examples/Post/integMomentPT_t1.py

@@ -3,6 +3,7 @@
 import Generator.PyTree as G
 import Post.PyTree as P
 import KCore.test as test
+test.TOLERANCE = 1.e-10
 
 ni = 11; nj = 11
 def f1(x,y): return 2*x + y

docs/doc/Examples/Post/integNormPT_m1.py

@@ -6,6 +6,7 @@
 import Post.Mpi as Pmpi
 import Converter.Mpi as Cmpi
 import KCore.test as test
+test.TOLERANCE = 1.e-10
 
 def f1(x,y): return 2*x + y
 def f2(x,y): return 3*x*y + 4