#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
Created on Tue May 21 15:09:52 2019

@author: mfischer
"""

from netgen.meshing import Element2D, MeshPoint, FaceDescriptor, Mesh
from netgen.csg import Pnt

# Gitter
meshRef = Mesh()
meshRef.dim = 2

# Ecken zu mesh hinzufügen
pnumsMesh = []
pnumsMesh.append(meshRef.Add(MeshPoint(Pnt(0,0,0))))
pnumsMesh.append(meshRef.Add(MeshPoint(Pnt(1,0,0))))
pnumsMesh.append(meshRef.Add(MeshPoint(Pnt(0,1,0))))

# Flächen (also 2D-Elemente) hinzufügen
meshRef.Add(FaceDescriptor(surfnr=1, bc=1))
meshRef.Add(Element2D(1, [pnumsMesh[0], pnumsMesh[1], pnumsMesh[2]]))

from ngsolve import Mesh
mesh = Mesh(meshRef)

if __name__== '__main__':
    import netgen.gui
    from ngsolve import Draw
    Draw(mesh)