Пользовательские линейные операторы в ViennaCL

Question

Можно ли передавать пользовательские линейные операторы (например, http://cusplibrary.github.io/classcusp_1_1linear__operator.html) в решатели ViennaCL?

Благодарю.

Answer 1

Да, это возможно для итерационных решателей. Пример можно найти в examples/tutorial/matrix-free.cpp или здесь: http://viennacl.sourceforge.net/doc/matrix-free_8cpp-example.html