I am trying to initial stress in a 2D subduction zone with a triangle mesh. I put the gravity (using the example) in my subduction (same example than in the manual) and now I am trying to do the step8. I already re-run the step3 with commenting out the CellFilterAvg settings and run pylith gravity_initstress.cfg gravity_vardensity.cfg. But when I did ./generate_statedb.py, I have this error:
“Traceback (most recent call last):
File “./generate_statedb.py”, line 63, in
quadCoords = calcQuadCoords(vertices, cells, qpts)
File “./generate_statedb.py”, line 42, in calcQuadCoords
quadCoords[:,:,iDim] = numpy.dot(cellCoords[:,:,iDim], qpts.transpose())
ValueError: shapes (523,3) and (4,4) not aligned: 3 (dim 1) != 4 (dim 0)”
I understood that in the generate_statedb.py file the qpts array are the basis function evaluated at the quadrature points but I still don’t understand how to calculate this value? Furthermore, the example of 2D gravity is for lines and my model is for triangle so I need to say in generate_statedb.py the good size. But I don’t understand how?
Can you help me?
Thank you in advance,
Phd Student Institut des Sciences de la Terre de Paris