Hafif / kurulum gerektirmeyen C veya C ++ tabanlı yoğun doğrusal cebir çözücü için öneriler

Programlamamın çoğu, kendi kullanımım için C'de bir kerelik araştırma kodlarıdır. Yakın işbirlikçiler dışında hiçbir kod dağıtmadım. Bilimsel bir dergide yayınladığım bir algoritma geliştirdim. Makaleye çevrimiçi ek kaynak kodu ve belki de yürütülebilir kodu sağlamak istiyorum. Bir meslektaşım, C ++ (ack!) Yazmamı gerektiren ve küçük yoğun doğrusal sistemleri çözmem gereken algoritmaya genelleme yapmamı istedi. Algoritma için bir kullanıcı tabanı almayı başarırsam, kısmen giriş çubuğunun kullanımı düşük olduğu için (zeminde olduğu gibi) olacaktır. Potansiyel kullanıcılar, kodu kullanmak için kitaplık vb. Yüklemez. Kodun tamamen bağımsız olmasını ve herhangi bir lisansla sınırsız olmasını istiyorum. Golub ve van Loan'dan bir şey alarak kendi çözücümü yazabilirim ama orada başka birinin yazmış olduğu bir vanilya çözücüsünü kullanmayı tercih ederim. Öneriler takdir edildi. Teşekkürler!

Lapack arayüzünü ihtiyacınız olan fonksiyona tam olarak kopyalamanızı öneririm, muhtemelen ihtiyacınız var dgesv. Bu şekilde, Lapack'in kurulu olduğu insanlar sadece ona bağlanabilir ve sadece çalışır. Lapack yüklü olmayan kişiler için, bu işlevle ilgili kendi basit uygulamanızı sağlar veya muhtemelen diğerlerinin önerdiği gibi Eigen veya FLENS kullanarak uygularsınız.

Fortran topraklarında, Lapack kütüphanesi öyle bir standarttır ki, çoğu insan kendi uygulamalarını sağlamak yerine onu kullanır ve işte budur.

Birçok kişinin bilimsel hesaplama yapmaya başladığında yaptığı çok erken bir hata, tüm kodunuzu aynı dilde yazmanız gerektiğini varsayar. Bence bu, derlenmiş programların aynı derleyicinin bile sürümleri arasında birbirleriyle nasıl iletişim kuracağı net olmadığında, büyük ölçüde tarihsel nedenlerden kaynaklandığını düşünüyorum. Bununla birlikte, bu durumda, yine de C ++ kullanacaksanız, ihtiyaçlarınızı karşılayabilecek çok iyi birkaç C ++ yalnızca başlık şablon kütüphanesi vardır.

Kodunuzu akademik nedenlerle dağıttığınız ve yoğun bir lineer cebir çözücüyü kodunuza gömmek istediğiniz için, Eigen'i düşünmenizi şiddetle tavsiye ederim . Eigen, Mozilla Kamu Lisansı kapsamında lisanslanmıştır ve yalnızca başlıklı bir kütüphanedir. Bu, Eigen'i kodunuzla kaynak formda dağıtabileceğiniz anlamına gelir (bu, kodunuza herhangi bir lisans kısıtlaması getirmez) ve son derece verimli yoğun doğrusal çözücüler de dahil olmak üzere genel özelliklerine erişim elde edeceğiniz anlamına gelir. GertVdE'den bahsedildiği gibi, başka seçenekleriniz de var .

Doğrusal denklem sistemleri için güvenilir bir çözümleyici istiyorsanız FLENS'i tavsiye ederim . LAPACK'in tam olarak yeniden uygulanmasını içerir (tek iş parçacıklı bir BLAS uygulaması kullanılıyorsa, LAPACK ile aynı yuvarlaklaştırma hatalarını bile üretir). Bu, tüm FLENS-LAPACK işlevleri için geçerlidir (yardımcı program işlevleriyle birlikte yaklaşık 100 rutin).

FLENS bir BSD Lisansı altındadır ve bu nedenle tescilli ürünlere dahil edilmesine izin verir.

FLENS yalnızca başlıktır ve yalnızca bir FLENS alt kümesine ihtiyacınız varsa, size yalnızca ihtiyacınız olan işlevleri içeren soyulmuş bir sürüm verebilirim. FLENS kendi referans BLAS uygulaması ile birlikte gelir. Ancak isteğe bağlı olarak kullanıcılarınız ATLAS, OpenBLAS veya GotoBALS gibi optimize edilmiş BLAS kitaplıklarına bağlanabilir. Büyük matrisler için bu, Eigen'e kıyasla yaklaşık% 40'lık bir performans kazancı sağlar.

Ve evet, Eigen de sonuçlarını kontrol etmek için LAPACK test paketini kullanıyor. Bunu 3 işlev için yaparlar (Lu, Cholesky ve Özdeğerler / - simetrik bir matrisin vektörleri). Bununla birlikte, simetrik olmayan bir matrisin özdeğer / vektörlerini hesaplamaları LAPACK test takımında başarısız olacaktır.

Feragatname: Evet, FLENS benim bebeğim! Bu yaklaşık% 95 kodladım ve her kod satırı değdi demektir.