smesh/test/netgen_runner_2D.py

#!/usr/bin/env python
import sys
import salome

salome.salome_init()

from os import path
import tempfile
import subprocess

import GEOM, SMESH, SALOMEDS

from salome.geom import geomBuilder
from salome.smesh import smeshBuilder

import math

smesh = smeshBuilder.New()
geompy = geomBuilder.New()

import medcoupling as mc

def create_param_file(param_file):
    """ Create a parameter file for runner """
    param = """1
34.641
0.14
16
0.15
1.5
0
0
1
5
1
1
-1
3
3
0.2
2
1
0
0
2
2
0

0
0
0"""
    with open(param_file, "w") as ffile:
        ffile.write(param)

def create_maxarea_param_file_2D(param_file,maxArea):
    """ Create a parameter file for runner """
    param = """1
{}
""".format(maxArea)
    with open(param_file, "w") as ffile:
        ffile.write(param)

def create_lenghtFromEdges_param_file_2D(param_file):
    """ Create a parameter file for runner """
    param = """1
"""

    with open(param_file, "w") as ffile:
        ffile.write(param)

def CommonFunction(netgen,case,numberOfEdges1D,hypo,maxArea):
    # Building geometry
    geometry = geompy.MakeBoxDXDYDZ(200, 200, 200)
    # Create 1D regular elements
    mesh_1d = smesh.Mesh(geometry, 'Maillage_1')
    mesh_1d.Segment().NumberOfSegments(numberOfEdges1D)
    isDone = mesh_1d.Compute()

    if not isDone:
        raise Exception("Error when computing Mesh")

    smesh.SetName(mesh_1d, 'Maillage_1')
    nb_points = 0
    nb_segments = 0
    nb_triangles = 0
    nb_tetras = 0
    with tempfile.TemporaryDirectory() as tmp_dir:
        mesh_file = path.join(tmp_dir, "mesh.med")
        shape_file = path.join(tmp_dir, "shape.brep")

        if ( hypo == 0 ): 
            param_file = path.join(tmp_dir, "param.txt")
            create_param_file(param_file )

        if ( hypo == 1 ):
            param_file = path.join(tmp_dir, "param_lenghtfromedge.txt")
            create_lenghtFromEdges_param_file_2D(param_file ) 

        if ( hypo == 2 ):
            param_file = path.join(tmp_dir, "param_maxarea.txt")
            create_maxarea_param_file_2D(param_file,maxArea)            
        
        bin_file    =  path.join(tmp_dir, "mesh.bin")
        output_mesh = path.join(tmp_dir, "meshout.med")

        print("Running in folder: ", tmp_dir)
                
        mesh_1d.ExportMED(mesh_file, 0, 41, 1, mesh_1d, 1, [], '', -1, 1)
        geompy.ExportBREP( geometry, shape_file )

        runner = path.join("${NETGENPLUGIN_ROOT_DIR}",
                            "bin",
                            "salome",
                            "NETGENPlugin_Runner")

        if sys.platform == 'win32':
            runner += ".exe"

        cmd = "{runner} {NETGEN} {mesh_file} {shape_file} "\
                "{param_file} NONE {bin_file} {output_mesh}"\
                .format(runner=runner,
                        NETGEN=netgen,
                        mesh_file=mesh_file,
                        shape_file=shape_file,
                        param_file=param_file,
                        bin_file=bin_file,
                        output_mesh=output_mesh)
        print(cmd)
        subprocess.check_call(cmd, shell=True)   
        mesh_read = mc.ReadUMeshFromFile(output_mesh, "MESH", 0)
        nb_points = mesh_read.getNumberOfNodes()
        if (case == 2):
            nb_triangles = mesh_read.getNumberOfCellsWithType(mc.NORM_TRI3)
            mesh_read    = mc.ReadUMeshFromFile(output_mesh, "MESH", -1)
            nb_segments  = mesh_read.getNumberOfCellsWithType(mc.NORM_SEG2)
            
    return [nb_points,nb_segments,nb_triangles,nb_tetras] 
    

def test_netgen2d():
    """ Test netgen2d mesher """
    [nb_points,nb_segments,nb_triangles,nb_tetras] = CommonFunction( "NETGEN2D", 2, 3, 0, 0.0 )
    
    print("Nb Triangles:", nb_triangles)
    print("Nb Segments:", nb_segments)
    print("Nb Points:", nb_points)

    assert nb_triangles > 12
    assert nb_points > 8
    assert nb_segments > 12

def test_netgen2dMaxArea():
    """ Test netgen2d mesher """
    [nb_points,nb_segments,nb_triangles,nb_tetras] = CommonFunction( "NETGEN2D", 2, 3, 2, 75.0 )
    
    print("Nb Triangles:", nb_triangles)
    print("Nb Segments:", nb_segments)
    print("Nb Points:", nb_points)

    assert nb_triangles > 12
    assert nb_points > 8
    assert nb_segments > 12

def test_netgen2dLenghtFromEdge():
    """ Test netgen2d mesher """
    [nb_points,nb_segments,nb_triangles,nb_tetras] = CommonFunction( "NETGEN2D", 2, 1, 1, 0.0 )
    
    print("Nb Triangles:", nb_triangles)
    print("Nb Segments:", nb_segments)
    print("Nb Points:", nb_points)

    assert nb_triangles == 12*2
    assert nb_points == 14
    assert nb_segments == 12

if __name__ == "__main__":
    if sys.platform == "win32":
        print("Disabled on windows")
        sys.exit(0)
    test_netgen2d()
    test_netgen2dMaxArea()
    test_netgen2dLenghtFromEdge()