Merge pull request #125 from sandialabs/eparish1/85-opinf-opinf

added rom-rom support to norma

Author: ikalash

examples/ahead/overlap/cuboid/dynamic-opinf-fom/make_op_inf_models_rom_rom_1.py

@@ -0,0 +1,26 @@
+import normaopinf
+import normaopinf.opinf
+import nnopinf
+import nnopinf.training
+import os
+import numpy as np
+
+if __name__ == '__main__':
+  settings = {}
+  settings['fom-yaml-file'] = "cuboid-1.yaml"
+  settings['training-data-directories'] = [os.getcwd()]
+  settings['model-type'] = 'linear'
+  settings['stop-training-time'] = 'end'
+  settings['training-skip-steps'] = 1
+  settings['forcing'] =  False
+  settings['truncation-type'] = 'energy'
+  settings['boundary-truncation-type'] =  'energy'
+  settings['regularization-parameter'] =  5.e-3
+  settings['trial-space-splitting-type'] = 'split'
+  settings['acceleration-computation-type'] = 'finite-difference'
+  snapshots_dict = normaopinf.opinf.get_processed_snapshots(settings)
+  settings['truncation-value'] = 1. - 1.e-5
+  settings['boundary-truncation-value'] = 1. - 1.e-5
+  settings['model-name'] = 'opinf-operator-1'
+  normaopinf.opinf.make_opinf_model_from_snapshots_dict(snapshots_dict,settings)
+

examples/ahead/overlap/cuboid/dynamic-opinf-fom/make_op_inf_models_rom_rom_2.py

@@ -0,0 +1,26 @@
+import normaopinf
+import normaopinf.opinf
+import nnopinf
+import nnopinf.training
+import os
+import numpy as np
+
+if __name__ == '__main__':
+  settings = {}
+  settings['fom-yaml-file'] = "cuboid-2.yaml"
+  settings['training-data-directories'] = [os.getcwd()]
+  settings['model-type'] = 'linear'
+  settings['stop-training-time'] = 'end'
+  settings['training-skip-steps'] = 1
+  settings['forcing'] =  False
+  settings['truncation-type'] = 'energy'
+  settings['boundary-truncation-type'] =  'energy'
+  settings['regularization-parameter'] =  5.e-3
+  settings['trial-space-splitting-type'] = 'split'
+  settings['acceleration-computation-type'] = 'finite-difference'
+  snapshots_dict = normaopinf.opinf.get_processed_snapshots(settings)
+  settings['truncation-value'] = 1. - 1.e-5
+  settings['boundary-truncation-value'] = 1. - 1.e-5
+  settings['model-name'] = 'opinf-operator-2'
+  normaopinf.opinf.make_opinf_model_from_snapshots_dict(snapshots_dict,settings)
+

examples/ahead/overlap/cuboid/dynamic-opinf-rom-rom/cuboid-1.g

@@ -0,0 +1,41 @@
+type: single
+input mesh file: cuboid-1.g
+output mesh file: cuboid-1.e
+model:
+  type: linear opinf rom
+  model-file: opinf-operator-1.npz
+  material:
+    blocks:
+      fine: hyperelastic
+    hyperelastic:
+      model: neohookean
+      elastic modulus: 1.0e+09
+      Poisson's ratio: 0.25
+      density: 1000.0
+time integrator:
+  type: Newmark
+  β: 0.25
+  γ: 0.5
+boundary conditions:
+  Dirichlet:
+    - node set: nsx-
+      component: x
+      function: "0.0"
+    - node set: nsy-
+      component: y
+      function: "0.0"
+    - node set: nsz-
+      component: z
+      function: "0.0"
+  Schwarz overlap:
+    - side set: ssz+
+      source: cuboid-2
+      source block: coarse
+      source side set: ssz-
+solver:
+  type: Hessian minimizer
+  step: full Newton
+  minimum iterations: 1
+  maximum iterations: 16
+  relative tolerance: 1.0e-10
+  absolute tolerance: 1.0e-06

examples/ahead/overlap/cuboid/dynamic-opinf-rom-rom/cuboid-1.yaml

@@ -0,0 +1,35 @@
+type: single
+input mesh file: cuboid-2.g
+output mesh file: cuboid-2.e
+model:
+  type: linear opinf rom
+  model-file: opinf-operator-2.npz
+  material:
+    blocks:
+      coarse: hyperelastic
+    hyperelastic:
+      model: neohookean
+      elastic modulus: 1.0e+09
+      Poisson's ratio: 0.25
+      density: 1000.0
+time integrator:
+  type: Newmark
+  β: 0.25
+  γ: 0.5
+boundary conditions:
+  Dirichlet:
+    - node set: nsz+
+      component: z
+      function: "1.0 * t"
+  Schwarz overlap:
+    - side set: ssz-
+      source: cuboid-1
+      source block: fine
+      source side set: ssz+
+solver:
+  type: Hessian minimizer
+  step: full Newton
+  minimum iterations: 1
+  maximum iterations: 16
+  relative tolerance: 1.0e-10
+  absolute tolerance: 1.0e-06

examples/ahead/overlap/cuboid/dynamic-opinf-rom-rom/cuboid-2.g

@@ -0,0 +1,11 @@
+type: multi
+domains: ["cuboid-1.yaml", "cuboid-2.yaml"]
+Exodus output interval: 1
+CSV output interval: 1
+initial time: 0.0
+final time: 0.1
+time step: 0.01
+minimum iterations: 1
+maximum iterations: 16
+relative tolerance: 1.0e-12
+absolute tolerance: 1.0e-08

examples/ahead/overlap/cuboid/dynamic-opinf-rom-rom/cuboid-2.yaml

@@ -425,7 +425,7 @@ function apply_bc_detail(model::SolidMechanics, bc::SolidMechanicsNonOverlapSchw
 end
 
 function apply_bc_detail(model::OpInfModel, bc::SolidMechanicsCouplingSchwarzBoundaryCondition)
-    if bc.coupled_subsim.model isa SolidMechanics
+    if bc.coupled_subsim.model isa SolidMechanics || bc.coupled_subsim.model isa OpInfModel
         ## Apply BC to the FOM vector
         apply_bc_detail(model.fom_model, bc)