search

luchete80 / EPSol

0 stars 0 forks source link

FValues of Grid Element Example 2 crash #20

Closed luchete80 closed 4 years ago

luchete80 commented 4 years ago

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

luchete80 commented 4 years ago

even overridying CSR and METIS is another error

luchete80 commented 4 years ago

The problem is not assemblying

Element: 0 GLOBAL DOFS8 GLOBAL DOFS9 GLOBAL DOFS6 GLOBAL DOFS7 GLOBAL DOFS0 GLOBAL DOFS1 GLOBAL DOFS2 GLOBAL DOFS3 Element: 1 GLOBAL DOFS10 GLOBAL DOFS11 GLOBAL DOFS8 GLOBAL DOFS9 GLOBAL DOFS2 GLOBAL DOFS3 GLOBAL DOFS4 GLOBAL DOFS5

luchete80 commented 4 years ago

The problem is in : FluxSol::FEValues<2> fev(grid.Elem(e),grid); Assigning a test element works ok

luchete80 commented 4 years ago

Output is as following:

setting global grad set Creating Matrices J matrix cols2

rows4 0 2 0.5 1 0 0 0.5 0

luchete80 commented 4 years ago

[0,0,0] [0.5,0,0] [0,1,0] [0.5,1,0] [0,2,0] [0.5,2,0]