Add basis method for finding el num from nb side num and mon num, rename similar existing fe/face -> el num method for clarity.
Replace 'non_boundary' with 'nb' in fn name for consistency.
Fix critical boundary side determination inversion of logic, and enable tests that should have caught the problem but weren't attributed as tests.
Remove unnecessary casts related to unwrapped sidefaces.
Add num_els method and change interacting els estimate to count interacting pairs instead of triplets, which is now an exact count in most cases.
Unit tests for new interacting pairs estimation function of last commit.
Add parameter to interacting elements estimation function to only compute interacting elements in upper triangle.
Rename field of NBSideInclusions structure for consistency.
Drop unneeded second size dimension from upper triangle constructor, simplify impl, add lower triangle version.
Clarify function doc.
Add env variable so tests run sequentially.
Add R_NaN.
Add tensor type for multi store and fetch of reals.
Add sparse matrix type.
Add trait impls to allow range() on newtype wrapped numeric types.
Add variational bilinear form type.
Revert One/Add impls for newtype wrapped integral types, because range() used directly on these types is slower than for built-in integral types.
Remove basis parameter on methods, provide basis accessor instead.
Add accessor for weak_grad_ops and rename wgrad accessors.
Changed to use new format! for conversion to string.