tdegeus / goosefem Goto Github PK

Introduce nnodePeriodic for all meshes, sync with Mesh::Quad4::Regular

• Mesh::Quad4::Regular
• Mesh::Quad4::FineLayer
• Mesh::Quad8::Regular
• Mesh::Quad8::FineLayer
• Mesh::Tri3::Regular

Apply code style everywhere.

Checkout automated tool used in HighFive

• overdamped/Verlet
• overdamped/velocityVerlet
• eta=00.10__alpha=0.00
• eta=00.10__alpha=0.20
• overdamped

It may be wise to create a more transparent setup.py in which xtensor, xtl, xsimd are automatically downloaded (as is pyxtensor).

Simplify implementation using new T2s::Copy( ... ) with one argument.

Benefit:

• Forces good setup.py
• Allows extension to other OS.

• introduction.rst
• examples/statics_linear.rst
• examples/dynamics_inertia.rst
• examples/dynamics_inertia-damping.rst
• examples/dynamics_overdamped.rst
• conventions/terminology.rst
• conventions/storage.rst
• conventions/vector.rst
• conventions/matrix.rst
• conventions/allocation.rst
• details/macros.rst
• details/Mesh.rst
• details/MeshTri3.rst
• details/MeshQuad4.rst
• details/MeshHex8.rst
• details/Element.rst
• details/ElementQuad4.rst
• details/ElementHex8.rst
• details/Vector.rst
• details/Matrix.rst
• details/Tyings.rst
• details/ParaView.rst
• details/Iterate.rst
• compile.rst
• develop.rst

Benefit:

• Keeps examples up-to-date.
• Forces good CMakeLists.txt
• Makes sure that examples stay up-to-date with material model.
• Allows extension to other OS.

Use the new rank members from xtensor to pave the way for removing the copies in the Python interface

Implement:

• asDense
• asSparse
• set

(and add to python.cpp)

• docs/examples/statics/FixedDisplacements_LinearElastic
• docs/examples/statics/MixedPeriodic_LinearElastic
• docs/examples/statics_linear.rst

Change all:

and -> &&
or -> ||

E.g.

using T2 = xt::xtensor_fixed<double, xt::xshape<2,2>>;

In ElementQuad4.h

• MatrixDiagonal.h
• Vector.h

Rename m_tyings to m_nodaltyings

Also (temporarily) deprecate all support functions.

I'm getting

Traceback (most recent call last):
  File "convert.py", line 4, in <module>
    import GooseFEM as gf
ImportError: dlopen(/usr/local/lib/python3.7/site-packages/GooseFEM-0.2.3-py3.7-macosx-10.15-x86_64.egg/GooseFEM.cpython-37m-darwin.so, 2): Symbol not found: __ZNK8GooseFEM23MatrixPartitionedTyingsIN5Eigen14SimplicialLDLTINS1_12SparseMatrixIdLi0EiEELi1ENS1_11AMDOrderingIiEEEEE3nndEv
  Referenced from: /usr/local/lib/python3.7/site-packages/GooseFEM-0.2.3-py3.7-macosx-10.15-x86_64.egg/GooseFEM.cpython-37m-darwin.so
  Expected in: flat namespace
 in /usr/local/lib/python3.7/site-packages/GooseFEM-0.2.3-py3.7-macosx-10.15-x86_64.egg/GooseFEM.cpython-37m-darwin.so

Increase readability of CMakelists.txt, see https://github.com/tdegeus/ElastoPlasticQPot/blob/master/src/CMakeLists.txt

Notably all "static" example. Note that it would be good to add non-linear history dependent examples with fixed displacements.

Should point to https://goosefem.readthedocs.io

This way:

ce[m,i] = Ae[m,i,j,n] * b[n,j]

whereby Ae comes from

np.einsum('mi,ijkl,nl->mjnk', dNdx, A4, dNdx)

Note that the current definition is just

np.einsum('mi,ijkl,nl->mjnk', dNdx, A4, dNdx).reshape(nne*ndim, nne*ndim)

• ElementQuad4.h

Examples:

• https://github.com/conda-forge/eigen-feedstock
• https://gomss-nowcast.readthedocs.io/en/latest/conda_pkgs.html
• https://gist.github.com/zshaheen/fe76d1507839ed6fbfbccef6b9c13ed9

cd /path/to/GMatElastoPlasticFiniteStrainSimo
cd conda-recipes
conda config --append channels conda-forge
conda build --croot /tmp/cbtmp .

See xtensor-stack/xtensor#1901

See xtensor-stack/xtensor#1769

When default arguments were introduced to program gave a segmentation fault.