anchor.diff:                                  det <- det * x'.[i, i]
                           Some(det)
                       | _ -> None
  -            member o.ReshapeCopy_MRows_V(x) = // Non-BLAS
  +            // Non-BLAS
  +            member o.ReshapeCopy_MRows_V(x) =
                   if Array2D.isEmpty x then
                       Array.empty<float>
                   else
                       let r = Array.zeroCreate<float> x.Length
                       Buffer.BlockCopy(x, 0, r, 0, x.Length * sizeof<float>)
                       r
  -            member o.ReshapeCopy_V_MRows(m, x) = // Non-BLAS
  +            // Non-BLAS
  +            member o.ReshapeCopy_V_MRows(m, x) =
                   if Array.isEmpty x then
                       Array2D.empty<float>
                   else
  
positive.code:                     let x' = Array2D.copyFast x
                      BLAS.dscal'(alpha, x')
                      x'
              member o.Sub_M_M(x, y) =
                  if Array2D.isEmpty x then
                      (o :> Backend<float>).Mul_S_M(-1., y)
  
negative.code:                                 det <- det * x'.[i, i]
                          Some(det)
                      | _ -> None
              // Non-BLAS
              member o.ReshapeCopy_MRows_V(x) =
                  if Array2D.isEmpty x then
                      Array.empty<float32>
                  else
                      let r = Array.zeroCreate<float32> x.Length
                      Buffer.BlockCopy(x, 0, r, 0, x.Length * sizeof<float32>)
                      r
              // Non-BLAS
              member o.ReshapeCopy_V_MRows(m, x) =
                  if Array.isEmpty x then
                      Array2D.empty<float32>
                  else
  