Hari Sundar


So I have been passively following the development of Julia for a while now. It sounds very promising and I hope that I can finally replace matlab as my prototyping environment. In a series of blogs, I am going to port my matlab code for high-order multigrid to Julia and compare its performance. I shall also consider the effort to parallelize the code in Julia. I plan to make small modifications to the code that should not affect the performance.

In the first part of this series, we shall build a simple hexahedral meshing framework. The matlab version used a simple regular grid like mesh. In this version, we shall modify this to support unstructured quad and hex meshes. This is mainly to eventually support adaptive quadtree/octree meshes. We shall start by writing a class for the reference element.

Reference Element — fem.refel

Mesh — fem.mesh

Transform — fem.xform