search

ju-kreber / paraview-scripts

Collection of python scripts for interfacing with paraview for the Chair for Electromagnetic Theory at Saarland University
10 stars 2 forks source link

Wedge interior nodes are outside of element domain #1

Open sandervdbroek opened 1 year ago

sandervdbroek commented 1 year ago

I noticed that interior nodes of wedges seem to be generated outside of the domain of the element.

If I specify the vertices to be [-1,1,-1,-1,1,-1],[-1,-1,1,-1,-1,1],[-1,-1,-1,1,1,1]] it generates the following nodes in space: wedgenodes With [-1.33333333 -1.33333333 -1.33333333] [-1.33333333 -1.33333333 -0.66666667]] not being within the domain of the wedge. This can be manually adjusted, but should probably be fixed in the code

sandervdbroek commented 1 year ago

I believe it can be fixed by changing the pos_z to be zero initially - but I am unsure if this is a general solution

'# pos_z = corner_verts[0].copy()' '# pos_z = np.expand_dims(pos_z, axis=0)' 'pos_z = np.zeros([1,3])`

sandervdbroek commented 1 year ago

I believe one of the wedge faces also has the order of vertices different from what paraview uses

# quadrilateral_faces = [(0,1,4,3), (1,2,5,4), (0,2,5,3)]
quadrilateral_faces = [(0, 1, 4, 3), (1, 2, 5, 4), (2, 0, 3, 5)]