Some more clean up and fix boundary numbering (#303)

* CFD boundaries are reordered

* Formatting

* Formatting 2

* Formatting 3

* Create _plot_sizing.py

Author: bburke38

examples/fun3d_examples/_ssw-testing/1_panel_thickness.py

@@ -145,7 +145,9 @@
 # <----------------------------------------------------
 
 # make a funtofem scenario
-cruise = Scenario.steady("cruise_inviscid", steps=300, coupling_frequency=30, uncoupled_steps=0)
+cruise = Scenario.steady(
+    "cruise_inviscid", steps=300, coupling_frequency=30, uncoupled_steps=0
+)
 cruise.adjoint_steps = (
     100  # outer coupling iterations, total 5000 flow adjoints, 100 grid adjoints
 )
@@ -286,7 +288,13 @@
 manager.register_to_problem(opt_problem)
 
 # run an SNOPT optimization
-snoptimizer = SNOPT(options={"Verify level": 0, "Function precision": 1e-4, "Major Optimality tol" : 1e-4})
+snoptimizer = SNOPT(
+    options={
+        "Verify level": 0,
+        "Function precision": 1e-4,
+        "Major Optimality tol": 1e-4,
+    }
+)
 
 sol = snoptimizer(
     opt_problem,

examples/fun3d_examples/_ssw-testing/_oneway_sizing.py

@@ -13,7 +13,7 @@
 Local machine optimization for the panel thicknesses using nribs-1 OML panels and nribs-1 LE panels
 """
 
-#from pyoptsparse import SLSQP, Optimization
+# from pyoptsparse import SLSQP, Optimization
 from pyoptsparse import SNOPT, Optimization
 
 # script inputs

examples/fun3d_examples/_ssw-testing/design/_plot_sizing.py

@@ -0,0 +1,27 @@
+from funtofem import *
+
+# 1 : plot 1_sizing_opt_local.py results
+# ---------------------------------------------------------------------------
+plotter = PlotManager.from_hist_file(
+    "ssw-sizing1_design.txt",
+    accepted_names=["cruise_inviscid-mass", "cruise_inviscid-ksfailure"],
+    plot_names=["mass", "ks-stress"],
+    ignore_other_names=True,
+)
+
+# MAKE EACH PLOT FOR A DIFFERENT MODE
+functions = []
+Function.plot("ks-stress").optimize(scale=1.0, plot_name="ks-stress").register_to(
+    plotter
+)
+plotter.add_constraint(value=1.0, name="ks-constr")
+Function.plot("mass").optimize(scale=0.01, plot_name="mass").register_to(plotter)
+
+plotter(
+    plot_name="sizing-hist.png",
+    legend_frac=0.9,
+    yaxis_name="function values",
+    color_offset=0,
+    niceplots_style="doumont-light",
+    yscale_log=True,
+)

examples/fun3d_examples/_ssw-testing/geometry/_mesh_fun3d_egads.py

@@ -8,7 +8,7 @@
 # Set whether to build an inviscid or viscous mesh
 # ------------------------------------------------
 case = "inviscid"
-#case = "turbulent"
+# case = "turbulent"
 
 if case == "inviscid":
     project_name = "ssw-inviscid"

examples/fun3d_examples/ssw_meshdef_optimization/1_panel_thickness.py

@@ -80,7 +80,7 @@
 #         theta_init=0.6, theta_max=0.95, history_file=aitken_file, debug=True
 #     )
 # )
-wing = Body.aeroelastic("wing", boundary=3)
+wing = Body.aeroelastic("wing", boundary=2)
 
 # setup the material and shell properties
 aluminum = caps2tacs.Isotropic.aluminum().register_to(tacs_model)

examples/fun3d_examples/ssw_meshdef_optimization/2_aero_aoa.py

@@ -79,7 +79,7 @@
 # BODIES AND STRUCT DVs
 # <----------------------------------------------------
 
-wing = Body.aeroelastic("wing", boundary=3)
+wing = Body.aeroelastic("wing", boundary=2)
 
 # setup the material and shell properties
 aluminum = caps2tacs.Isotropic.aluminum().register_to(tacs_model)

examples/fun3d_examples/ssw_meshdef_optimization/3_geom_twist.py

@@ -128,7 +128,7 @@
 # BODIES AND STRUCT DVs
 # <----------------------------------------------------
 
-wing = Body.aeroelastic("wing", boundary=3)
+wing = Body.aeroelastic("wing", boundary=2)
 
 # setup the material and shell properties
 aluminum = caps2tacs.Isotropic.aluminum().register_to(tacs_model)

examples/fun3d_examples/ssw_meshdef_optimization/4_oml_shape.py

@@ -128,7 +128,7 @@
 # BODIES AND STRUCT DVs
 # <----------------------------------------------------
 
-wing = Body.aeroelastic("wing", boundary=3)
+wing = Body.aeroelastic("wing", boundary=2)
 
 # setup the material and shell properties
 aluminum = caps2tacs.Isotropic.aluminum().register_to(tacs_model)

examples/fun3d_examples/ssw_meshdef_optimization/_2_cfd_only.py

@@ -34,7 +34,7 @@
 # BODIES AND STRUCT DVs
 # <----------------------------------------------------
 
-wing = Body.aeroelastic("wing", boundary=3)
+wing = Body.aeroelastic("wing", boundary=2)
 
 # register the wing body to the model
 wing.register_to(f2f_model)

examples/fun3d_examples/ssw_meshdef_optimization/_oneway_sizing.py

@@ -75,7 +75,7 @@
 # BODIES AND STRUCT DVs
 # -------------------------------------------------
 
-wing = Body.aeroelastic("wing", boundary=3)
+wing = Body.aeroelastic("wing", boundary=2)
 
 # setup the material and shell properties
 aluminum = caps2tacs.Isotropic.aluminum().register_to(tacs_model)

examples/fun3d_examples/ssw_meshdef_optimization/_run_flow.py

@@ -20,7 +20,7 @@
 
 # Construct the FUNtoFEM model
 f2f_model = FUNtoFEMmodel("ssw_flow")
-wing = Body.aeroelastic("wing", boundary=3)
+wing = Body.aeroelastic("wing", boundary=2)
 wing.register_to(f2f_model)
 
 # Make a FUNtoFEM scenario

examples/fun3d_examples/ssw_meshdef_optimization/cfd/cruise/Flow/ssw-turb.mapbc

@@ -1,4 +1,4 @@
 3
-1 5000 Farfield
-2 6662 SymmetryY
-3 4000 wing
+1 5000 farfield
+2 4000 wing
+3 6662 symmetry

examples/fun3d_examples/ssw_meshdef_optimization/cfd/inv-cruise/Flow/ssw-turb.mapbc

@@ -1,4 +1,5 @@
 3
-1 5000 Farfield
-2 6662 SymmetryY
-3 3000 wing
+1 5000 farfield
+2 3000 wing
+3 6662 symmetry
+

funtofem/interface/test_solver.py

@@ -1083,7 +1083,11 @@ def finite_difference(
         f_functions = [func.value.real for func in model.get_functions(all=True)]
 
         finite_diff_TD = [
-            (f_functions[ifunc] - i_functions[ifunc]) / 2 / epsilon if central_diff else (f_functions[ifunc] - m_functions[ifunc])/epsilon
+            (
+                (f_functions[ifunc] - i_functions[ifunc]) / 2 / epsilon
+                if central_diff
+                else (f_functions[ifunc] - m_functions[ifunc]) / epsilon
+            )
             for ifunc in range(nfunctions)
         ]