Я пытаюсь использовать CUSP в качестве внешнего линейного решателя для Mathematica для использования мощности графического процессора.
Вот веб-страница проекта CUSP . Я прошу несколько предложений о том, как мы можем интегрировать CUSP с Mathematica. Я уверен, что многим из вас здесь будет интересно обсудить это. Я думаю, что написать входную матрицу и затем передать ее в программу CUSP - не самый лучший способ. Использование Mathematica LibrarayFunctionLoad будет лучшим способом передачи входной матрицы в решатель на основе графического процессора на лету. Каким будет способ поставки матрицы и правой матрицы непосредственно из Mathematica?
Вот фрагмент кода CUSP.
#include <cusp/hyb_matrix.h>
#include <cusp/io/matrix_market.h>
#include <cusp/krylov/cg.h>
int main(void)
{
// create an empty sparse matrix structure (HYB format)
cusp::hyb_matrix<int, float, cusp::device_memory> A;
// load a matrix stored in MatrixMarket format
cusp::io::read_matrix_market_file(A, "5pt_10x10.mtx");
// allocate storage for solution (x) and right hand side (b)
cusp::array1d<float, cusp::device_memory> x(A.num_rows, 0);
cusp::array1d<float, cusp::device_memory> b(A.num_rows, 1);
// solve the linear system A * x = b with the Conjugate Gradient method
cusp::krylov::cg(A, x, b);
return 0;
}
Этот вопрос дает нам возможность обсудить возможности компиляции Mathematica 8. Также можно вызвать тему интерфейса mathlink MMA. Я надеюсь, что люди здесь находят эту проблему достойной и достаточно интересной, чтобы обдумать.
BR