elasticipy 6.0.1__tar.gz → 6.1.0__tar.gz

elasticipy-6.1.0/.github/FUNDING.yml

@@ -0,0 +1 @@
+ko_fi: doriandepriester

{elasticipy-6.0.1 → elasticipy-6.1.0}/.github/workflows/Codecov.yml

@@ -34,7 +34,7 @@ jobs:
         env:
           MP_API_KEY: ${{ secrets.MP_TOKEN }}
 
-        run: pytest tests
+        run: pytest --cov=elasticipy --cov-report xml:coverage.xml --cov-report term tests
 
       - name: Upload results to Codecov
         uses: codecov/codecov-action@v5

{elasticipy-6.0.1/src/elasticipy.egg-info → elasticipy-6.1.0}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: elasticipy
-Version: 6.0.1
+Version: 6.1.0
 Summary: A Python library for elasticity tensor computations
 Author-email: Dorian Depriester <dorian.dep@gmail.com>
 License: MIT
@@ -25,6 +25,7 @@ Requires-Dist: matplotlib
 Requires-Dist: qtpy
 Requires-Dist: pyqt5
 Requires-Dist: pandas
+Requires-Dist: orix
 Provides-Extra: dev
 Requires-Dist: pytest; extra == "dev"
 Requires-Dist: pytest-cov; extra == "dev"
@@ -44,7 +45,7 @@ Dynamic: license-file
 [![codecov](https://codecov.io/gh/DorianDepriester/Elasticipy/graph/badge.svg?token=VUZPEUPBH1)](https://codecov.io/gh/DorianDepriester/Elasticipy)
 ![PyPI - Python Version](https://img.shields.io/pypi/pyversions/Elasticipy)
 [![DOI](https://joss.theoj.org/papers/10.21105/joss.07940/status.svg)](https://doi.org/10.21105/joss.07940)
-[![Binder](https://mybinder.org/badge_logo.svg)](https://mybinder.org/v2/gh/DorianDepriester/Elasticipy/HEAD?urlpath=%2Fdoc%2Ftree%2FElasticipy_for_the_Impatient.ipynb)
+[![Binder](https://mybinder.org/badge_logo.svg)](https://mybinder.org/v2/gh/DorianDepriester/elasticipy-notebook/HEAD?urlpath=%2Fdoc%2Ftree%2FElasticipy.ipynb)
 
 
 # ![Elasticipy](https://raw.githubusercontent.com/DorianDepriester/Elasticipy/refs/heads/main/docs/source/logo/logo_text.svg)
@@ -82,7 +83,7 @@ Tutorials and full documentation are available on [ReadTheDoc](https://elasticip
 
 ## ⏱️ Elasticipy in a nutshell
 Take a 5-minute tour through Elasticipy's main features by running the online Jupyter Notebook, hosted on 
-[Binder](https://mybinder.org/v2/gh/DorianDepriester/Elasticipy/HEAD?urlpath=%2Fdoc%2Ftree%2FElasticipy_for_the_Impatient.ipynb).
+[Binder](https://mybinder.org/v2/gh/DorianDepriester/elasticipy-notebook/HEAD?urlpath=%2Fdoc%2Ftree%2FElasticipy.ipynb).
  
 
 ## 🔍 Sources
@@ -113,7 +114,8 @@ You can use the following BibTeX entry:
     number = {115}, 
     pages = {7940}, 
     author = {Depriester, Dorian and Kubler, Régis}, 
-    title = {elasticipy: A Python package for linear elasticity and tensor analysis}, 
+    title = {Elasticipy: A Python package for linear elasticity and tensor analysis}, 
     journal = {Journal of Open Source Software}
 }
 ````
+Alternatively, you can by me a coffee on [![ko-fi](https://ko-fi.com/img/githubbutton_sm.svg)](https://ko-fi.com/W7W11S5TCH)

{elasticipy-6.0.1 → elasticipy-6.1.0}/README.md

@@ -9,7 +9,7 @@
 [![codecov](https://codecov.io/gh/DorianDepriester/Elasticipy/graph/badge.svg?token=VUZPEUPBH1)](https://codecov.io/gh/DorianDepriester/Elasticipy)
 ![PyPI - Python Version](https://img.shields.io/pypi/pyversions/Elasticipy)
 [![DOI](https://joss.theoj.org/papers/10.21105/joss.07940/status.svg)](https://doi.org/10.21105/joss.07940)
-[![Binder](https://mybinder.org/badge_logo.svg)](https://mybinder.org/v2/gh/DorianDepriester/Elasticipy/HEAD?urlpath=%2Fdoc%2Ftree%2FElasticipy_for_the_Impatient.ipynb)
+[![Binder](https://mybinder.org/badge_logo.svg)](https://mybinder.org/v2/gh/DorianDepriester/elasticipy-notebook/HEAD?urlpath=%2Fdoc%2Ftree%2FElasticipy.ipynb)
 
 
 # ![Elasticipy](https://raw.githubusercontent.com/DorianDepriester/Elasticipy/refs/heads/main/docs/source/logo/logo_text.svg)
@@ -47,7 +47,7 @@ Tutorials and full documentation are available on [ReadTheDoc](https://elasticip
 
 ## ⏱️ Elasticipy in a nutshell
 Take a 5-minute tour through Elasticipy's main features by running the online Jupyter Notebook, hosted on 
-[Binder](https://mybinder.org/v2/gh/DorianDepriester/Elasticipy/HEAD?urlpath=%2Fdoc%2Ftree%2FElasticipy_for_the_Impatient.ipynb).
+[Binder](https://mybinder.org/v2/gh/DorianDepriester/elasticipy-notebook/HEAD?urlpath=%2Fdoc%2Ftree%2FElasticipy.ipynb).
  
 
 ## 🔍 Sources
@@ -78,7 +78,8 @@ You can use the following BibTeX entry:
     number = {115}, 
     pages = {7940}, 
     author = {Depriester, Dorian and Kubler, Régis}, 
-    title = {elasticipy: A Python package for linear elasticity and tensor analysis}, 
+    title = {Elasticipy: A Python package for linear elasticity and tensor analysis}, 
     journal = {Journal of Open Source Software}
 }
 ````
+Alternatively, you can by me a coffee on [![ko-fi](https://ko-fi.com/img/githubbutton_sm.svg)](https://ko-fi.com/W7W11S5TCH)

{elasticipy-6.0.1 → elasticipy-6.1.0}/docs/source/API/API.rst

@@ -4,12 +4,13 @@ API documentation
 .. toctree::
    :maxdepth: 5
 
+   Elasticipy.crystal_texture
+   Elasticipy.interfaces.PRISMS
+   Elasticipy.interfaces.FEPX
+   Elasticipy.plasticity
+   Elasticipy.spherical_function
    Elasticipy.tensors.elasticity
    Elasticipy.tensors.fourth_order
    Elasticipy.tensors.second_order
    Elasticipy.tensors.stress_strain
    Elasticipy.tensors.thermal_expansion
-   Elasticipy.plasticity
-   Elasticipy.spherical_function
-   Elasticipy.interfaces.PRISMS
-   Elasticipy.interfaces.FEPX

elasticipy-6.1.0/docs/source/API/Elasticipy.crystal_texture.rst

@@ -0,0 +1,8 @@
+Elasticipy.crystal_texture
+======================================
+
+.. automodule:: elasticipy.crystal_texture
+   :members:
+   :undoc-members:
+   :show-inheritance:
+   :inherited-members:

{elasticipy-6.0.1 → elasticipy-6.1.0}/docs/source/Tutorials/FourthOrderTensors.rst

@@ -189,7 +189,7 @@ representation. Indeed, one can check that their components as full representati
 
 Actually, a more direct way to check this equality is:
 
-    >>> I_kelvin == I_voigt
+    >>> print(I_kelvin == I_voigt)
     True
 
 As said above, the two tensors above can `mixed`, although they are not represented the same way. E.g.:

{elasticipy-6.0.1 → elasticipy-6.1.0}/docs/source/Tutorials/Tutorial_AveragingMethods.rst

@@ -59,6 +59,7 @@ As a comparison, let's see how the underlying assumption impair the Young moduli
     >>> print(E_Voigt.mean(), E_Reuss.mean(), E_Hill.mean())
     145.66862361382908 76.13802022165395 111.67690532486962
 
+.. _rotations:
 
 Finite number of orientations
 =============================
@@ -66,6 +67,11 @@ If one wants to consider a finite set of orientations, or consider that these or
 over the SO3 (i.e. if we have a crystallographic texture), the first steps consists in estimating the mean of all the
 rotated stiffnesses/compliances (depending on the model).
 
+.. note::
+
+    Elasticipy now handles *textures* to compute the "exact" averages (i.e. based on sampling). See :ref:`textures`.
+
+
 As an example, we assume that the material displays a perfect fiber texture along **z**. In terms of Bunge-Euler angles,
 it means that we have Phi=0 for each orientation:
 
@@ -93,6 +99,11 @@ Let's see how the Young modulus is distributed over space:
 
 It is thus clear that the fiber along the *z* axis results in a transverse-isotropic behavior in the X-Y plane.
 
+Actually, a quicker way to compute the average from a finite set of orientations is to pass them as the averaging
+function, e.g.:
+
+    >>> Chill = C.Hill_average(C, orientations=rotations)
+
 Plotting Voigt, Reuss and Hill averages at once
 ===============================================
 
@@ -105,22 +116,19 @@ sections (see :ref:`here<plotting>` for details) for each of the aforementioned
     from scipy.spatial.transform import Rotation
     import numpy as np
     from elasticipy.tensors.elasticity import StiffnessTensor
+    from elasticipy.crystal_texture import FibreTexture
+
     C = StiffnessTensor.monoclinic(phase_name='TiNi',
                                     C11=231, C12=127, C13=104,
                                     C22=240, C23=131, C33=175,
                                     C44=81, C55=11, C66=85,
                                     C15=-18, C25=1, C35=-3, C46=3)
-    phi1 = np.random.random(10000)*2*np.pi # Random sampling from 0 to 2pi
-    Phi = phi2 = np.zeros(10000)
-    Euler_angles = np.array([phi1, Phi, phi2]).T
-    rotations = Rotation.from_euler('ZXZ', Euler_angles)    # Bunge-Euler angles
-    C_rotated = C * rotations
+    texture = FibreTexture.from_Euler(Phi=0, phi2=0)
     for method in ['Reuss', 'Hill', 'Voigt']:
-         C_avg = C_rotated.average(method)
+         C_avg = C.average(method, orientations = texture)
          if method == 'Reuss':
              fig, axs = C_avg.Young_modulus.plot_xyz_sections(label='Reuss') # Create fig and axes (sections)
          else:
              fig, axs = C_avg.Young_modulus.plot_xyz_sections(fig=fig, axs=axs, label=method) # Use existing axes
     axs[-1].legend()
-    fig.tight_layout()
     fig.show()

{elasticipy-6.0.1 → elasticipy-6.1.0}/docs/source/Tutorials/Tutorial_StiffnessTensor.rst

@@ -40,7 +40,7 @@ shear and bulk moduli
     Hyperspherical function
     Min=83999.99999999991, Max=84000.00000000007
 
-    >>> C.bulk_modulus
+    >>> print(C.bulk_modulus)
     140000.0
 
 .. note::
@@ -51,7 +51,7 @@ shear and bulk moduli
 One can check that both approaches yield the same tensor:
 
     >>> C2 = StiffnessTensor.isotropic(G=84e3, K=140e3)
-    >>> C2 == C
+    >>> print(C2 == C)
     True
 
 
@@ -84,7 +84,7 @@ let's see its value along the x, y and z directions:
     >>> Ey = E.eval([0,1,0])
     >>> Ez = E.eval([0,0,1])
     >>> print((Ex, Ey, Ez))
-    (124.52232440357189, 120.92120854784433, 96.13750721721384)
+    (np.float64(124.52232440357189), np.float64(120.92120854784433), np.float64(96.13750721721384))
 
 .. note::
 
@@ -113,7 +113,7 @@ deviation of the Young modulus:
 
 Another way to evidence anisotropy is to use the universal anisotropy factor [Ranganathan]_:
 
-    >>> C.universal_anisotropy
+    >>> print(C.universal_anisotropy)
     5.141009551641412
 
 Shear moduli and Poisson ratios

{elasticipy-6.0.1 → elasticipy-6.1.0}/docs/source/Tutorials/Tutorial_StressStrain.rst

@@ -141,7 +141,7 @@ that the maximum shear stress is around 0.849. Those can be checked by:
 
     >>> s.principal_stresses()
     array([ 1.51167769,  0.3519979 , -0.18566326])
-    >>> s.Tresca() / 2
+    >>> print(s.Tresca() / 2)
     0.848670477704235
 
 

elasticipy-6.1.0/docs/source/Tutorials/texture.rst

@@ -0,0 +1,188 @@
+.. _textures:
+
+Working with crystallographic textures
+---------------------------------------
+
+With the help of `orix <https://orix.readthedocs.io/en/stable/index.html>`_, Elasticipy allows to compute averages based
+on crystallographic texture (in addition to single orientations; see :ref:`rotations`).
+
+Define and compose textures
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+Discrete textures
+=================
+A series of "usual" texture components (e.g. Goss, cube, Brass etc.) are already implemented in Elasticipy:
+
+.. doctest::
+
+    >>> from elasticipy.crystal_texture import DiscreteTexture
+    >>> goss = DiscreteTexture.Goss()
+    >>> goss
+    Crystallographic texture
+    φ1=0.00°, ϕ=45.00°, φ2=0.00°
+
+This texture actually consists in a single orientation (as opposed to fibre textures, see below). It can be used to
+rotate a stiffness tensor as follows:
+
+    >>> from elasticipy.tensors.elasticity import StiffnessTensor
+    >>> C = StiffnessTensor.cubic(C11=186, C12=134, C44=77) # Copper, mp-30
+    >>> C * goss
+    Stiffness tensor (in Voigt mapping):
+    [[ 1.86000000e+02  1.34000000e+02  1.34000000e+02  0.00000000e+00
+       0.00000000e+00  0.00000000e+00]
+     [ 1.34000000e+02  2.37000000e+02  8.30000000e+01  1.42108547e-14
+       0.00000000e+00  0.00000000e+00]
+     [ 1.34000000e+02  8.30000000e+01  2.37000000e+02 -7.10542736e-15
+       0.00000000e+00  0.00000000e+00]
+     [ 0.00000000e+00  0.00000000e+00  0.00000000e+00  2.60000000e+01
+       0.00000000e+00  0.00000000e+00]
+     [ 0.00000000e+00  0.00000000e+00  0.00000000e+00  0.00000000e+00
+       7.70000000e+01  0.00000000e+00]
+     [ 0.00000000e+00  0.00000000e+00  0.00000000e+00  0.00000000e+00
+       0.00000000e+00  7.70000000e+01]]
+
+Fibre textures
+==============
+Fibre textures are defined as uniformly distributed orientations around a given axis, leading to a single line when
+plotting an Orientation Distribution Function (ODF), hence the name. Therefore, there are three ways to define such
+textures.
+
+The first one consists in defining the *uvw* direction to align with a given axis (related to the sample coordinate
+system). For instance, let's consider the texture such that the <111> direction is aligned with the Z axis:
+
+    >>> from elasticipy.crystal_texture import FibreTexture
+    >>> from orix.crystal_map import Phase
+    >>> from orix.vector.miller import Miller
+    >>>
+    >>> phase = Phase(point_group='m-3m')   # Cubic symmetry
+    >>> m = Miller(uvw=[1,1,1], phase=phase)
+    >>> fibre_111 = FibreTexture.from_Miller_axis(m, [0,0,1])
+    >>> fibre_111
+    Fibre texture
+    <1. 1. 1.> || [0, 0, 1]
+
+Actually, this texture is usually referred to as the γ fibre; therefore, a quicker way to define it is to use:
+
+    >>> gamma = FibreTexture.gamma()
+
+This kind of shortcut also exists for alpha and epsilon textures
+(see `here <http://pajarito.materials.cmu.edu/lectures/L23-Metal_Textures-Torsion_BCC-part2.pdf>`_) for details).
+
+The third way is to use two (out of the three) Bunge-Euler angles to define the fix angles (assuming that the
+orientations are uniformly distributed over the remaining angle). E.g.:
+
+    >>> fibre_phi2 = FibreTexture.from_Euler(phi1=0., Phi=0.)
+    >>> fibre_phi2
+    Fibre texture
+    φ1= 0.0°, ϕ= 0.0°
+
+Then, the average stiffness resulting from such orientation distribution can be estimated from Voigt, Reuss or Hill
+method by integration over the all possible rotations. E.g.:
+
+    >>> Chill = C.Hill_average(orientations=fibre_phi2)
+
+One can check that this results in a transversely isotropic behaviour by plotting the values of the Young modulus in
+planar sections. The full code is:
+
+.. plot::
+
+    >>> from elasticipy.tensors.elasticity import StiffnessTensor
+    >>> from elasticipy.crystal_texture import FibreTexture
+    >>>
+    >>> C = StiffnessTensor.cubic(C11=186, C12=134, C44=77)
+    >>> fibre_phi2 = FibreTexture.from_Euler(phi1=0., Phi=0.)
+    >>> Chill = C.Hill_average(orientations=fibre_phi2)
+    >>> E = Chill.Young_modulus
+    >>> fig, ax = E.plot_xyz_sections()
+    >>> fig.show()
+
+Composite textures
+==================
+Different textures can be composed together to create a ``CompositeTexture`` object. For instance, if we consider a
+a cubic material which exhibits 40% of uniform texture, 30% of <111>||[0,0,1] and balanced copper texture, we can do the
+following:
+
+    >>> from elasticipy.crystal_texture import UniformTexture
+    >>> t = 0.4 * UniformTexture() + 0.3 * fibre_111 +  0.3 * DiscreteTexture.copper()
+    >>> t
+    Mixture of crystallographic textures
+     Wgt.  Type      Component
+     ------------------------------------------------------------
+     0.40  uniform   Uniform distribution over SO(3)
+     0.30  fibre     <1. 1. 1.> || [0, 0, 1]
+     0.30  discrete  φ1=90.00°, ϕ=35.26°, φ2=45.00°
+
+Again, the Hill average can be computed as follows:
+
+    >>> C.Hill_average(orientations=t)
+    Stiffness tensor (in Voigt mapping):
+    [[ 2.29927005e+02  1.15004605e+02  1.09068390e+02  8.38850529e-16
+       1.68295444e-01 -1.42501837e-15]
+     [ 1.15004605e+02  2.24466802e+02  1.14528593e+02 -2.58759922e-15
+       6.02903899e+00  9.94342227e-16]
+     [ 1.09068390e+02  1.14528593e+02  2.30403016e+02 -2.87690846e-16
+      -6.19733443e+00 -2.35261749e-15]
+     [ 9.65499494e-16 -2.35273493e-15 -1.02889997e-15  3.85692073e+01
+      -7.04774156e-17  5.34392791e+00]
+     [ 1.68295444e-01  6.02903899e+00 -6.19733443e+00 -5.84888929e-16
+       3.50396705e+01 -8.30517047e-16]
+     [ 5.75037319e-16  3.35378009e-15 -2.73847961e-16  5.34392791e+00
+       2.89814218e-16  3.85229349e+01]]
+
+Random sampling
+~~~~~~~~~~~~~~
+Sample of orientations can be drawn from all kind of textures. While a "random" sample from a single discrete texture
+does not make any sense, it can be of great interest for fibres or composite textures.
+
+For instance, a sample of 10 orientations can be drawn from the composite texture defined above as follows:
+
+    >>> sample = t.sample(num=10, seed=123) # Seed here is used to ensure reproducibility
+    >>> sample
+    Orientation (10,) 1
+    [[ 0.1157  0.5901  0.2194 -0.7683]
+     [ 0.3954 -0.6714 -0.421   0.4644]
+     [ 0.9571 -0.1986 -0.2022  0.0609]
+     [ 0.0802 -0.7012  0.3947 -0.5883]
+     [ 0.1952 -0.8854  0.3814  0.1802]
+     [-0.4814  0.4494  0.097   0.7463]
+     [ 0.3647 -0.2798 -0.1159 -0.8805]
+     [ 0.3647 -0.2798 -0.1159 -0.8805]
+     [ 0.3647 -0.2798 -0.1159 -0.8805]
+     [ 0.3647 -0.2798 -0.1159 -0.8805]]
+
+One may note that here, the 4 last orientations are the same. This is because it relates to the copper texture:
+
+    >>> DiscreteTexture.copper().orientation
+    Orientation (1,) 1
+    [[ 0.3647 -0.2798 -0.1159 -0.8805]]
+
+Plotting
+~~~~~~~~
+Pole figures can be drawn to evidence how each texture works. For instance, the (pure) Goss texture results in the
+following pole figure:
+
+.. plot::
+
+    >>> from elasticipy.crystal_texture import DiscreteTexture
+    >>> from orix.crystal_map import Phase
+    >>> from orix.vector.miller import Miller
+    >>>
+    >>> phase = Phase(point_group='m-3m')   # Cubic symmetry
+    >>> m = Miller(uvw=[1,0,0], phase=phase)
+    >>> Goss = DiscreteTexture.Goss()
+    >>> fig, ax = Goss.plot_as_pole_figure(m.symmetrise(unique=True))
+    >>> fig.show()
+
+Fibre textures can be drawn in a simular way. E.g.:
+
+.. plot::
+
+    >>> from elasticipy.crystal_texture import FibreTexture
+    >>> from orix.crystal_map import Phase
+    >>> from orix.vector.miller import Miller
+    >>>
+    >>> gamma = FibreTexture.gamma()
+    >>> point_group = 'm3m' # BCC symmetry
+    >>> m = Miller(uvw=[1,0,0], phase=Phase(point_group=point_group))
+    >>> fig, ax = gamma.plot_as_pole_figure(m.symmetrise(unique=True))
+    >>> fig.show()

{elasticipy-6.0.1 → elasticipy-6.1.0}/docs/source/Tutorials.rst

@@ -13,4 +13,5 @@ Tutorials
     ./Tutorials/Tutorial_MultiplePhases.rst
     ./Tutorials/Tutorial_ReadWriteFiles.rst
     ./Tutorials/Plasticity.rst
-    ./Tutorials/FourthOrderTensors.rst
+    ./Tutorials/FourthOrderTensors.rst
+    ./Tutorials/texture.rst

{elasticipy-6.0.1 → elasticipy-6.1.0}/docs/source/index.rst

@@ -3,7 +3,7 @@ Elasticipy Documentation
 
 .. |binder| image:: https://mybinder.org/badge_logo.svg
    :alt: Binder
-   :target: https://mybinder.org/v2/gh/DorianDepriester/Elasticipy/HEAD?urlpath=%2Fdoc%2Ftree%2FElasticipy_for_the_Impatient.ipynb
+   :target: https://mybinder.org/v2/gh/DorianDepriester/elasticipy-notebook/HEAD?urlpath=%2Fdoc%2Ftree%2FElasticipy.ipynb
 
 .. |joss| image:: https://joss.theoj.org/papers/10.21105/joss.07940/status.svg
   :alt: doi:10.21105/joss.07940
@@ -49,7 +49,7 @@ Elasticipy Documentation
 
 |pypi| |pypi-downloads| |conda| |conda-downloads|
 |github-commit| |github-stars| |license|
-|doc| |coverage| |pypi-version| |joss|
+|doc| |coverage| |pypi-version| |joss| |binder|
 
 
 Purpose of this package

{elasticipy-6.0.1 → elasticipy-6.1.0}/pyproject.toml

@@ -26,7 +26,8 @@ dependencies = [
     "matplotlib",
     "qtpy",
     "pyqt5",
-    "pandas"
+    "pandas",
+    "orix"
 ]
 requires-python = ">=3.9"
 
@@ -56,15 +57,8 @@ addopts = "--cov=src/Elasticipy --cov-report=xml"
 
 [tool.coverage.run]
 omit = [
-    "src/elasticipy/gui.py",
+    "src/elasticipy/gui/*",
 ]
 
-[tool.poetry]
-name = "jupyter-notebook-tutorial"
-version = "0.1.0"
-description = "Short tutorial for Elasticipy"
-authors = ["Dorian DEPRIESTER <dorian.dep@gmail.com>"]
-readme = "README.md"
-
 [project.gui-scripts]
-elasticipy-gui = "elasticipy.gui:crystal_elastic_plotter"
+elasticipy-gui = "elasticipy.gui:crystal_elastic_plotter"

elasticipy-6.1.0/requirements.txt

@@ -0,0 +1,18 @@
+numpy
+scipy
+matplotlib
+setuptools
+pandas
+qtpy
+PyQt5
+orix
+pytest
+pymatgen
+requests
+seaborn
+scikit-learn
+emmet-core
+ase
+matscipy
+paraqus
+sympy