All of these are based on the setup and changes described in my previous email.  lapack::gesv is working, from last email.<div><br></div><div>lapack::getrf seems to be working as well!<br></div><div><br></div><div><br></div>
<div>lapack::getri is failing to compile.  Looks like it is the &quot;workspace&quot; issue I had with an older routine.  Is this implemented?  Is there a workaround?</div><div><br></div><div>Code:</div><div><div>ublas::matrix&lt;double, column_major&gt; A(2,2);<br>
</div><div>...</div><div>std::vector&lt;int&gt; ipiv(2);<br></div><div>lapack::getri(A, ipiv);</div><div><br></div><div>Compile errors:</div><div><div>1&gt;C:\working\etk_lib\test_suite\etk\test_ublas.cpp(102): error: no instance of function template &quot;boost::numeric::bindings::lapack::getri&quot; matches the argument list</div>
<div>1&gt;            argument types are: (boost::numeric::ublas::matrix&lt;double, boost::numeric::ublas::column_major, boost::numeric::ublas::unbounded_array&lt;double, std::allocator&lt;QuantLib::Real&gt;&gt;&gt;, std::vector&lt;ptrdiff_t, std::allocator&lt;int&gt;&gt;)</div>
<div>1&gt;  <span class="Apple-tab-span" style="white-space:pre">	</span>lapack::getri(A, ipiv);</div><div>1&gt;  <span class="Apple-tab-span" style="white-space:pre">	</span>^</div></div><div><br></div><div><br></div><div>
<br></div><div><br></div><div><br></div><div>Also, when I change the code to:</div><div>lapack::getri(A, ipiv, lapack::optimal_workspace());<br></div><div><br></div><div>Iit creates the following results, which look like traits problems.:</div>
<div><br></div><div><div>1&gt;C:\working\libraries\boost\boost/numeric/bindings/lapack/computational/getri.hpp(84): error: no operator &quot;&gt;=&quot; matches these operands</div><div>1&gt;            operand types are: boost::numeric::bindings::traits::detail::array&lt;double&gt; &gt;= integer_t</div>
<div>1&gt;          assert( traits::vector_size(work.select(real_type()) &gt;= min_size_work(</div><div>1&gt;          ^</div><div>1&gt;          detected during:</div><div>1&gt;            instantiation of &quot;void boost::numeric::bindings::lapack::getri_impl&lt;ValueType, boost::enable_if&lt;boost::numeric::bindings::traits::is_real&lt;ValueType&gt;, void&gt;::type&gt;::compute(MatrixA &amp;, VectorIPIV &amp;, integer_t={int} &amp;, boost::numeric::bindings::lapack::detail::workspace1&lt;WORK&gt;) [with ValueType=double, MatrixA=boost::numeric::ublas::matrix&lt;double, boost::numeric::ublas::column_major, boost::numeric::ublas::unbounded_array&lt;double, std::allocator&lt;QuantLib::Real={double}&gt;&gt;&gt;,</div>
<div>1&gt;                      VectorIPIV=std::vector&lt;ptrdiff_t={int}, std::allocator&lt;int&gt;&gt;, WORK=boost::numeric::bindings::traits::detail::array&lt;double&gt;]&quot; at line 112</div><div>1&gt;            instantiation of &quot;void boost::numeric::bindings::lapack::getri_impl&lt;ValueType, boost::enable_if&lt;boost::numeric::bindings::traits::is_real&lt;ValueType&gt;, void&gt;::type&gt;::compute(MatrixA &amp;, VectorIPIV &amp;, integer_t={int} &amp;, boost::numeric::bindings::lapack::optimal_workspace) [with ValueType=double, MatrixA=boost::numeric::ublas::matrix&lt;double, boost::numeric::ublas::column_major, boost::numeric::ublas::unbounded_array&lt;double, std::allocator&lt;QuantLib::Real={double}&gt;&gt;&gt;,</div>
<div>1&gt;                      VectorIPIV=std::vector&lt;ptrdiff_t={int}, std::allocator&lt;int&gt;&gt;]&quot; at line 180</div><div>1&gt;            instantiation of &quot;integer_t={int} boost::numeric::bindings::lapack::getri(MatrixA &amp;, VectorIPIV &amp;, Workspace) [with MatrixA=boost::numeric::ublas::matrix&lt;double, boost::numeric::ublas::column_major, boost::numeric::ublas::unbounded_array&lt;double, std::allocator&lt;QuantLib::Real={double}&gt;&gt;&gt;, VectorIPIV=std::vector&lt;ptrdiff_t={int}, std::allocator&lt;int&gt;&gt;, Workspace=boost::numeric::bindings::lapack::optimal_workspace]&quot; at line 102 of &quot;C:\working\etk_lib\test_suite\etk\test_ublas.cpp&quot;</div>
<div>1&gt;</div><div>1&gt;C:\working\libraries\boost\boost/numeric/bindings/traits/vector_traits.hpp(51): error: name followed by &quot;::&quot; must be a class or namespace name</div><div>1&gt;    struct vector_detail_traits : default_vector_traits&lt;VType, typename VType::value_type &gt; </div>
<div>1&gt;                                                                        ^</div><div>1&gt;          detected during:</div><div>1&gt;            instantiation of class &quot;boost::numeric::bindings::traits::vector_detail_traits&lt;VIdentifier, VType&gt; [with VIdentifier=&lt;error-type&gt;, VType=&lt;error-type&gt;]&quot; at line 59</div>
<div>1&gt;            instantiation of class &quot;boost::numeric::bindings::traits::vector_traits&lt;V&gt; [with V=&lt;error-type&gt;]&quot; at line 77</div><div>1&gt;            instantiation of &quot;ptrdiff_t={int} boost::numeric::bindings::traits::vector_size(V &amp;) [with V=&lt;error-type&gt;]&quot; at line 84 of &quot;C:\working\libraries\boost\boost/numeric/bindings/lapack/computational/getri.hpp&quot;</div>
<div>1&gt;            instantiation of &quot;void boost::numeric::bindings::lapack::getri_impl&lt;ValueType, boost::enable_if&lt;boost::numeric::bindings::traits::is_real&lt;ValueType&gt;, void&gt;::type&gt;::compute(MatrixA &amp;, VectorIPIV &amp;, integer_t={int} &amp;, boost::numeric::bindings::lapack::detail::workspace1&lt;WORK&gt;) [with ValueType=double, MatrixA=boost::numeric::ublas::matrix&lt;double, boost::numeric::ublas::column_major, boost::numeric::ublas::unbounded_array&lt;double, std::allocator&lt;QuantLib::Real={double}&gt;&gt;&gt;,</div>
<div>1&gt;                      VectorIPIV=std::vector&lt;ptrdiff_t={int}, std::allocator&lt;int&gt;&gt;, WORK=boost::numeric::bindings::traits::detail::array&lt;double&gt;]&quot; at line 112 of &quot;C:\working\libraries\boost\boost/numeric/bindings/lapack/computational/getri.hpp&quot;</div>
<div>1&gt;            instantiation of &quot;void boost::numeric::bindings::lapack::getri_impl&lt;ValueType, boost::enable_if&lt;boost::numeric::bindings::traits::is_real&lt;ValueType&gt;, void&gt;::type&gt;::compute(MatrixA &amp;, VectorIPIV &amp;, integer_t={int} &amp;, boost::numeric::bindings::lapack::optimal_workspace) [with ValueType=double, MatrixA=boost::numeric::ublas::matrix&lt;double, boost::numeric::ublas::column_major, boost::numeric::ublas::unbounded_array&lt;double, std::allocator&lt;QuantLib::Real={double}&gt;&gt;&gt;,</div>
<div>1&gt;                      VectorIPIV=std::vector&lt;ptrdiff_t={int}, std::allocator&lt;int&gt;&gt;]&quot; at line 180 of &quot;C:\working\libraries\boost\boost/numeric/bindings/lapack/computational/getri.hpp&quot;</div>
<div>1&gt;            instantiation of &quot;integer_t={int} boost::numeric::bindings::lapack::getri(MatrixA &amp;, VectorIPIV &amp;, Workspace) [with MatrixA=boost::numeric::ublas::matrix&lt;double, boost::numeric::ublas::column_major, boost::numeric::ublas::unbounded_array&lt;double, std::allocator&lt;QuantLib::Real={double}&gt;&gt;&gt;, VectorIPIV=std::vector&lt;ptrdiff_t={int}, std::allocator&lt;int&gt;&gt;, Workspace=boost::numeric::bindings::lapack::optimal_workspace]&quot; at line 102 of &quot;C:\working\etk_lib\test_suite\etk\test_ublas.cpp&quot;</div>
</div></div>