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