LUFactorized Implemented?
Chưa có ai nhận issue này.
Đánh giá
- Độ khó
- 4/5
- Thời gian dự kiến
- 3-5 ngày
- Mức phù hợp với người mới
- 25/100
Hướng nghiên cứu
Start with the LUFactorized::factorize_into and solve_inplace implementations shown in the issue, then trace the A::lu and A::solve entry points. Confirm whether the matrix passed to A::lu is mutated into factors and whether the stored ipiv is used; done when the behavior is explained or corrected with a regression test.
Do mô hình lập chỉ mục viết ra từ nội dung của issue.
Mô tả
I was having a look at the code behind LUFactorized and it appears that there is no LU factorization done at all. For example, the factorize_into source is
fn factorize_into(mut self) -> Result<LUFactorized<S>> {
let ipiv = A::lu(self.layout()?, self.as_allocated_mut()?)?;
Ok(LUFactorized { a: self, ipiv })
}
where self here is the original matrix. While it is fine to be storing the original matrix A (perhaps the LU decomposition is lazily evaluated), the solve_inplace code indicates that this is actually not the case at all. Instead it just uses the standard solve function without actually doing any LU decomposition at any point.
fn solve_inplace<'a, Sb>(
&self,
rhs: &'a mut ArrayBase<Sb, Ix1>,
) -> Result<&'a mut ArrayBase<Sb, Ix1>>
where
Sb: DataMut<Elem = A>,
{
assert_eq!(
rhs.len(),
self.a.len_of(Axis(1)),
"The length of `rhs` must be compatible with the shape of the factored matrix.",
);
A::solve(
self.a.square_layout()?,
Transpose::No,
self.a.as_allocated()?,
&self.ipiv,
rhs.as_slice_mut().unwrap(),
)?;
Ok(rhs)
}
Am I missing something? Or is this a temporary implementation until an LU factorization can be properly implemented?
- Ngôn ngữ chính
- Rust
- Star
- 452
- Fork
- 95
- Chỉ số merge pull request
- Không có pull request nào được merge trong 30 ngày
Hướng dẫn đóng góp
Chưa lập chỉ mục được hướng dẫn đóng góp cho kho mã nguồn này
Bắt đầu từ đâu
- Đọc hết issue, rồi đọc hướng dẫn đóng góp của dự án.
- Bình luận trên issue rằng bạn sẽ nhận — tránh hai người làm cùng một việc.
- Fork repository và làm thay đổi trên một nhánh.
- Mở pull request có tham chiếu số hiệu của issue.
Issue khác của rust-ndarray/ndarray-linalg
-
Thin SVD Đang mở
Độ khó 5/5 Hơn một tuần Mức phù hợp với người mới 38/100
rust-ndarray/ndarray-linalg#414 ·
-
Độ khó 5/5 Hơn một tuần Mức phù hợp với người mới 30/100
rust-ndarray/ndarray-linalg#413 · 1 bình luận ·
-
Cyclically Tridiagonal Matrices? Đang mở
Độ khó 5/5 Hơn một tuần Mức phù hợp với người mới 25/100
rust-ndarray/ndarray-linalg#404 ·
-
SIGSEGV on qr decomposition Đang mở
Độ khó 4/5 3-5 ngày Mức phù hợp với người mới 35/100
rust-ndarray/ndarray-linalg#402 · 1 bình luận ·
-
Độ khó 4/5 3-5 ngày Mức phù hợp với người mới 35/100
rust-ndarray/ndarray-linalg#401 · 2 reaction ·
Tất cả issue của rust-ndarray/ndarray-linalg
Issue tương tự
-
Độ khó 2/5 1-3 giờ Mức phù hợp với người mới 75/100
TheLarkInn/aipm#2413 ·
-
documentation
Độ khó 1/5 Dưới một giờ Mức phù hợp với người mới 90/100
alexgorbatchev/simple-ptt#15 ·
-
tooling
Độ khó 2/5 1-3 giờ Mức phù hợp với người mới 75/100
-
todo:ticket
Độ khó 2/5 1-3 giờ Mức phù hợp với người mới 70/100
-
Độ khó 2/5 1-3 giờ Mức phù hợp với người mới 75/100
taikoxyz/taiko-mono#22168 · 1 bình luận ·