Skip to content
  • Mark OLESEN's avatar
    ENH: robuster handling of inv() of fields with singular tensors (#2724) · a96dcf70
    Mark OLESEN authored
    - replace the "one-size-fits-all" approach of tensor field inv()
      with individual 'failsafe' inverts.
    
      The inv() field function historically just checked the first entry
      to detect 2D cases and adjusted/readjusted *all* tensors accordingly
      (to avoid singularity tensors and/or noisy inversions).
    
      This seems to have worked reasonably well with 3D volume meshes, but
      breaks down for 2D area meshes, which can be axis-aligned
      differently on different sections of the mesh.
    a96dcf70