'ip_simple' mimetic inner product. Coincides with the results from
'computeMimeticIP' to (roughly) the order of round-off on a
non-trivial test case.
The calling interface is
BI = mex_ip_simple(G, rock)
and some latitude has been extended towards multiple data types in
the various G fields (e.g., G.cells.facePos may be an 'int32').
Additional clean-up and optimisation is likely.
- mim_ip_span_nullspace()
Constructs orthogonal basis for mimetic IP null space.
- mim_ip_linpress_exact()
Adds term to impose exactness for linear pressure fields.
Re-implement mim_ip_simple() in terms of these procedures.
Current status
--------------
Given vectors ZCORN, COORD and ACTNUM as well as the Cartesian
dimensions these vectors implicitly refer to, the code is
currently capable of
* Identify unique points along each pillar
* Assign point numbers for each point specified in ZCORN
* Compute face topology, i.e., the corners that define the geometry
of the faces as well as the cells that are connected through the face.
* Identify and compute intesections that occur in the processing of
face topology.
What remains is
* Handle the face geometry of boundary faces. (simple)
* Compute point coordinates of the final point list.
* Put all pieces together in a tidy manner.
