Boost logo

Ublas :

From: Antonio Martino (antonio.martino_at_[hidden])
Date: 2005-02-06 06:36:55


Here is an example of the unexpected behaviour of project():

CODE (uncomment to get the working versions!):

#include "boost/numeric/ublas/matrix.hpp"

namespace UBLAS = boost::numeric::ublas;

//after run m is filled with the cholesky factorization
template <typename M>
inline
std::size_t choleskyT (M &U)
{
    
   //counter
   std::size_t i = 0;
            
        bool Positive = true;
        std::size_t size = U.size2();
        for (i = 0; i < size; ++ i)
        {
                typename M::value_type t(U (i, i));
        
          /*
      //working version:
      (U.row (i)) (UBLAS::range (i + 1, size)) -= UBLAS::prod (
         (UBLAS::column (U, i)) (UBLAS::range (0, i)),
         UBLAS::conj (U (UBLAS::range (0, i), UBLAS::range (i + 1, size)))
         );
          */
                 
      //compiler error:
      UBLAS::project(U.row (i))(UBLAS::range (i + 1, size)) -= UBLAS::prod (
         UBLAS::project (UBLAS::column (U, i), UBLAS::range (0, i)),
         UBLAS::conj (UBLAS::project (U, UBLAS::range (0, i), UBLAS::range
(i + 1, size)))
         );

                //although it is using UBLAS::project it works properly
(worth noting
                // that there is no assignment to projection!):
                t -= UBLAS::inner_prod (
               UBLAS::conj (UBLAS::project (UBLAS::column (U, i),
UBLAS::range (0, i))),
               UBLAS::project (UBLAS::column (U, i), UBLAS::range (0, i))
            );
                
                if ( t <= 0 ){
                        Positive = false;
                        break;
                }
                
                t = sqrt (t);
                U (i, i) = t;
                
      //working version:
      //(U.row (i)) (UBLAS::range (i + 1, size)) /= t;
      //compiler error:
          UBLAS::project(U.row (i), UBLAS::range (i + 1, size)) /= t;
   }
        
        return i;
}

int main(int argc, char *argv[])
{
    boost::numeric::ublas::matrix<double> A(3,3);
    
    A(0,0) = 3; A(0, 1) = 1; A(0, 2) = 0;
    A(1,0) = 1; A(1, 1) = 3; A(1, 2) = 1;
    A(2,0) = 0; A(2, 1) = 1; A(2, 2) = 3;
    
    choleskyT(A);
}

COMPILATION LOG for VC71 (sorry but my compiler is in italian!)

------ Inizio generazione: Progetto: Chol, Configurazione: Debug Win32
------

Compilazione in corso...
main.cpp
f:\antonio\chol\main.cpp(29) : error C2780:
"boost::numeric::ublas::matrix_indirect<M,boost::numeric::ublas::indirect_ar
ray<A>> boost::numeric::ublas::project(M &,const
boost::numeric::ublas::indirect_array<A> &,const
boost::numeric::ublas::indirect_array<A> &)": argomenti previsti: 3,
forniti: 1
        f:\antonio\boost_1_32_0\boost\numeric\ublas\matrix_proxy.hpp(5394):
vedere la dichiarazione di "boost::numeric::ublas::project"
        f:\antonio\chol\main.cpp(67): vedere il riferimento all'istanza
"size_t choleskyT<boost::numeric::ublas::matrix<T>>(M &)" del template di
funzione in corso di compilazione
        with
        [
            T=double,
            M=boost::numeric::ublas::matrix<double>
        ]
f:\antonio\chol\main.cpp(29) : error C2780: "const
boost::numeric::ublas::matrix_slice<M> boost::numeric::ublas::project(const
boost::numeric::ublas::matrix_slice<M> &,const matrix_slice<M>::slice_type
&,const matrix_slice<M>::slice_type &)": argomenti previsti: 3, forniti: 1
        f:\antonio\boost_1_32_0\boost\numeric\ublas\matrix_proxy.hpp(4442):
vedere la dichiarazione di "boost::numeric::ublas::project"
f:\antonio\chol\main.cpp(29) : error C2780:
"boost::numeric::ublas::matrix_slice<M>
boost::numeric::ublas::project(boost::numeric::ublas::matrix_slice<M>
&,const matrix_slice<M>::slice_type &,const matrix_slice<M>::slice_type &)":
argomenti previsti: 3, forniti: 1
        f:\antonio\boost_1_32_0\boost\numeric\ublas\matrix_proxy.hpp(4437):
vedere la dichiarazione di "boost::numeric::ublas::project"
f:\antonio\chol\main.cpp(29) : error C2780: "const
boost::numeric::ublas::matrix_slice<const M>
boost::numeric::ublas::project(const M &,const matrix_slice<M>::slice_type
&,const matrix_slice<M>::slice_type &)": argomenti previsti: 3, forniti: 1
        f:\antonio\boost_1_32_0\boost\numeric\ublas\matrix_proxy.hpp(4419):
vedere la dichiarazione di "boost::numeric::ublas::project"
f:\antonio\chol\main.cpp(29) : error C2780:
"boost::numeric::ublas::matrix_slice<M> boost::numeric::ublas::project(M
&,const matrix_slice<M>::slice_type &,const matrix_slice<M>::slice_type &)":
argomenti previsti: 3, forniti: 1
        f:\antonio\boost_1_32_0\boost\numeric\ublas\matrix_proxy.hpp(4407):
vedere la dichiarazione di "boost::numeric::ublas::project"
f:\antonio\chol\main.cpp(29) : error C2780: "const
boost::numeric::ublas::matrix_range<M> boost::numeric::ublas::project(const
boost::numeric::ublas::matrix_range<M> &,const matrix_range<M>::range_type
&,const matrix_range<M>::range_type &)": argomenti previsti: 3, forniti: 1
        f:\antonio\boost_1_32_0\boost\numeric\ublas\matrix_proxy.hpp(3484):
vedere la dichiarazione di "boost::numeric::ublas::project"
f:\antonio\chol\main.cpp(29) : error C2780:
"boost::numeric::ublas::matrix_range<M>
boost::numeric::ublas::project(boost::numeric::ublas::matrix_range<M>
&,const matrix_range<M>::range_type &,const matrix_range<M>::range_type &)":
argomenti previsti: 3, forniti: 1
        f:\antonio\boost_1_32_0\boost\numeric\ublas\matrix_proxy.hpp(3479):
vedere la dichiarazione di "boost::numeric::ublas::project"
f:\antonio\chol\main.cpp(29) : error C2780: "const
boost::numeric::ublas::matrix_range<const M>
boost::numeric::ublas::project(const M &,const matrix_range<M>::range_type
&,const matrix_range<M>::range_type &)": argomenti previsti: 3, forniti: 1
        f:\antonio\boost_1_32_0\boost\numeric\ublas\matrix_proxy.hpp(3473):
vedere la dichiarazione di "boost::numeric::ublas::project"
f:\antonio\chol\main.cpp(29) : error C2780:
"boost::numeric::ublas::matrix_range<M> boost::numeric::ublas::project(M
&,const matrix_range<M>::range_type &,const matrix_range<M>::range_type &)":
argomenti previsti: 3, forniti: 1
        f:\antonio\boost_1_32_0\boost\numeric\ublas\matrix_proxy.hpp(3461):
vedere la dichiarazione di "boost::numeric::ublas::project"
f:\antonio\chol\main.cpp(29) : error C2780: "const
boost::numeric::ublas::vector_indirect<V,boost::numeric::ublas::indirect_arr
ay<A>> boost::numeric::ublas::project(const
boost::numeric::ublas::vector_indirect<V,boost::numeric::ublas::indirect_arr
ay<A>> &,const boost::numeric::ublas::indirect_array<A> &)": argomenti
previsti: 2, forniti: 1
        f:\antonio\boost_1_32_0\boost\numeric\ublas\vector_proxy.hpp(1711):
vedere la dichiarazione di "boost::numeric::ublas::project"
f:\antonio\chol\main.cpp(29) : error C2780:
"boost::numeric::ublas::vector_indirect<V,boost::numeric::ublas::indirect_ar
ray<A>>
boost::numeric::ublas::project(boost::numeric::ublas::vector_indirect<V,boos
t::numeric::ublas::indirect_array<A>> &,const
boost::numeric::ublas::indirect_array<A> &)": argomenti previsti: 2,
forniti: 1
        f:\antonio\boost_1_32_0\boost\numeric\ublas\vector_proxy.hpp(1706):
vedere la dichiarazione di "boost::numeric::ublas::project"
f:\antonio\chol\main.cpp(29) : error C2780: "const
boost::numeric::ublas::vector_indirect<V,IA>
boost::numeric::ublas::project(const
boost::numeric::ublas::vector_indirect<V,IA> &,const
vector_indirect<V,IA>::slice_type &)": argomenti previsti: 2, forniti: 1
        f:\antonio\boost_1_32_0\boost\numeric\ublas\vector_proxy.hpp(1701):
vedere la dichiarazione di "boost::numeric::ublas::project"
f:\antonio\chol\main.cpp(29) : error C2780:
"boost::numeric::ublas::vector_indirect<V,IA>
boost::numeric::ublas::project(boost::numeric::ublas::vector_indirect<V,IA>
&,const vector_indirect<V,IA>::slice_type &)": argomenti previsti: 2,
forniti: 1
        f:\antonio\boost_1_32_0\boost\numeric\ublas\vector_proxy.hpp(1696):
vedere la dichiarazione di "boost::numeric::ublas::project"
f:\antonio\chol\main.cpp(29) : error C2780: "const
boost::numeric::ublas::vector_indirect<V,IA>
boost::numeric::ublas::project(const
boost::numeric::ublas::vector_indirect<V,IA> &,const
vector_indirect<V,IA>::range_type &)": argomenti previsti: 2, forniti: 1
        f:\antonio\boost_1_32_0\boost\numeric\ublas\vector_proxy.hpp(1691):
vedere la dichiarazione di "boost::numeric::ublas::project"
f:\antonio\chol\main.cpp(29) : error C2780:
"boost::numeric::ublas::vector_indirect<V,IA>
boost::numeric::ublas::project(boost::numeric::ublas::vector_indirect<V,IA>
&,const vector_indirect<V,IA>::range_type &)": argomenti previsti: 2,
forniti: 1
        f:\antonio\boost_1_32_0\boost\numeric\ublas\vector_proxy.hpp(1686):
vedere la dichiarazione di "boost::numeric::ublas::project"
f:\antonio\chol\main.cpp(29) : error C2780: "const
boost::numeric::ublas::vector_indirect<const
V,boost::numeric::ublas::indirect_array<A>>
boost::numeric::ublas::project(const V &,const
boost::numeric::ublas::indirect_array<A> &)": argomenti previsti: 2,
forniti: 1
        f:\antonio\boost_1_32_0\boost\numeric\ublas\vector_proxy.hpp(1680):
vedere la dichiarazione di "boost::numeric::ublas::project"
f:\antonio\chol\main.cpp(29) : error C2780:
"boost::numeric::ublas::vector_indirect<V,boost::numeric::ublas::indirect_ar
ray<A>> boost::numeric::ublas::project(V &,const
boost::numeric::ublas::indirect_array<A> &)": argomenti previsti: 2,
forniti: 1
        f:\antonio\boost_1_32_0\boost\numeric\ublas\vector_proxy.hpp(1668):
vedere la dichiarazione di "boost::numeric::ublas::project"
f:\antonio\chol\main.cpp(29) : error C2780: "const
boost::numeric::ublas::vector_slice<V> boost::numeric::ublas::project(const
boost::numeric::ublas::vector_slice<V> &,const vector_slice<V>::slice_type
&)": argomenti previsti: 2, forniti: 1
        f:\antonio\boost_1_32_0\boost\numeric\ublas\vector_proxy.hpp(1102):
vedere la dichiarazione di "boost::numeric::ublas::project"
f:\antonio\chol\main.cpp(29) : error C2780:
"boost::numeric::ublas::vector_slice<V>
boost::numeric::ublas::project(boost::numeric::ublas::vector_slice<V>
&,const vector_slice<V>::slice_type &)": argomenti previsti: 2, forniti: 1
        f:\antonio\boost_1_32_0\boost\numeric\ublas\vector_proxy.hpp(1097):
vedere la dichiarazione di "boost::numeric::ublas::project"
f:\antonio\chol\main.cpp(29) : error C2780: "const
boost::numeric::ublas::vector_slice<const V>
boost::numeric::ublas::project(const V &,const vector_slice<V>::slice_type
&)": argomenti previsti: 2, forniti: 1
        f:\antonio\boost_1_32_0\boost\numeric\ublas\vector_proxy.hpp(1091):
vedere la dichiarazione di "boost::numeric::ublas::project"
f:\antonio\chol\main.cpp(29) : error C2780:
"boost::numeric::ublas::vector_slice<V> boost::numeric::ublas::project(V
&,const vector_slice<V>::slice_type &)": argomenti previsti: 2, forniti: 1
        f:\antonio\boost_1_32_0\boost\numeric\ublas\vector_proxy.hpp(1079):
vedere la dichiarazione di "boost::numeric::ublas::project"
f:\antonio\chol\main.cpp(29) : error C2780: "const
boost::numeric::ublas::vector_range<V> boost::numeric::ublas::project(const
boost::numeric::ublas::vector_range<V> &,const vector_range<V>::range_type
&)": argomenti previsti: 2, forniti: 1
        f:\antonio\boost_1_32_0\boost\numeric\ublas\vector_proxy.hpp(560):
vedere la dichiarazione di "boost::numeric::ublas::project"
f:\antonio\chol\main.cpp(29) : error C2780:
"boost::numeric::ublas::vector_range<V>
boost::numeric::ublas::project(boost::numeric::ublas::vector_range<V>
&,const vector_range<V>::range_type &)": argomenti previsti: 2, forniti: 1
        f:\antonio\boost_1_32_0\boost\numeric\ublas\vector_proxy.hpp(555):
vedere la dichiarazione di "boost::numeric::ublas::project"
f:\antonio\chol\main.cpp(29) : error C2780: "const
boost::numeric::ublas::vector_range<const V>
boost::numeric::ublas::project(const V &,const vector_range<V>::range_type
&)": argomenti previsti: 2, forniti: 1
        f:\antonio\boost_1_32_0\boost\numeric\ublas\vector_proxy.hpp(549):
vedere la dichiarazione di "boost::numeric::ublas::project"
f:\antonio\chol\main.cpp(29) : error C2780:
"boost::numeric::ublas::vector_range<V> boost::numeric::ublas::project(V
&,const vector_range<V>::range_type &)": argomenti previsti: 2, forniti: 1
        f:\antonio\boost_1_32_0\boost\numeric\ublas\vector_proxy.hpp(537):
vedere la dichiarazione di "boost::numeric::ublas::project"
f:\antonio\chol\main.cpp(52) : error C2678: "/=" binario: non stato
trovato alcun operatore che accetti un operando sinistro di tipo "const
boost::numeric::ublas::vector_range<V>". anche possibile che non vi siano
conversioni accettabili
        with
        [
            V=const
boost::numeric::ublas::matrix_row<boost::numeric::ublas::matrix<double>>
        ]

Il log di generazione stato salvato in
"file://f:\antonio\chol\Vc\Chol\Debug\BuildLog.htm"
Chol - 26 errore/i, 0 avviso/i

COMPILATION LOG FOR GCC 3.4.2

Compilatore: dev-cpp-3.4.2
Building Makefile: "F:\antonio\chol\Makefile.win"
Esecuzione di make...
mingw32-make.exe -f "F:\antonio\chol\Makefile.win" all
g++.exe -c main.cpp -o main.o -I"C:/Dev-Cpp/include/c++/3.4.2"
-I"C:/Dev-Cpp/include/c++/3.4.2/mingw32"
-I"C:/Dev-Cpp/include/c++/3.4.2/backward" -I"C:/Dev-Cpp/include"
-I..\boost_1_32_0

main.cpp: In function `size_t choleskyT(M&) [with M =
boost::numeric::ublas::matrix<double, boost::numeric::ublas::row_major,
boost::numeric::ublas::unbounded_array<double, std::allocator<double> > >]':
main.cpp:67: instantiated from here
main.cpp:29: error: no matching function for call to
`project(boost::numeric::ublas::matrix_row<boost::numeric::ublas::matrix<dou
ble, boost::numeric::ublas::row_major,
boost::numeric::ublas::unbounded_array<double, std::allocator<double> > >
>)'

main.cpp:67: instantiated from here
main.cpp:52: error: passing `const boost::numeric::ublas::vector_range<const
boost::numeric::ublas::matrix_row<boost::numeric::ublas::matrix<double,
boost::numeric::ublas::row_major,
boost::numeric::ublas::unbounded_array<double, std::allocator<double> > > >
>' as `this' argument of `boost::numeric::ublas::vector_range<E>&
boost::numeric::ublas::vector_range<E>::operator/=(const AT&) [with AT =
double, V = const
boost::numeric::ublas::matrix_row<boost::numeric::ublas::matrix<double,
boost::numeric::ublas::row_major,
boost::numeric::ublas::unbounded_array<double, std::allocator<double> > >
>]' discards qualifiers

../boost_1_32_0/boost/numeric/ublas/functional.hpp: In static member
function `static void boost::numeric::ublas::scalar_divides_assign<T1,
T2>::apply(typename boost::numeric::ublas::scalar_binary_assign_functor<T1,
T2>::argument1_type, typename
boost::numeric::ublas::scalar_binary_assign_functor<T1, T2>::argument2_type)
[with T1 = const double&, T2 = double]':
../boost_1_32_0/boost/numeric/ublas/vector_assign.hpp:126: instantiated
from `void boost::numeric::ublas::indexing_vector_assign_scalar(F, V&, const
T&) [with F = boost::numeric::ublas::scalar_divides_assign<const double&,
double>, V = boost::numeric::ublas::vector_range<const
boost::numeric::ublas::matrix_row<boost::numeric::ublas::matrix<double,
boost::numeric::ublas::row_major,
boost::numeric::ublas::unbounded_array<double, std::allocator<double> > > >
>, T = double]'
../boost_1_32_0/boost/numeric/ublas/vector_assign.hpp:140: instantiated
from `void boost::numeric::ublas::vector_assign_scalar(F, V&, const T&,
boost::numeric::ublas::dense_proxy_tag) [with F =
boost::numeric::ublas::scalar_divides_assign<const double&, double>, V =
boost::numeric::ublas::vector_range<const
boost::numeric::ublas::matrix_row<boost::numeric::ublas::matrix<double,
boost::numeric::ublas::row_major,
boost::numeric::ublas::unbounded_array<double, std::allocator<double> > > >
>, T = double]'
../boost_1_32_0/boost/numeric/ublas/vector_assign.hpp:182: instantiated
from `void boost::numeric::ublas::vector_assign_scalar(F, V&, const T&)
[with F = boost::numeric::ublas::scalar_divides_assign<const double&,
double>, V = boost::numeric::ublas::vector_range<const
boost::numeric::ublas::matrix_row<boost::numeric::ublas::matrix<double,
boost::numeric::ublas::row_major,
boost::numeric::ublas::unbounded_array<double, std::allocator<double> > > >
>, T = double]'

../boost_1_32_0/boost/numeric/ublas/vector_proxy.hpp:203: instantiated
from `boost::numeric::ublas::vector_range<E>&
boost::numeric::ublas::vector_range<E>::operator/=(const AT&) [with AT =
double, V = const
boost::numeric::ublas::matrix_row<boost::numeric::ublas::matrix<double,
boost::numeric::ublas::row_major,
boost::numeric::ublas::unbounded_array<double, std::allocator<double> > >
>]'
main.cpp:52: instantiated from `size_t choleskyT(M&) [with M =
boost::numeric::ublas::matrix<double, boost::numeric::ublas::row_major,
boost::numeric::ublas::unbounded_array<double, std::allocator<double> > >]'
main.cpp:67: instantiated from here
../boost_1_32_0/boost/numeric/ublas/functional.hpp:272: error: assignment of
read-only reference `t1'

mingw32-make.exe: *** [main.o] Error 1

Esecuzione terminata

The problem is reported twice as you may see. Note that there is no test of
matrix projection in assignment!

I have used ublas boost 1.32 distribution for both tests

Bye
Antonio