EasyFEA 1.2.2__tar.gz → 1.2.4__tar.gz

{easyfea-1.2.2 → easyfea-1.2.4}/EasyFEA/Simulations.py

@@ -5,7 +5,7 @@
 """This module contains all available simulation classes."""
 
 from .simulations._simu import _Simu, Load_Simu
-from .simulations._utils import Save_Force_Displacement, Load_Force_Displacement
+from .simulations._utils import Save_Force_Displacement, Load_Force_Displacement, Save_pickle, Load_pickle
 
 # ----------------------------------------------
 # Elastic

{easyfea-1.2.2 → easyfea-1.2.4}/EasyFEA/fem/_gmsh_interface.py

@@ -441,7 +441,7 @@ class Mesher:
                 p1 = factory.addPoint(*line.pt1.coord, line.meshSize) # start
                 p2 = factory.addPoint(*line.pt2.coord, line.meshSize) # end
                 p3 = factory.addPoint(*line.pt3.coord, line.meshSize) # mid
-                geom_line = factory.addCircleArc(p1, p3, p2)
+                geom_line = factory.addCircleArc(p1, p3, p2, center=False)
             else:
                 raise Exception("You need to give lines or arcs.")
             factory.fragment(ents, [(1, geom_line)], removeTool=False)
@@ -879,6 +879,7 @@ class Mesher:
             elif isinstance(crack, Contour):  # 2D CRACK
 
                 loop, lines, points, openLns, openPts = self._Loop_From_Contour(crack)
+                
                 try:
                     surf = self._Surface_From_Loops([loop])
                 except Exception:
@@ -925,7 +926,7 @@ class Mesher:
         newEntities.extend([(2, surf) for surf in entities2D])
 
         if factory == gmsh.model.occ:
-            o, m = gmsh.model.occ.fragment(entities, newEntities)
+            o, m = gmsh.model.occ.fragment(entities, newEntities, removeTool=False)
 
         self._Synchronize() # mandatory
 
@@ -1315,6 +1316,8 @@ class Mesher:
             Geometric objects to refine de background mesh
         meshSize : float
             size of elements outside the domain
+        extrude : list[float]
+            extrusion vector, by default [0,0,1]
         """
 
         # See t10.py in the gmsh tutorials
@@ -1323,8 +1326,8 @@ class Mesher:
         assert isinstance(refineGeoms, Iterable), "refineGeoms must be a list of geometric objects."
 
         if refineGeoms is None or len(refineGeoms) == 0: return
-
-        fields = []
+        
+        fields = list(gmsh.model.mesh.field.list())
 
         for geom in refineGeoms:
 
@@ -1466,6 +1469,7 @@ class Mesher:
         """
         
         self._Set_mesh_algorithm(elemType)
+        # self._factory.removeAllDuplicates()
         self._Synchronize() # mandatory
 
         tic = Tic()

{easyfea-1.2.2 → easyfea-1.2.4}/EasyFEA/materials/_elastic.py

@@ -4,6 +4,7 @@
 
 from abc import ABC, abstractmethod
 from typing import Union
+from scipy.linalg import sqrtm
 
 # utilities
 from .. import np
@@ -106,6 +107,7 @@ class _Elas(_IModel, ABC):
         shape = (3, 3) if self.dim == 2 else (6, 6)
         assert array.shape[-2:] == shape, f"With dim = {self.dim} array must be a {shape} matrix"
         self.__C = array
+        self.__sqrt_C = None # dont remove
 
     @property
     def isHeterogeneous(self) -> bool:
@@ -129,6 +131,7 @@ class _Elas(_IModel, ABC):
         shape = (3, 3) if self.dim == 2 else (6, 6)
         assert array.shape[-2:] == shape, f"With dim = {self.dim} array must be a {shape} matrix"
         self.__S = array
+        self.__sqrt_S = None # dont remove
 
     @abstractmethod
     def Walpole_Decomposition(self) -> tuple[np.ndarray, np.ndarray]:
@@ -136,6 +139,62 @@ class _Elas(_IModel, ABC):
         C = sum(ci * Ei).\n        
         returns ci, Ei"""
         return np.array([]), np.array([])
+    
+    def Get_sqrt_C_S(self) -> tuple[np.ndarray, np.ndarray]:
+        """Returns the Matrix square root of C and S."""
+
+        C = self.C
+
+        try:
+            self.__sqrt_C is None
+            self.__sqrt_S is None
+        except AttributeError:
+            # init
+            self.__sqrt_C = None
+            self.__sqrt_S = None
+
+        if self.__sqrt_C is None:
+
+            if self.isHeterogeneous:
+
+                shape = C.shape
+                
+                assert len(shape) == 3, "This function is not currently implemented for heterogeneous matrices where material properties are defined on Gauss points."
+                
+                uniq_C, inverse = np.unique(C, return_inverse=True, axis=0)
+
+                sqrtC = np.zeros_like(C, dtype=float)
+                sqrtS = np.zeros_like(C, dtype=float)
+
+                for i, C in enumerate(uniq_C):
+
+                    elems = np.where(inverse == i)[0]
+
+                    sqrtmC = sqrtm(C)
+                    sqrtC[elems] = sqrtmC
+
+                    sqrtmS = np.linalg.inv(sqrtmC)
+                    sqrtS[elems] = sqrtmS
+                
+            else:
+                sqrtC = sqrtm(C)
+                
+                sqrtS = np.linalg.inv(sqrtC) # faster than sqrtm(self.S)
+                # sqrtS = sqrtm(self.S)
+
+                # # give the same results !!!
+                # test = np.linalg.norm(sqrtS - sqrtm(self.S))/np.linalg.norm(sqrtS)
+                # assert test < 1e-12
+
+            self.__sqrt_C = sqrtC
+            self.__sqrt_S = sqrtS
+        
+        else:
+
+            sqrtC = self.__sqrt_C.copy()
+            sqrtS = self.__sqrt_S.copy()
+
+        return sqrtC, sqrtS
 
 # ----------------------------------------------
 # Isotropic
@@ -747,8 +806,17 @@ class Elas_Anisot(_Elas):
         shape = C.shape
         assert (shape[-2], shape[-1]) in [(3,3), (6,6)], 'C must be a (3,3) or (6,6) matrix'        
         dim = 3 if C.shape[-1] == 6 else 2
-        testSym = np.linalg.norm(C.T - C)/np.linalg.norm(C)
-        assert testSym <= 1e-12, "The matrix is not symmetrical."
+        if len(C.shape) == 2:
+            Ct = C.T
+        elif len(C.shape) == 3:
+            Ct = np.transpose(C, (0,2,1))
+        elif len(C.shape) == 4:
+            Ct = np.transpose(C, (0,1,3,2))
+        else:
+            raise ValueError("This matrix must be of dimensions (dim, dim), (Ne, dim, dim) or (Ne, nPg, dim, dim).")
+
+        testSym = np.linalg.norm(Ct - C, axis=(-2, -1))/np.linalg.norm(C, axis=(-2, -1))
+        assert np.max(testSym) <= 1e-12, "The matrix is not symmetrical."
 
         if useVoigtNotation:
             C_mandel = KelvinMandel_Matrix(dim, C)

{easyfea-1.2.2 → easyfea-1.2.4}/EasyFEA/materials/_phasefield.py

@@ -5,8 +5,6 @@
 from typing import Union
 from enum import Enum
 
-from scipy.linalg import sqrtm
-
 # utilities
 from .. import np
 from ..utilities import Tic, Numba_Interface
@@ -737,22 +735,28 @@ class PhaseField(_IModel):
         # Here the material is supposed to be homogeneous
         material = self.__material
 
-        C = material.C        
-        
-        assert not material.isHeterogeneous, "He decomposition has not been implemented for heterogeneous materials"
-        # for heterogeneous materials how to make sqrtm ?
-        # TODO obtain the eigenvalues and eigenvectors of the 2D case with the 3D spectral decomposition?
-        sqrtC = sqrtm(C)
+        C = material.C
+
+        tic = Tic()
+        sqrtC, inv_sqrtC = material.Get_sqrt_C_S()
+        # inv(sqrtC) = sqrtS
+        tic.Tac("Split",f"sqrt C and S", False)
         
-        if verif :
+        if verif:
             # checks that C^1/2 * C^1/2 = C
-            testC = sqrtC @ sqrtC - C
-            assert np.linalg.norm(testC)/np.linalg.norm(C) < 1e-12
-
-        inv_sqrtC = np.linalg.inv(sqrtC)
+            diff_C = sqrtC @ sqrtC - C
+            test_C = np.linalg.norm(diff_C, axis=(-2,-1))/np.linalg.norm(C, axis=(-2,-1))
+            assert np.max(test_C) < 1e-12
 
         # computes new "strain" field
-        Epsilont_e_pg = np.einsum('ij,epj->epi', sqrtC, Epsilon_e_pg, optimize='optimal')
+        
+        ind = ""
+        if len(C.shape) == 3:
+            ind = 'e'
+        elif len(C.shape) == 4:
+            ind = 'ep'
+
+        Epsilont_e_pg = np.einsum(f'{ind}ij,epj->epi', sqrtC, Epsilon_e_pg, optimize='optimal')
 
         # Computes projectors
         projPt_e_pg, projMt_e_pg = self.__Spectral_Decomposition(Epsilont_e_pg, verif)
@@ -763,18 +767,23 @@ class PhaseField(_IModel):
         # projM_e_pg = inv_sqrtC @ (projMt_e_pg @ sqrtC)
         
         # faster
-        projP_e_pg = np.einsum('ij,epjk,kl->epil', inv_sqrtC, projPt_e_pg, sqrtC, optimize='optimal')
-        projM_e_pg = np.einsum('ij,epjk,kl->epil', inv_sqrtC, projMt_e_pg, sqrtC, optimize='optimal')
+        projP_e_pg = np.einsum(f'{ind}ij,epjk,{ind}kl->epil', inv_sqrtC, projPt_e_pg, sqrtC, optimize='optimal')
+        projM_e_pg = np.einsum(f'{ind}ij,epjk,{ind}kl->epil', inv_sqrtC, projMt_e_pg, sqrtC, optimize='optimal')
 
         tic.Tac("Split",f"proj Tild to proj", False)
 
-        cP_e_pg = sqrtC @ (projPt_e_pg @ sqrtC)
-        cM_e_pg = sqrtC @ (projMt_e_pg @ sqrtC)
-        # cP_e_pg = C @ projP_e_pg
-        # cM_e_pg = C @ projM_e_pg
+        if material.isHeterogeneous:
+            cP_e_pg = np.einsum(f'{ind}ij,epjk->epik', C, projP_e_pg, optimize='optimal')
+            cM_e_pg = np.einsum(f'{ind}ij,epjk->epik', C, projM_e_pg, optimize='optimal')
+
+        else:
+            # cP_e_pg = sqrtC @ (projPt_e_pg @ sqrtC)
+            # cM_e_pg = sqrtC @ (projMt_e_pg @ sqrtC)
+            cP_e_pg = C @ projP_e_pg
+            cM_e_pg = C @ projM_e_pg
 
-        # test_cP = sqrtC @ (projPt_e_pg @ sqrtC) - C @ projP_e_pg
-        # test_cP = sqrtC @ (projMt_e_pg @ sqrtC) - C @ projM_e_pg
+            # test_cP = sqrtC @ (projPt_e_pg @ sqrtC) - C @ projP_e_pg
+            # test_cP = sqrtC @ (projMt_e_pg @ sqrtC) - C @ projM_e_pg
 
         tic.Tac("Split",f"cP_e_pg and cM_e_pg", False)
 

{easyfea-1.2.2 → easyfea-1.2.4}/EasyFEA/simulations/_dic.py

@@ -189,12 +189,12 @@ class DIC(_IObserver):
         tic.Tac("DIC", "ROI", self._verbosity)
     
     @property
-    def roi(self) -> np.ndarray[bool]:
+    def roi(self) -> np.ndarray:
         """roi as a vector."""
         return self.__roi.copy()
 
     @property
-    def ROI(self) -> np.ndarray[bool]:
+    def ROI(self) -> np.ndarray:
         """roi as a matrix."""
         return self.roi.reshape(self.shape)
 

{easyfea-1.2.2 → easyfea-1.2.4}/EasyFEA/simulations/_phasefield.py

@@ -210,12 +210,12 @@ class PhaseFieldSimu(_Simu):
         Returns
         -------
         np.ndarray, np.ndarray, csr_matrix, bool
-            u_np1, d_np1, Kglob, converged
+            u_np1, d_np1, Ku, converged
 
             such that:\n
             - u_np1: displacement vector field\n
             - d_np1: damage scalar field\n
-            - Kglob: displacement stiffness matrix\n
+            - Ku: displacement stiffness matrix\n
             - converged: the solution has converged\n
         """
 
@@ -311,9 +311,9 @@ class PhaseFieldSimu(_Simu):
         self.__convIter = convIter
         self.__timeIter = timeIter
 
-        Kglob = self.__Ku.copy()
+        Ku = self.__Ku.copy()
             
-        return u_np1, d_np1, Kglob, converged
+        return u_np1, d_np1, Ku, converged
 
     # ------------------------------------------- Elastic problem -------------------------------------------
 
@@ -509,7 +509,7 @@ class PhaseFieldSimu(_Simu):
         # Calculating elementary matrix
         Kd_e, Fd_e = self.__Construct_Damage_Matrix()
 
-        # Assemblage
+        # Assembly
         tic = Tic()        
 
         self.__Kd = sparse.csr_matrix((Kd_e.ravel(), (linesScalar_e, columnsScalar_e)), shape = (Ndof, Ndof))
@@ -669,16 +669,16 @@ class PhaseFieldSimu(_Simu):
         
         return self.Results_Reshape_values(values, nodeValues)
 
-    def __indexResult(self, resultat: str) -> int:
+    def __indexResult(self, result: str) -> int:
 
         dim = self.dim
 
-        if len(resultat) <= 2:
-            if "x" in resultat:
+        if len(result) <= 2:
+            if "x" in result:
                 return 0
-            elif "y" in resultat:
+            elif "y" in result:
                 return 1
-            elif "z" in resultat:
+            elif "z" in result:
                 return 1
 
     def _Calc_Psi_Elas(self) -> float:

{easyfea-1.2.2 → easyfea-1.2.4}/EasyFEA/simulations/_utils.py

@@ -12,7 +12,7 @@ import pickle
 # ----------------------------------------------
 
 def Save_Force_Displacement(force: np.ndarray, displacement: np.ndarray, folder:str):
-    """Save the values of load and displacements in the folder"""
+    """Saves the values of force and displacement in the folder"""
     
     folder_PythonEF = Folder.Get_Path(Folder.Get_Path())
     filename = Folder.Join(folder, "force-displacement.pickle")
@@ -31,7 +31,7 @@ def Save_Force_Displacement(force: np.ndarray, displacement: np.ndarray, folder:
     Display.MyPrint(f'Saved:\n{filename.replace(folder_PythonEF,"")}\n','green')
     
 def Load_Force_Displacement(folder:str):
-    """Loads forces and displacements.
+    """Loads force and displacement.
 
     Parameters
     ----------
@@ -56,4 +56,41 @@ def Load_Force_Displacement(folder:str):
     
     Display.MyPrint(f'Loaded:\n{shortName}\n','green')
 
-    return force, displacement
+    return force, displacement
+
+# ----------------------------------------------
+# Save obj in pickle file
+# ----------------------------------------------
+
+def Save_pickle(obj, folder: str, filename: str) -> None:
+    """Saves the object in folder/filename.pickle."""
+    
+    folder_PythonEF = Folder.Get_Path(Folder.Get_Path())
+    
+    file = Folder.Join(folder, f"{filename}.pickle")
+
+    if not Folder.os.path.exists(folder):
+        Folder.os.makedirs(folder)
+
+    with open(file, "wb") as f:
+        pickle.dump(obj, f)
+    
+    Display.MyPrint(f'Saved:\n{file.replace(folder_PythonEF,"")}\n','green')
+
+def Load_pickle(folder:str, filename: str):
+    """Returns folder/filename.pickle object."""
+
+    folder_PythonEF = Folder.Get_Path(Folder.Get_Path())
+
+    file = Folder.Join(folder, f"{filename}.pickle")
+
+    shortName = file.replace(folder_PythonEF,'') 
+    error = f"{shortName} does not exist"
+    assert Folder.Exists(file), error
+
+    with open(file, 'rb') as f:
+        obj = pickle.load(f)
+    
+    Display.MyPrint(f'Loaded:\n{shortName}\n','green')
+
+    return obj

{easyfea-1.2.2 → easyfea-1.2.4}/EasyFEA/utilities/Display.py

@@ -43,7 +43,7 @@ def Plot_Result(obj, result: Union[str,np.ndarray], deformFactor=0.0, coef=1.0,
         object containing the mesh
     result : str or np.ndarray
         result you want to display.\n
-        Must be included in simu.Get_Results() or be a numpy array of size of (Ne, Ne).
+        Must be included in simu.Get_Results() or be a numpy array of size of (Nn, Ne).
     deformFactor : float, optional
         factor used to display the deformed solution (0 means no deformations), default 0.0
     coef : float, optional
@@ -133,7 +133,7 @@ def Plot_Result(obj, result: Union[str,np.ndarray], deformFactor=0.0, coef=1.0,
         levels = np.linspace(min, max, ncolors)
     
     if ncolors != 256:
-        norm = colors.BoundaryNorm(boundaries=np.linspace(min, max, 11), ncolors=256)
+        norm = colors.BoundaryNorm(boundaries=levels, ncolors=256)
     else:
         norm = None
 
@@ -178,7 +178,7 @@ def Plot_Result(obj, result: Union[str,np.ndarray], deformFactor=0.0, coef=1.0,
             pc.set_clim(min, max)
             pc.set_array(values)
             ax.add_collection(pc)
-            ticks = None if ncolors != 11 else ticks
+            # ticks = None if ncolors != 11 else ticks
 
         # Plots node values
         elif mesh.Nn == len(values):

{easyfea-1.2.2 → easyfea-1.2.4}/EasyFEA/utilities/Folder.py

@@ -104,7 +104,10 @@ def PhaseField_Folder(folder: str, material: str,
 
     if nL != 0:
         assert nL > 0
-        name = f"{name} nL={nL}"
+        if isinstance(nL, float):
+            name = f"{name} nL={nL:.2f}"
+        else:
+            name = f"{name} nL={nL}"
 
     workFolder = New_File(folder, results=True)
     path = workFolder.split(folder)[0]
@@ -114,9 +117,4 @@ def PhaseField_Folder(folder: str, material: str,
     else:
         workFolder = Join(folder, name)
 
-    path_PythonEf = Get_Path(Get_Path())
-    
-    from . import Display
-    Display.MyPrint(workFolder.replace(path_PythonEf,""), 'green')
-
     return Join(path, workFolder)

{easyfea-1.2.2 → easyfea-1.2.4}/EasyFEA/utilities/PyVista_Interface.py

@@ -2,7 +2,8 @@
 # This file is part of the EasyFEA project.
 # EasyFEA is distributed under the terms of the GNU General Public License v3 or later, see LICENSE.txt and CREDITS.md for more information.
 
-"""Module providing an interface with PyVista (https://docs.pyvista.org/version/stable/)."""
+"""Module providing an interface with PyVista (https://docs.pyvista.org/version/stable/).\n
+https://docs.pyvista.org/api/plotting/plotting.html"""
 
 from typing import Union, Callable
 from cycler import cycler
@@ -21,7 +22,7 @@ def Plot(obj, result: Union[str,np.ndarray]=None, deformFactor=0.0, coef=1.0, no
                 color=None, show_edges=False, edge_color='k', line_width=None,
                 show_vertices=False, point_size=None, opacity=1.0,
                 style='surface', cmap="jet", n_colors=256, clim=None,
-                plotter: pv.Plotter=None, show_grid=False, verticalColobar=True, **kwargs):
+                plotter: pv.Plotter=None, show_grid=False, colorbarTitle=None, verticalColobar=True, **kwargs):
     """Plots the object obj that can be either a simu, mesh, MultiBlock, PolyData.\n
     If you want to plot the solution use plotter.show().
 
@@ -64,6 +65,8 @@ def Plot(obj, result: Union[str,np.ndarray]=None, deformFactor=0.0, coef=1.0, no
         The pyvista plotter, by default None and create a new Plotter instance
     show_grid : bool, optionnal
         Show the grid, by default False
+    colorbarTitle: str, optionnal
+        colorbar title, by default None
     verticalColobar : bool, optionnal
         color bar is vertical, by default True
     **kwargs:
@@ -95,6 +98,8 @@ def Plot(obj, result: Union[str,np.ndarray]=None, deformFactor=0.0, coef=1.0, no
     name = None if pvMesh.n_arrays == 0 else name
     if name != None:
         pvMesh[name] *= coef
+    
+    colorbarTitle = name if colorbarTitle is None else colorbarTitle
 
     if plotter is None:        
         plotter = _Plotter()
@@ -120,7 +125,7 @@ def Plot(obj, result: Union[str,np.ndarray]=None, deformFactor=0.0, coef=1.0, no
                         show_vertices=show_vertices, point_size=point_size,
                         opacity=opacity,
                         style=style, cmap=cmap, n_colors=n_colors, clim=clim,
-                        scalar_bar_args={'title': name, 'vertical': verticalColobar, pos: val},
+                        scalar_bar_args={'title': colorbarTitle, 'vertical': verticalColobar, pos: val},
                         **kwargs)
 
     _setCameraPosition(plotter, inDim)
@@ -161,7 +166,7 @@ def Plot_Mesh(obj, deformFactor=0.0, opacity=1.0, color='cyan', edge_color='blac
     return plotter
 
 def Plot_Nodes(obj, nodes: np.ndarray=None, showId=False, deformFactor=0, color='red',
-               folder="", plotter: pv.Plotter=None):
+               folder="", label=None, plotter: pv.Plotter=None):
     """Plots mesh's nodes.
 
     Parameters
@@ -175,7 +180,9 @@ def Plot_Nodes(obj, nodes: np.ndarray=None, showId=False, deformFactor=0, color=
     deformFactor : float, optional
         Factor used to display the deformed solution (0 means no deformations), default 0.0
     color : str, optional
-        color, default 'red'    
+        color, default 'red'
+    label : str, optional
+        label, by default None
     plotter : pv.Plotter, optional
         The pyvista plotter, by default None and create a new Plotter instance
 
@@ -215,11 +222,13 @@ def Plot_Nodes(obj, nodes: np.ndarray=None, showId=False, deformFactor=0, color=
         pvData["myLabels"] = myLabels
         plotter.add_point_labels(pvData, "myLabels", point_color=color, render_points_as_spheres=True)
     else:
-        plotter.add_mesh(pvData, color=color, render_points_as_spheres=True)
+        plotter.add_mesh(pvData, color=color, label=label, render_points_as_spheres=True)
 
     return plotter
 
-def Plot_Elements(obj, nodes: np.ndarray=None, dimElem: int=None, showId=False, deformFactor=0, opacity=1.0, color='red', edge_color='black', line_width=None, plotter: pv.Plotter=None):
+def Plot_Elements(obj, nodes: np.ndarray=None, dimElem: int=None, showId=False,
+                  deformFactor=0, opacity=1.0, color='red', edge_color='black',
+                  line_width=None, label=None, plotter: pv.Plotter=None):
     """Plots the mesh elements corresponding to the given nodes.
 
     Parameters
@@ -242,6 +251,8 @@ def Plot_Elements(obj, nodes: np.ndarray=None, dimElem: int=None, showId=False,
         color used to display segments, by default 'black'    
     line_width : float, optional
         Thickness of lines, by default None
+    label : str, optional
+        label, by default None
     plotter : pv.Plotter, optional
         The pyvista plotter, by default None and create a new Plotter instance
 
@@ -264,7 +275,8 @@ def Plot_Elements(obj, nodes: np.ndarray=None, dimElem: int=None, showId=False,
             return
 
     if plotter == None:
-        plotter = Plot(obj, deformFactor=deformFactor, style='wireframe', color=edge_color, line_width=line_width)
+        # plotter = Plot(obj, deformFactor=deformFactor, style='wireframe', color=edge_color, line_width=line_width)
+        plotter = _Plotter()
     
     for groupElem in mesh.Get_list_groupElem(dimElem):
 
@@ -281,7 +293,7 @@ def Plot_Elements(obj, nodes: np.ndarray=None, dimElem: int=None, showId=False,
 
         pvGroup = _pvGrid(newGroupElem)
 
-        Plot(pvGroup, opacity=opacity, color=color, edge_color=edge_color, plotter=plotter, line_width=line_width)
+        Plot(pvGroup, opacity=opacity, color=color, edge_color=edge_color, plotter=plotter, line_width=line_width, label=label)
 
         if showId:
             centers = np.mean(coordo[groupElem.connect[elements]], axis=1)
@@ -622,11 +634,11 @@ def _Plotter(off_screen=False, add_axes=True, shape=(1,1), linkViews=True):
     plotter.subplot(0,0)
     return plotter
 
-def _setCameraPosition(plotter: pv.Plotter, inDim: int):
+def _setCameraPosition(plotter: pv.Plotter, inDim: int, elevation=25, azimuth=10):
     plotter.camera_position = 'xy'
     if inDim == 3:
-        plotter.camera.elevation += 25
-        plotter.camera.azimuth += 10
+        plotter.camera.elevation = elevation
+        plotter.camera.azimuth = azimuth
         plotter.camera.reset_clipping_range()
     # if inDim == 3:        
     #     plotter.camera.elevation += 15

{easyfea-1.2.2 → easyfea-1.2.4}/EasyFEA.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: EasyFEA
-Version: 1.2.2
+Version: 1.2.4
 Summary: User-friendly Python library that simplifies finite element analysis.
 Author-email: Matthieu Noel <matthieu.noel7@gmail.com>
 License:                     GNU GENERAL PUBLIC LICENSE

{easyfea-1.2.2 → easyfea-1.2.4}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: EasyFEA
-Version: 1.2.2
+Version: 1.2.4
 Summary: User-friendly Python library that simplifies finite element analysis.
 Author-email: Matthieu Noel <matthieu.noel7@gmail.com>
 License:                     GNU GENERAL PUBLIC LICENSE

{easyfea-1.2.2 → easyfea-1.2.4}/pyproject.toml

@@ -4,7 +4,7 @@ build-backend = "setuptools.build_meta"
 
 [project]
 name = "EasyFEA"
-version = "1.2.2"
+version = "1.2.4"
 description = "User-friendly Python library that simplifies finite element analysis."
 authors = [
     {name = "Matthieu Noel", email = "matthieu.noel7@gmail.com"}