#!/usr/bin/env python3
# -*- coding: utf-8 -*-

# Gebiet Omega für Aufgabe 25

from netgen.geom2d import * 

Omega = SplineGeometry()
pnts = [(0,0), (2,0), (2,1), (1,1), (1,2), (0,2)]
pnums = [Omega.AppendPoint(*p) for p in pnts]

Omega.Append(['line', pnums[0], pnums[1]], bc='Gamma1')
Omega.Append(['line', pnums[1], pnums[2]], bc='Gamma2')
Omega.Append(['spline3', pnums[2], pnums[3], pnums[4]], bc='Gamma1')
Omega.Append(['line', pnums[4], pnums[5]], bc='Gamma2')
Omega.Append(['line', pnums[5], pnums[0]], bc='Gamma1')

if __name__== '__main__':
    import netgen.gui
    from ngsolve import *
    mesh = Mesh(Omega.GenerateMesh(maxh=0.2))
    Draw(mesh)