Re: [ublas] Triangular Matrix Question

From: Manoj Rajagopalan (rmanoj_at_[hidden])
Date: 2006-03-30 10:30:04

I have often come across such errors in compiling my code. Often, 
including one or more other headers solves the problem.

try including one or more of the following:
expression_types.hpp
matrix_proxy.hpp
matrix_expression.hpp
traits.hpp

you could try other header files that you suspect contain typedefs of 
closure_type (there are so many!).

This is only a hack though - the best I have come up with.

cheers
Manoj

  James N. Knight wrote:
> Hello all.  I'm trying to initialize a triangular matrix using the 
> following code.
> 
>      #include <boost/numeric/ublas/matrix.hpp>
>      #include <boost/numeric/ublas/triangular.hpp>
> 
>      using namespace boost::numeric::ublas;
> 
>      int main(){
> 
>         scalar_matrix<double> sm(5,5,1.0);
> 
>         triangular_adaptor<scalar_matrix<double>,upper> tma(sm);
> 
>         triangular_matrix<double,upper> tm(tma);
>      }
> 
> I'm compiling against the latest (3/23/06) cvs head of boost using 
> gcc4.0 on a linux machine.  A number of compiler errors are generated. 
> Here are the first few (I've modified them a bit for easier reading).
> 
> .../ublas/triangular.hpp: In instantiation of 
> ‘ublas::triangular_adaptor<ublas::scalar_matrix<double>, 
> ublas::basic_upper<long unsigned int> >’:
> test_tria.cpp:10:   instantiated from here
> .../ublas/triangular.hpp:933: error: no type named ‘closure_type’ in 
> ‘class ublas::scalar_matrix<double>’
> .../ublas/triangular.hpp:965: error: no type named ‘closure_type’ in 
> ‘class ublas::scalar_matrix<double>’
> .../ublas/triangular.hpp:969: error: no type named ‘closure_type’ in 
> ‘class ublas::scalar_matrix<double>’
> .../ublas/triangular.hpp: In instantiation of 
> ‘ublas::triangular_adaptor<ublas::scalar_matrix<double>, 
> ublas::basic_upper<long unsigned int> >’:
> test_tria.cpp:10:   instantiated from here
> .../ublas/triangular.hpp:1771: error: no type named ‘closure_type’ in 
> ‘class ublas::scalar_matrix<double>’
> .../ublas/triangular.hpp: In constructor ‘ublas::triangular_adaptor<M, 
> TRI>::triangular_adaptor(M&) [with M = ublas::scalar_matrix<double>, TRI 
> = ublas::basic_upper<long unsigned int>]’:
> test_tria.cpp:10:   instantiated from here
> .../ublas/triangular.hpp:947: error: using invalid field 
> ‘ublas::triangular_adaptor<M, TRI>::data_’
> 
> 
> It goes on from there.  Any ideas on why this is happening? Is this an 
> invalid use of triangular_adaptor?  I've also done it this way.
> 
> triangular_matrix<double,upper> tm(5,5);
> tm.assign(scalar_matrix<double>(1.0));
> 
> This compiles, but produces an external_logic error (as expected) when I 
> compile and run without -DNDEBUG.  It does appear to work when I compile 
> without debugging though.
> 
> James
> _______________________________________________
> ublas mailing list
> ublas_at_[hidden]
> http://lists.boost.org/mailman/listinfo.cgi/ublas
>