CHOLMOD warning and Segmentation fault
Hi,
I use the PeriodicFlowEngine to simulate the erosion, I get such warning:
CHOLMOD warning: matrix not positive definite. file: ../Supernodal/
something went wrong in Cholesky factorization, use LDLt as fallback this time1
Segmentation fault (core dumped)
I am looking forward to your help ? Thanks a lot !
My code is shown as follow:
##______________ First section, generate sample_________
from __future__ import print_function
from yade import pack, qt, plot
from math import *
nRead=readParam
## model parameters
## material parameters
young=2e8,
poisson=.2,
alphaKtw=0,
competaRoll=.22,
etaTwist=0,
## fluid parameters
## control parameters
damp=0,
## output specifications
)
from yade.params.table import *
O.periodic=True
O.cell.
# create materials for spheres
#shear strength is the sum of friction and adhesion, so the momentRotationL
O.materials.
# generate particles packing
sp=pack.
sp.makeCloud(
sp.toSimulation
O.engines=[
),
# specify target values and whether they are strains or stresses
# type of servo-control, the strain rate isn't determined, it shloud check the unbalanced force
# wait until the unbalanced force goes below this value
),
]
import sys
def compactionFinis
# after sample preparation, save the state
O.save(
print(
# next time, called python command
triax.
O.pause()
O.run()
O.wait()
#B. Activate flow engine
flow.dead=0
flow.defToleran
flow.meshUpdate
flow.useSolver=3
flow.permeabili
flow.viscosity=10
flow.gradP=
flow.duplicateT
flow.updateTria
O.run(2,1)
csdList=
print(len(
flow.dead=1
print(O.
for i in range(10,80):
O.bodies.
print(len(
O.run(1000,True)
O.engines=
O.step()
O.step()
O.engines=
flow.defToleran
flow.meshUpdate
flow.useSolver=3
flow.permeabili
flow.viscosity=10
flow.gradP=
flow.duplicateT
flow.updateTria
O.run(2,1)
csdList1=
print(len(
Question information
- Language:
- English Edit question
- Status:
- Expired
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Last query:
- Last reply: