| 62 let xinit = unsafe { |
62 let xinit = unsafe { |
| 63 //core::intrinsics::assert_inhabited::<[[F; K]; M]>(); |
63 //core::intrinsics::assert_inhabited::<[[F; K]; M]>(); |
| 64 (&x as *const _ as *const [[F; K]; M]).read() |
64 (&x as *const _ as *const [[F; K]; M]).read() |
| 65 }; |
65 }; |
| 66 |
66 |
| 67 std::mem::forget(x); |
67 //std::mem::forget(x); |
| 68 xinit |
68 xinit |
| 69 } |
69 } |
| 70 |
70 |
| 71 /// Gaussian elimination for $Ax=b$, where $A$ and $b$ are both stored in `ab`, |
71 /// Gaussian elimination for $Ax=b$, where $A$ and $b$ are both stored in `ab`, |
| 72 /// $A \in \mathbb{R}^{M \times M}$ and $x, b \in \mathbb{R}^M$. |
72 /// $A \in \mathbb{R}^{M \times M}$ and $x, b \in \mathbb{R}^M$. |