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.
fem.refel
fem.mesh
fem.xform