DG Method From Scratch (in Python) - Mapping
From reference space to physical space

last update: Feb 3, 2025

Article Motivation

  • How do I map my nodes from the reference element to the mesh?
  • How do I find the Jacobians and surface Jacobians for each cell and cell face in my mesh?
  • How do I find the normal vectors for each face in each cell?
  • Code it in python

Introduction

In the last post on extracting information from a Gmsh mesh we created numpy arrays that held the necessary connectivity data about our mesh, telling us things like which cells were adjacent and which faces were touching.

Now we need to go one step further and map our nodal basis from the reference tetrahedron to each cell in our mesh. This means we need to create numpy arrays for the coordinates for each node in each cell in our mesh.

References