commit
d6cca60837
|
@ -294,7 +294,8 @@ impl<N: Scalar + Ring, D: DimName, S: StorageMut<N, D, D>> SquareMatrix<N, D, S>
|
||||||
{
|
{
|
||||||
for i in 0..D::dim() {
|
for i in 0..D::dim() {
|
||||||
for j in 0..D::dim() - 1 {
|
for j in 0..D::dim() - 1 {
|
||||||
self[(j, i)] += shift[j] * self[(D::dim() - 1, i)];
|
let add = shift[j] * self[(D::dim() - 1, i)];
|
||||||
|
self[(j, i)] += add;
|
||||||
}
|
}
|
||||||
}
|
}
|
||||||
}
|
}
|
||||||
|
|
|
@ -58,7 +58,7 @@ impl<N: Scalar + Zero, R: Dim, C: Dim, S: Storage<N, R, C>> Matrix<N, R, C, S> {
|
||||||
for j in 0..ncols.value() {
|
for j in 0..ncols.value() {
|
||||||
// FIXME: use unchecked column indexing
|
// FIXME: use unchecked column indexing
|
||||||
let mut res = res.column_mut(j);
|
let mut res = res.column_mut(j);
|
||||||
let mut src = self.column(j);
|
let src = self.column(j);
|
||||||
|
|
||||||
for (destination, source) in irows.clone().enumerate() {
|
for (destination, source) in irows.clone().enumerate() {
|
||||||
unsafe {
|
unsafe {
|
||||||
|
|
Loading…
Reference in New Issue