CHOLMOD warning: matrix not positive definite with PeriodicFlowEngine
Dear All,
I have tried to run the following code (simple shear test upon a saturated granular sample under constant pressure conditions):
from yade import pack,qt,
from yade import export
from yade import timing
from yade import plot
from builtins import range
import numpy
import time
from math import *
O.load('prova.xml')
num_spheres=5000
young=750.e6
finalFricDegree = 0.35
graindensity=2600.0
poisson=0.17
en=0.9
sigmaIso=-150.e3
radiusR=0.001/2.
D50=radiusR*2.0
#O.materials.
triax=PeriTriax
dynCell=True,
mass=0.2,
goal=[
stressMask=7,
maxStrainRate=
label='triax')
newton=
O.engines=[
ForceResetter(),
InsertionSortC
InteractionLoop(
),
PeriodicFlowEn
GlobalStiffnes
triax,
newton
]
while 1:
O.run(1000, True)
unb=unbalanced
if unb<0.00001:
break
triax.dead=0
flow.dead=0
flow.pressureFo
flow.viscousShe
flow.viscousNor
flow.viscousShe
flow.normalLubr
flow.shearLubri
flow.defToleran
flow.meshUpdate
flow.useSolver=3
flow.permeabili
flow.viscosity=
flow.bndCondIsP
flow.bndCondVal
GlobalStiffness
O.dt=min(
O.dynDt=False
O.run(10,1)
while 1:
O.run(5000, True)
unb=unbalanced
if unb<0.0001:
if O.cell.
O.cell.
O.dt=
if O.cell.
break
with "prova.xml" a dry granular sample properly generated in a periodic cell.
Even if I switch off both the PeriTriaxialCon
CHOLMOD warning: matrix not positive definite. file: ../Supernodal/
something went wrong in Cholesky factorization, use LDLt as fallback this time1
Would you mind helping me with this error?
Question information
- Language:
- English Edit question
- Status:
- Solved
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Solved by:
- Katia Boschi
- Solved:
- Last query:
- Last reply: