On 05/25/10 23:33, Logan Sorenson wrote:
Hey Robert,
On Mon, May 24, 2010 at 4:34 AM, Robert Cimrman<cimr...@ntc.zcu.cz> wrote:
Hi Logan,
On 05/22/10 18:59, Logan Sorenson wrote:
Hi all,
I added to the existing Comsol meshio faciliites to implement full read/write for SfePy suported types. Hope someone can find it useful. At the very least, it should make comparisons between Comsol and SfePy easier. The code is on the comsol_mesh branch at [1].
I have tried several meshes with different kinds of elements, e.g.
./script/convert_mesh.py meshes/3d/cylinder.mesh cyli.txt ./postproc.py cyli.txt
and it worked perfectly!
I tried to run the mesh generated by Salome and converted using convert_mesh.py to Comsol format. It looks like Comsol can read the mesh okay, but when I run the Comsol solution it gives some errors about inverted elements and the solution looks strange compared to a mesh generated directly in Comsol. Therefore, I think that something is not quite right in the implementation of the Comsol write code. I know from reading and writing a Comsol generated mesh file that some of the elements get reordered by SfePy, so maybe this is causing the error. I will look more into it, but just wanted to post a warning if anyone is planning to use this to create Comsol meshes.
FYI: The orientation of elements (= ordering of vertices in connectivity) is depicted in [2].
cheers, r.
[1] http://github.com/logansorenson/sfepy/tree/comsol_mesh [2] http://docs.sfepy.org/doc-devel/users_guide.html#fields