Can not read elements when using ComponentElmer

Dear Experts

I calculated an electric field map using Gmsh and ElmerSolver, and want to import it into Garfield++. I got an error when reading the elements but I don’t know what does it mean.
My code:

fCompE_field = new ComponentElmer(
            "fieldmap/prototype/mesh.header", "fieldmap/prototype/mesh.elements", "fieldmap/prototype/mesh.nodes",
            "fieldmap/prototype/dielectrics.dat", "fieldmap/prototype/prototype.result", "mm");

And I got such output:

ComponentElmer::Initialise:
    Read 568725 nodes and 3529037 elements from file fieldmap/prototype/mesh.header.
ComponentElmer::Initialise:
    Set material 0 of 10 to eps 1e+10.
ComponentElmer::Initialise:
    Set material 1 of 10 to eps 1e+10.
ComponentElmer::Initialise:
    Set material 2 of 10 to eps 1e+10.
ComponentElmer::Initialise:
    Set material 3 of 10 to eps 1e+10.
ComponentElmer::Initialise:
    Set material 4 of 10 to eps 1e+10.
ComponentElmer::Initialise:
    Set material 5 of 10 to eps 1e+10.
ComponentElmer::Initialise:
    Set material 6 of 10 to eps 1e+10.
ComponentElmer::Initialise:
    Set material 7 of 10 to eps 1e+10.
ComponentElmer::Initialise:
    Set material 8 of 10 to eps 1e+10.
ComponentElmer::Initialise:
    Set material 9 of 10 to eps 1.
ComponentElmer::Initialise:
    Error reading file fieldmap/prototype/mesh.elements (line 0).

Can you provide the field map files you are using?

@hschindl It’s just too large…(100MB after compress) so I put on google drive.

Your mesh.elements file is for a linear tetrahedral mesh (four nodes per element). ComponentElmer expects a second-order tetrahedral mesh (ten nodes per element).

This topic was automatically closed 14 days after the last reply. New replies are no longer allowed.