Rewrite LAPACKE by Rust, call LAPACK directly

[termoshtt]

Resolve #111, meta development branch

Goal

Drop lapacke dependency

Changes

• The order of eigenvalue for general matrices changes
• Implement C/F-continuous array handling
• Create temporal work vector in Rust (replacement of LAPACKE_malloc)
• Drop unsafe from lax traits

TODO

• cholesky Cholesky factorization by LAPACK #225
• eig Eigenvalue problem for general matrices based on LAPACK #212
• eigh (in this PR, partially Split eigenvalue and eigenvector tests #217)
• least_square least square by LAPACK #220
• opnorm 5f6a2c3
• qr Rewrite QR using LAPACK #224
• solve LU decomposition based algorithms using LAPACK #213
• solveh Triangular factorization methods for real-symmetric, Hermitian matrix using LAPACK #216
• svd SVD using LAPACK #218
• svddc SVD divide-and-conquer by LAPACK #219
• triangular e9c3481
• tridiagonal Impl tridiagonal by LAPACK #235

[termoshtt]

https://github.com/rust-ndarray/ndarray-linalg/runs/911395642

failures:

---- least_squares_f32_exact_af_bf stdout ----
[ndarray-linalg\tests\least_squares_nrhs.rs:12] &result = LeastSquaresResult {
    singular_values: [1.2196491, 0.4119807, 0.0000654747], shape=[3], strides=[1], layout=CF (0x3), const ndim=1,
    solution: [[-919.2319, 6852.095],
     [235.25002, -1752.6882],
     [-25.910263, 201.41443]], shape=[3, 2], strides=[1, 3], layout=F (0x2), const ndim=2,
    rank: 3,
    residual_sum_of_squares: Some(
        [0.0, 0.0], shape=[2], strides=[1], layout=CF (0x3), const ndim=1,
    ),
}
==== Assetion Failed ====
Expected:
[[0.5019542, 0.5917609],
 [0.7649941, 0.07200937],
 [0.0060762856, 0.9174205]]
Actual:
[[0.5019411, 0.59185815],
 [0.7650014, 0.07195431],
 [0.0059943795, 0.91799533]]
thread 'least_squares_f32_exact_af_bf' panicked at 'Too large deviation in L2-norm: 0.00041470854 > 0.0001 ', /rustc/c7087fe00d2ba919df1d813c040a5d47e43b0fe7\src\libstd\macros.rs:16:9