Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Comsol 5.3 Meshfile not readable #2

Open
christoph33 opened this issue Dec 10, 2017 · 13 comments
Open

Comsol 5.3 Meshfile not readable #2

christoph33 opened this issue Dec 10, 2017 · 13 comments

Comments

@christoph33
Copy link

Hello
I tried a 2D Comsol mphtxt to convert into unv or msh files. Neither the -l listing is working.
I get errors like:
INFO: Reading mphtxt file ...
INFO: Reading piece 1 ...
INFO: Element type: Node
INFO: Element type: Edge lagrange P1
At line 117 of file source/mphtxt/module_manage_mphtxt.f90
Fortran runtime error: Index '0' of dimension 1 of array 'minelindx' below lower bound of 1

Error termination. Backtrace:

Program received signal SIGSEGV: Segmentation fault - invalid memory reference.

Backtrace for this error:

Or
Loading COMSOL mesh file...
INFO: Reading mphtxt file ...
INFO: Reading piece 1 ...
INFO: Element type: Edge lagrange P1
ERROR: mphtxt_file/object/etype, #-1
STOP 1

I tried with triangle, quad and mapped meshes and toggled the four options in export comsol:
Data elements
Boundaries elements
Vertex elements
Geometric entity information

Nothing works.
Any help?
feconv.exe compiled with gfortran in Cygwin.

Regards
Christoph

@fran-pena
Copy link
Collaborator

Dear Christoph,

it seems that the type of the finite element is not being properly recognized. In order to give you a more accurate answer I need the mesh file to debug the program because the tests I conducted with my COMSOL meshes were Ok.

Best regards,
Fran Pena

@christoph33
Copy link
Author

christoph33 commented Dec 11, 2017 via email

@christoph33
Copy link
Author

christoph33 commented Dec 12, 2017 via email

@fran-pena
Copy link
Collaborator

Dear Christoph,

Sorry I couldn't reply yesterday. I do not find your files neither in my GitHub account nor in the e-mail.

Best,
Fran

@christoph33
Copy link
Author

Dear Fran,
the missing files:
Transducer_1000.zip
Transducer_1000_tri.zip

@fran-pena
Copy link
Collaborator

Dear Christoph,

Now I have your meshes! Let me have some time to study the problem. I'll give you an answer.

Best,
Fran

@fran-pena
Copy link
Collaborator

Dear Christoph,

I think I found the reason why your mesh is not properly parsed. The Mesh type specifications are explained in the page 657 of the COMSOL Reference Guide ; see mesh-type.pdf.
In your mesh Transducter_1000_tri.mphtxt, the type mesh does not contain the parameter sets, domain labels and up/down boundaries.
Can you confirm me the COMSOL version you used to generate the mesh?

Best,
Fran

@fran-pena
Copy link
Collaborator

fran-pena commented Dec 15, 2017

Dear Christoph,
I have made a new commit that avoid that the MPHTXT files require to contain parameters or domains for Mesh type. In such a case, the element references does not exist and they are set to zero. I am not sure that to have missing parameters or domains is permitted by the MPHTXT standard. I am still interested into know how you created the mesh (COMSOL version, options used to save the mesh,...)

I have transformed Transducer_1000_tri.mphtxt and the result can be seen below: attached image.

Best,
Fran

@christoph33
Copy link
Author

christoph33 commented Dec 15, 2017 via email

@fran-pena
Copy link
Collaborator

Dear Christoph,
I cannot see the attached image.
Best,
Fran

@christoph33
Copy link
Author

christoph33 commented Dec 15, 2017 via email

@fran-pena
Copy link
Collaborator

Sorry, I still do not see any link to an image. Have you tried to click on "Attach files by dragging &, dropping, selecting them, or pasting from the clipboard"?
capture

@christoph33
Copy link
Author

christoph33 commented Dec 15, 2017 via email

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants