PyANSYS Improvements and Demos (#106)

<!--
Please use a PR title that conforms to *conventional commits*:
"<commit_type>: Describe your change"; for example: "fix: prevent race
condition". Some other commit types are: fix, feat, ci, doc, refactor...
For a full list of commit types visit
https://www.conventionalcommits.org/en/v1.0.0/
-->

#### Relevant issue or PR
<!-- If the changes resolve an issue or follow some other PR, link to
them here. Only link something if it is directly relevant. -->

#### Description of changes
<!-- Add a high-level description of changes, focusing on the *what* and
*why*. -->

#### Testing done
<!-- Describe how the changes were tested; e.g., "CI passes", "Tested
manually in stagingrepo#123", screenshots of a terminal session that
verify the changes, or any other evidence of testing the changes. -->

Co-authored-by: Andrin Rehmann <[email protected]>

Author: NolanBlack

examples/ansys/optim_bars_pymapdl.ipynb

@@ -53,10 +53,10 @@ class HexMesh(BaseModel):
         description="Array of hexahedral faces defined by indices into the points array."
     )
     n_points: Int32 = Field(
-        default=0, description="Number of valid points in the points array."
+        default = 0, description="Number of valid points in the points array."
     )
     n_faces: Int32 = Field(
-        default=0, description="Number of valid faces in the faces array."
+        default = 0, description="Number of valid faces in the faces array."
     )
 
 
@@ -158,12 +158,15 @@ def __init__(self, inputs: InputSchema, mapdl: Mapdl) -> None:
 
         # Extract input parameters
         self.rho = inputs.rho
-        self.hex_mesh = inputs.hex_mesh
         self.dirichlet_mask = inputs.dirichlet_mask
         self.dirichlet_values = inputs.dirichlet_values
         self.van_neumann_mask = inputs.van_neumann_mask
         self.van_neumann_values = inputs.van_neumann_values
         self.hex_mesh = inputs.hex_mesh
+        if self.hex_mesh.n_faces == 0:
+            self.hex_mesh.n_faces = self.hex_mesh.faces.shape[0]
+        if self.hex_mesh.n_points == 0:
+            self.hex_mesh.n_points = self.hex_mesh.points.shape[0]
         self.E0 = inputs.E0
         self.rho_min = inputs.rho_min
         self.p = inputs.p
@@ -669,4 +672,9 @@ def vector_jacobian_product(
 
 def abstract_eval(abstract_inputs):
     """Calculate output shape of apply from the shape of its inputs."""
-    return {"compliance": ShapeDType(shape=(), dtype="float32")}
+    n_faces = abstract_inputs.hex_mesh.faces.shape[0]
+    return {
+        "compliance": ShapeDType(shape=(), dtype="float32"),
+        "strain_energy": ShapeDType(shape=(n_faces,), dtype="float32"),
+        "sensitivity": ShapeDType(shape=(n_faces,), dtype="float32"),
+    }