Hello.<br><br>I have incorporated cublas bindings on top of blas bindings, such that you can use blas and cublas at once:<br><a href="http://code.google.com/p/asadchev/source/browse/#svn%2Ftrunk%2Fprojects%2Fboost%2Fnumeric%2Fbindings%2Fcublas">http://code.google.com/p/asadchev/source/browse/#svn%2Ftrunk%2Fprojects%2Fboost%2Fnumeric%2Fbindings%2Fcublas</a><br>
actual levelN functions are auto-parsed from blas bindings<br><br>Also several other things are include, such as copy semantics, exceptions, and adaptors<br><br>examples:<br><br>cublas::matrix&lt;double&gt; C1(m, k), C2(n, k);<br>
C1 = H1; // copy from host<br>cublas::matrix_adaptor&lt;double, layout&gt; C3(m, n, dev_ptr);<br>cublas::gemm(1, C1, cublas::trans(C2), 0, C3);<br>cublas::host(H) = C3;  // copy to host<br>C3.clear() // zero<br>cublas::row(C3, 0).clear() // zero row<br>
<br>Would you be willing to incorporate this?<br><br><br>