Giải pháp hệ thống tuyến tính nhanh cho D?

Question

Tôi có thể nhận được giải pháp hệ thống tuyến tính nhanh được viết bằng D ở đâu? Nó sẽ có thể để có một ma trận vuông Một và một vector b và giải quyết các phương trình Axe = b cho b và lý tưởng, cũng thực hiện đảo ngược rõ ràng về Một. Tôi có một bản thân mình viết, nhưng nó khá chậm, có lẽ vì nó hoàn toàn là bộ nhớ cache ngây thơ. Tuy nhiên, đối với trường hợp sử dụng của tôi, tôi cần một cái gì đó với tuyệt đối không thể thương lượng yêu cầu, sau tức là nếu nó không đáp ứng các, sau đó tôi không có thoả chăm sóc như thế nào tốt, nếu không thì là:

1. Phải là miền công cộng được cấp phép, Giấy phép nâng cao hoặc một số giấy phép cho phép tương tự. Lý tưởng nhất là không nên yêu cầu phân bổ trong các tệp nhị phân (tức là không phải BSD), mặc dù điểm này có phần thương lượng.

2. Phải được viết bằng D thuần túy hoặc dễ dàng dịch sang thuần túy D. Mã Fortran không thể sửa chữa (tức là LAPACK) không phải là câu trả lời hay cho dù nó có nhanh như thế nào.

3. Phải được tối ưu hóa cho các hệ thống lớn (nghĩa là 1000). Tôi không muốn một cái gì đó được thiết kế cho các lập trình viên trò chơi để giải quyết các ma trận 4x4 thực sự, rất nhanh.

4. Không được liên kết chặt chẽ với một thư viện khổng lồ mà tôi không cần.

Chỉnh sửa: Lý do cho các yêu cầu dường như điên rồ này là tôi cần mã này cho thư viện nguồn mở được cấp phép được phép mà tôi không muốn có bất kỳ phụ thuộc bên thứ ba nào.

Answer 1

Nếu bạn không thích mã Fortran, một thư viện ma trận dày đặc C++ nhanh có hỗ trợ đa lõi khiêm tốn, mã được viết tốt và giao diện người dùng tốt là Eigen. Nó nên được đơn giản để dịch mã của nó để D (hoặc để có một số thuật toán từ nó).

Và bây giờ là "suy nghĩ về yêu cầu của bạn": có một lý do tại sao "mọi người" (Mathematica, Matlab, Maple, SciPy, GSL, R, ...) sử dụng ATLAS/LAPACK, UMFPACK, PARDISO, CHOLMOD v.v ... Đó là công việc khó khăn để viết nhanh, đa luồng, bộ nhớ hiệu quả, xách tay và số ổn định ma trận giải quyết (tin tưởng tôi, tôi đã cố gắng). Rất nhiều công việc khó khăn này đã đi vào ATLAS và phần còn lại.

Vì vậy, cách tiếp cận của tôi sẽ là viết các ràng buộc cho thư viện có liên quan tùy thuộc vào loại ma trận của bạn và liên kết từ D với giao diện C. Có lẽ các ràng buộc trong multiarray là đủ (tôi đã không cố gắng). Nếu không, tôi khuyên bạn nên xem xét một thư viện C++ khác, cụ thể là uBlas và bindings tương ứng cho các ý tưởng.