luchete80/epsol

epsol's Issues

Create Residuals from both forms

Create grid Info Member

Build Shape and Derivatives matrices (N and Bs) for all forms

Implement stability SUPG parameter tau

MPI ERROR

---

[2]PETSC ERROR:
[3]PETSC ERROR:
[0]PETSC ERROR:
[1]PETSC ERROR:
Caught signal number 11 SEGV: Segmentation Violation, probably memory access out of range
Caught signal number 11 SEGV: Segmentation Violation, probably memory access out of range
Caught signal number 11 SEGV: Segmentation Violation, probably memory access out of range
Caught signal number 11 SEGV: Segmentation Violation, probably memory access out of range

Write All DOFs outputs to vtk file

Cleaning all debug output

undefined reference to `MPI_XXXX'

text+0x231): undefined reference to MPI_Bcast' C:/msys64/mingw64/bin/../lib/gcc/x86_64-w64-mingw32/9.2.0/../../. text+0x2df): undefined reference to MPI_Bcast'
C:/msys64/mingw64/bin/../lib/gcc/x86_64-w64-mingw32/9.2.0/../../.
toseq.c:(.text+0x2e6): undefined reference to MPI_Comm_rank' C:/msys64/mingw64/bin/../lib/gcc/x86_64-w64-mingw32/9.2.0/../../. c:(.text+0x268): undefined reference to MPI_Win_fence'
C:/msys64/mingw64/bin/../lib/gcc/x86_64-w64-mingw32/9.2.0/../../.
c:(.text+0x2bd): undefined reference to MPI_Win_complete' C:/msys64/mingw64/bin/../lib/gcc/x86_64-w64-mingw32/9.2.0/../../. c:(.text+0x30e): undefined reference to MPI_Win_wait'
C:/msys64/mingw64/bin/../lib/gcc/x86_64-w64-mingw32/9.2.0/../../.
c:(.text+0x3b1): undefined reference to MPI_Win_free' C:/msys64/mingw64/bin/../lib/gcc/x86_64-w64-mingw32/9.2.0/../../. c:(.text+0x6c4): undefined reference to MPI_Win_fence'
C:/msys64/mingw64/bin/../lib/gcc/x86_64-w64-mingw32/9.2.0/../../.
c:(.text+0x965): undefined reference to MPI_Type_dup' C:/msys64/mingw64/bin/../lib/gcc/x86_64-w64-mingw32/9.2.0/../../. c:(.text+0xa8d): undefined reference to MPI_Type_create_indexed_
C:/msys64/mingw64/bin/../lib/gcc/x86_64-w64-mingw32/9.2.0/../../.
c:(.text+0xafd): undefined reference to MPI_Type_create_indexed_ C:/msys64/mingw64/bin/../lib/gcc/x86_64-w64-mingw32/9.2.0/../../. c:(.text+0xb5b): undefined reference to MPI_Type_commit'
C:/msys64/mingw64/bin/../lib/gcc/x86_64-w64-mingw32/9.2.0/../../.
c:(.text+0xbb9): undefined reference to MPI_Type_commit' C:/msys64/mingw64/bin/../lib/gcc/x86_64-w64-mingw32/9.2.0/../../. c:(.text+0xd3b): undefined reference to MPI_Type_get_extent'
C:/msys64/mingw64/bin/../lib/gcc/x86_64-w64-mingw32/9.2.0/../../.
c:(.text+0xd98): undefined reference to MPI_Type_get_true_extent C:/msys64/mingw64/bin/../lib/gcc/x86_64-w64-mingw32/9.2.0/../../. c:(.text+0xf69): undefined reference to MPI_Win_create'
C:/msys64/mingw64/bin/../lib/gcc/x86_64-w64-mingw32/9.2.0/../../.
c:(.text+0x100d): undefined reference to MPI_Win_fence' C:/msys64/mingw64/bin/../lib/gcc/x86_64-w64-mingw32/9.2.0/../../. c:(.text+0x10d6): undefined reference to MPI_Win_post'
C:/msys64/mingw64/bin/../lib/gcc/x86_64-w64-mingw32/9.2.0/../../.
c:(.text+0x113d): undefined reference to MPI_Win_start' C:/msys64/mingw64/bin/../lib/gcc/x86_64-w64-mingw32/9.2.0/../../. c:(.text+0x138f): undefined reference to MPI_Win_lock'
C:/msys64/mingw64/bin/../lib/gcc/x86_64-w64-mingw32/9.2.0/../../.
c:(.text+0x1425): undefined reference to MPI_Accumulate' C:/msys64/mingw64/bin/../lib/gcc/x86_64-w64-mingw32/9.2.0/../../. c:(.text+0x1482): undefined reference to MPI_Win_unlock'
C:/msys64/mingw64/bin/../lib/gcc/x86_64-w64-mingw32/9.2.0/../../.
c:(.text+0x1688): undefined reference to MPI_Win_lock' C:/msys64/mingw64/bin/../lib/gcc/x86_64-w64-mingw32/9.2.0/../../. c:(.text+0x1719): undefined reference to MPI_Get'
C:/msys64/mingw64/bin/../lib/gcc/x86_64-w64-mingw32/9.2.0/../../.
c:(.text+0x1776): undefined reference to MPI_Win_unlock' C:/msys64/mingw64/bin/../lib/gcc/x86_64-w64-mingw32/9.2.0/../../. c:(.text+0x1933): undefined reference to MPI_Type_free'
C:/msys64/mingw64/bin/../lib/gcc/x86_64-w64-mingw32/9.2.0/../../.
c:(.text+0x19a3): undefined reference to MPI_Type_free' C:/msys64/mingw64/bin/../lib/gcc/x86_64-w64-mingw32/9.2.0/../../. c:(.text+0x19f9): undefined reference to MPI_Type_free'
C:/msys64/mingw64/bin/../lib/gcc/x86_64-w64-mingw32/9.2.0/../../.
c:(.text+0x1bb7): undefined reference to MPI_Win_free' C:/msys64/mingw64/bin/../lib/gcc/x86_64-w64-mingw32/9.2.0/../../. c:(.text+0x2141): undefined reference to MPI_Win_lock'
C:/msys64/mingw64/bin/../lib/gcc/x86_64-w64-mingw32/9.2.0/../../.

Single Element viscoplastic example 12 (c++ & py)

Calculate Rate of almansi strain derivative (dE/dU)

Calculate Fe form 2

Correct example 2 to 100 elements and petsc

FValues of Grid Element Example 2 crash

Creating derivatives
gi0
setting global grad
set
Creating Matrices
Creating Elemental Stiffness Matrix, GaussPoints: 4
B matrix cols8
rows3
nan 0 nan 0 nan 0 nan 0
0 nan 0 nan 0 nan 0 nan
nan nan nan nan nan nan nan nan

cols8
rows8
nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan
B matrix terminate called after throwing an instance of 'std::bad_alloc'
what(): std::bad_alloc

Include several element plane strain example

Instantiate output

Order 3 integration fails

Option to build DOFHandler without METIS and CSR

Create cilindrical mesh for radial flow example

Create matrix operator(i,j) for i/o

Instantiate FEValues.h

Set up Newton Raphson Mehod in python and c++(Petsc)

Coorect vtk output from cell to element and fields

• Write a single solution
• Store output fields (Integration points and nodes)

Resolve "AtGaussPoint" for Rv

Gauss Full Matrices outstr() not working

Crash in FEValues in runtime

Program stuck in derivative matrices creation

Create definitions for including PETSC or LASPACK

https://stackoverflow.com/questions/9017573/define-preprocessor-macro-through-cmake

Reemplace matrices indices integers per variables

Creating dofxnod, dofxelem
eg:
16--> nodxel*dof+dofxnod

Example 14: die

Build elastic strain rate derivatives (Ee)

For both formulations

DOFHandler CSR ERROR

D:\Luciano\Numerico\epsol\bin_examples\ex2>ex2
pushing back elems
pushing back elems
Total Grid Nodes6
Initializing DoF Handler...
Num Elements: 2
Element field degree:
4259376Field Dim:
Mesh degree: 4259376
vertex dof size6
field dim 2DoF Number: 12
adding adj dofs number
Element number 2
adding adj dofs number
Creating csr
Loop through nodes...
csr xadj 0 = 0
csr xadj 1 = 2
csr xadj 2 = 5
csr xadj 3 = 7
csr xadj 4 = 9
csr xadj 5 = 12
csr xadj 6 = 14

Subdivide shared libraries

Assemble element tangent matrix

Example 5 crashes

Add vtk output to examples 1 a/o 2

Measure time in example 2

Working Example 2

Still not working, carshing at solving

Calculate gauss points calculated variables (v,F,Fvp, sig,s)

Disturbed Flow Model to compare 2 formulations

From Qin chapter 4 example.

Prevent to create incompatible inverted Jacobian

Calculate viscoplastic strain rate tensor dD/dUx DERIVATIVES matrices for Formulation 2

Create and check direction of Plastic Flow Tensor Nij

Example 13:radial flow

• Apply bcs

• Compare results

Calculate viscosity in form 2

" no match for operator[]" ^"

D:/Luciano/Numerico/epsol/src/FiniteElement/FEValues/FEValues.h:112:123: error: no match for 'operator[]' (operand types are 'const FluxSol::GaussFullMatrices' and 'const int') 112 | inline const double & shape_value_component (const int &fn, const int &gaussp) const{return this->shape_value_components[gaussp][0][fn];}; //Like in dealii |

Calculate dE/dU for formulation 2

Create elemental tangent Matrix

Create incremental stiffness f(s,sig) for material

Possibility of choose solver tolerances from constructors

Calculate Second Piola-Kirchoff Tensor

Set element charact. length for SUPG