search

luchete80 / EPSol

0 stars 0 forks source link

More than 1 iteration shows dEdU singular matrix in disturbed flow #107

Closed luchete80 closed 3 years ago

luchete80 commented 3 years ago

This happens in 2 x 1 elements example

luchete80 commented 3 years ago

File "D:\Luciano\Numerico\epsol\src\Examples\Ex12\ex12_1_generic_uvFFvp.py", line 776, in dEdU=deriv.calc_dEdU(Fd,Fvpd,NsigF,NFvp) File "D:\Luciano\Numerico\epsol\src\Examples\Ex12\deriv.py", line 76, in calc_dEdU Fvpt_inv=linalg.inv(Fvpt) File "<__array_function__ internals>", line 5, in inv File "C:\Users\Luciano\AppData\Local\Programs\Python\Python39\lib\site-packages\numpy\linalg\linalg.py", line 546, in inv ainv = _umath_linalg.inv(a, signature=signature, extobj=extobj) File "C:\Users\Luciano\AppData\Local\Programs\Python\Python39\lib\site-packages\numpy\linalg\linalg.py", line 88, in _raise_linalgerror_singular

luchete80 commented 3 years ago

In one iterations shows the following: v,Uglob[v] 0 1.0 v,Uglob[v] 1 0.0 v,Uglob[v] 2 1.001 v,Uglob[v] 3 0.0 v,Uglob[v] 4 1.0 v,Uglob[v] 5 0.0 v,Uglob[v] 6 1.0 v,Uglob[v] 7 0.0 v,Uglob[v] 8 1.001 v,Uglob[v] 9 0.0 v,Uglob[v] 10 1.0 v,Uglob[v] 11 0.0 v,Uglob[v] 12 1.0 v,Uglob[v] 13 0.0 v,Uglob[v] 14 0.0 v,Uglob[v] 15 1.0 v,Uglob[v] 16 1.0004157816192465 v,Uglob[v] 17 0.0 v,Uglob[v] 18 -2.0789254205137395e-06 v,Uglob[v] 19 1.0 v,Uglob[v] 20 1.0002921836121481 v,Uglob[v] 21 0.0 v,Uglob[v] 22 -1.4609302350144738e-06 v,Uglob[v] 23 1.0 v,Uglob[v] 24 1.0 v,Uglob[v] 25 0.0 v,Uglob[v] 26 0.0 v,Uglob[v] 27 1.0 v,Uglob[v] 28 0.998570369763282 v,Uglob[v] 29 0.0 v,Uglob[v] 30 7.1482107516815635e-06 v,Uglob[v] 31 1.0 v,Uglob[v] 32 0.9992626900965899 v,Uglob[v] 33 0.0 v,Uglob[v] 34 3.686580238119408e-06 v,Uglob[v] 35 1.0

luchete80 commented 3 years ago

In second iteration Fvpt is:

Fvpt [[0. 0. 0.] [0. 0. 0.] [0. 0. 1.]]

Problem is with Fvpd