Nabble - MTL Dev

traversal and update of sparse matrix

2009-10-18T09:33:37Z

Hi there,

I want to traverse a sparse matrix and update its elements. I have slightly modified the iteration example from the tutorial, and the traversal works fine. But if I add an inserter, the programm terminates with

Assertion failed: ma.starts[major+1] <= ma.my_nnz, file c:\programme\boost\boost_1_39\boost\numeric\mtl\matrix\compressed2d.hpp, line 209

If I put the inserter in an own subroutine, everything works fine, but the update process is very slow (probably due to the repeated construction of the inserter).
Has someone an idea how to solve this problem?

The source code is as follows:

template <typename Matrix>
void f(Matrix& m)
{
using traits::range_generator; using traits::range::min;

// Choose between row and column traversal
typedef typename min<range_generator<tag::row, Matrix>, range_generator<tag::col, Matrix> >::type range_type;
range_type my_range;

// Type of outer cursor
typedef typename range_type::type c_type;
// Type of inner cursor
typedef typename traits::range_generator<tag::nz, c_type>::type ic_type;

// Define the property maps
typename traits::row<Matrix>::type row(m);
typename traits::col<Matrix>::type col(m);
typename traits::const_value<Matrix>::type value(m);

// Type of m's elements
typedef typename Collection<Matrix>::value_type value_type;
mtl::matrix::inserter<Matrix, update_store<value_type> > ins(m);

for (c_type cursor(my_range.begin(m)), cend(my_range.end(m)); cursor != cend; ++cursor)
for (ic_type icursor(begin<tag::nz>(cursor)), icend(end<tag::nz>(cursor)); icursor != icend; ++icursor){
int r=row(*icursor);
int c=col(*icursor);
double insval=6.789; // or something else
ins[r][c] << insval;
}
}

Thanks in advance,

Burkhard

Re: MTL: compressed2D insertion

2009-10-16T08:09:54Z

Hi Elmar,

thanks for your help!

Since I am a C++ novice, this will take its time... but I am working on it! Since the adjacency matrix is of a graph that consits of cliques only, the building of the blocks will come in a natural way.

Best,
Burkhard

Hess, Elmar wrote:

Burkhard

I think you should use block-wise insertion (see http://www.osl.iu.edu/research/mtl/mtl4/doc/matrix_insertion.html ). This is much faster.
Nevertheless since you are reading your matrix data from a file you will need to think about how to block it first.

Regards,
Elmar

-----Original Message-----
From: mtl-devel-bounces@osl.iu.edu [mailto:mtl-devel-bounces@osl.iu.edu] On Behalf Of Burkhard Möller
Sent: 16 October 2009 11:52
To: mtl-devel@osl.iu.edu
Subject: MTL: compressed2D insertion

Hi there,

I have an adjacency matrix of about 30000x30000 with ~200000 nonzero
entries and I would like to use it with MTL. The matrix is in a file in
matrix market format and I have slightly modified the first Matrix
Insertion example to read it. The insertion of the first 100 elements
already takes over 15 sec time, what am I doing wrong? I am using Visual
Studio 2008 Express.

The source code is as follows:

// File: insert.cpp

#include <iostream>
#include <boost/numeric/mtl/mtl.hpp>

using namespace mtl;

template <typename Matrix>
void insval(Matrix& m, int index1, int index2, double val)
{
// Matrices are not initialized by default
m= 0.0;

// Create inserter for matrix m
matrix::inserter<Matrix> ins(m);

// Insert value
ins[index1][index2] << val;

}

int main(int argc, char* argv[])
{
compressed2D<double> A(30000, 30000);

FILE* matrixfile;

matrixfile=fopen("A_cpp.txt","r");

int count=0, index1=0, index2=0;
double value=0;

clock_t start, finish;
double duration;

start = clock();

while(fscanf(matrixfile, "%d %d %Lg", &index1, &index2, &value)!=EOF
&& count<100){// insert first 100 values only
count++;
insval(A, index1-40360, index2-40360, value); // smallest index
is 40360...
}

finish = clock();

duration = (double)(finish - start) / CLOCKS_PER_SEC;
printf( "%2.1f seconds\n", duration );

return 0;
}

Thanks for your time!

Rechnan
_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/

Forward this message to spam@mailcontrol.com to report this email as spam.

Visit our website at http://www.halcrow.com

------------------------------------------------------------------------
The contents of this email are confidential, for the sole use
of the intended recipient at the email address to which it has
been addressed and do not give rise to any binding legal
obligation upon Halcrow companies unless subsequently confirmed
on headed business notepaper sent by fax, letter or as an email
attachment. Whilst reasonable care has been taken to avoid virus
transmission, no responsibility for viruses is taken and it is
your responsibility to carry out such checks as you feel
appropriate. Emails supplied are as found and there's no
guarantee that the messages contained within the body of the
email have not been edited after receipt. If you receive this
email in error, please contact the sender immediately and delete
the message from your system.
Thank you.
-------------------------------------------------------------------------

_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/

Re: MTL: compressed2D insertion

2009-10-16T04:27:48Z

Burkhard

I think you should use block-wise insertion (see http://www.osl.iu.edu/research/mtl/mtl4/doc/matrix_insertion.html ). This is much faster.
Nevertheless since you are reading your matrix data from a file you will need to think about how to block it first.

Regards,
Elmar

-----Original Message-----
From: mtl-devel-bounces@... [mailto:mtl-devel-bounces@...] On Behalf Of Burkhard Möller
Sent: 16 October 2009 11:52
To: mtl-devel@...
Subject: MTL: compressed2D insertion

Hi there,

I have an adjacency matrix of about 30000x30000 with ~200000 nonzero
entries and I would like to use it with MTL. The matrix is in a file in
matrix market format and I have slightly modified the first Matrix
Insertion example to read it. The insertion of the first 100 elements
already takes over 15 sec time, what am I doing wrong? I am using Visual
Studio 2008 Express.

The source code is as follows:

// File: insert.cpp

#include <iostream>
#include <boost/numeric/mtl/mtl.hpp>

using namespace mtl;

template <typename Matrix>
void insval(Matrix& m, int index1, int index2, double val)
{
// Matrices are not initialized by default
m= 0.0;

// Create inserter for matrix m
matrix::inserter<Matrix> ins(m);

// Insert value
ins[index1][index2] << val;

}

int main(int argc, char* argv[])
{
compressed2D<double> A(30000, 30000);

FILE* matrixfile;

matrixfile=fopen("A_cpp.txt","r");

int count=0, index1=0, index2=0;
double value=0;

clock_t start, finish;
double duration;

start = clock();

while(fscanf(matrixfile, "%d %d %Lg", &index1, &index2, &value)!=EOF
&& count<100){// insert first 100 values only
count++;
insval(A, index1-40360, index2-40360, value); // smallest index
is 40360...
}

finish = clock();

duration = (double)(finish - start) / CLOCKS_PER_SEC;
printf( "%2.1f seconds\n", duration );

return 0;
}

Thanks for your time!

Rechnan
_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/

Forward this message to spam@... to report this email as spam.

Visit our website at http://www.halcrow.com

------------------------------------------------------------------------
The contents of this email are confidential, for the sole use
of the intended recipient at the email address to which it has
been addressed and do not give rise to any binding legal
obligation upon Halcrow companies unless subsequently confirmed
on headed business notepaper sent by fax, letter or as an email
attachment. Whilst reasonable care has been taken to avoid virus
transmission, no responsibility for viruses is taken and it is
your responsibility to carry out such checks as you feel
appropriate. Emails supplied are as found and there's no
guarantee that the messages contained within the body of the
email have not been edited after receipt. If you receive this
email in error, please contact the sender immediately and delete
the message from your system.
Thank you.
-------------------------------------------------------------------------

_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/

MTL: compressed2D insertion

2009-10-16T03:52:10Z

Hi there,

I have an adjacency matrix of about 30000x30000 with ~200000 nonzero
entries and I would like to use it with MTL. The matrix is in a file in
matrix market format and I have slightly modified the first Matrix
Insertion example to read it. The insertion of the first 100 elements
already takes over 15 sec time, what am I doing wrong? I am using Visual
Studio 2008 Express.

The source code is as follows:

// File: insert.cpp

#include <iostream>
#include <boost/numeric/mtl/mtl.hpp>

using namespace mtl;

template <typename Matrix>
void insval(Matrix& m, int index1, int index2, double val)
{
// Matrices are not initialized by default
m= 0.0;

// Create inserter for matrix m
matrix::inserter<Matrix> ins(m);

// Insert value
ins[index1][index2] << val;

}

int main(int argc, char* argv[])
{
compressed2D<double> A(30000, 30000);

FILE* matrixfile;

matrixfile=fopen("A_cpp.txt","r");

int count=0, index1=0, index2=0;
double value=0;

clock_t start, finish;
double duration;

start = clock();

while(fscanf(matrixfile, "%d %d %Lg", &index1, &index2, &value)!=EOF
&& count<100){// insert first 100 values only
count++;
insval(A, index1-40360, index2-40360, value); // smallest index
is 40360...
}

finish = clock();

duration = (double)(finish - start) / CLOCKS_PER_SEC;
printf( "%2.1f seconds\n", duration );

return 0;
}

Thanks for your time!

Rechnan
_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/

compressed2D insertion

2009-10-16T03:48:30Z

Hi there,

I have an adjacency matrix of about 30000x30000 with ~200000 nonzero entries and I would like to use it with MTL. The matrix is in a file in matrix market format and I have slightly modified the first Matrix Insertion example to read it. The insertion of the first 100 elements already takes over 15 sec time, what am I doing wrong? I am using Visual Studio 2008 Express.

The source code is as follows:

// File: insert.cpp

#include <iostream>
#include <boost/numeric/mtl/mtl.hpp>

using namespace mtl;

template <typename Matrix>
void insval(Matrix& m, int index1, int index2, double val)
{
// Matrices are not initialized by default
m= 0.0;

// Create inserter for matrix m
matrix::inserter<Matrix> ins(m);

// Insert value
ins[index1][index2] << val;

}

int main(int argc, char* argv[])
{
compressed2D<double> A(30000, 30000);

FILE* matrixfile;

matrixfile=fopen("A_cpp.txt","r");

int count=0, index1=0, index2=0;
double value=0;

clock_t start, finish;
double duration;

start = clock();

while(fscanf(matrixfile, "%d %d %Lg", &index1, &index2, &value)!=EOF && count<100){// insert first 100 values only
count++;
insval(A, index1-40360, index2-40360, value); // smallest index is 40360...
}

finish = clock();

duration = (double)(finish - start) / CLOCKS_PER_SEC;
printf( "%2.1f seconds\n", duration );

return 0;
}

Thanks for your time!

Rechnan

MTL: matrix library for relational joins

2009-01-28T06:22:12Z

I'm looking for a lib to use to do relational joins (cross products
with restrictions on a subset of the columns) and was wondering if MTL
was appropriate.

While poking through the faq, I noticed a couple broken fragments:
<A NAME="#howtoinvert"></A>
<A NAME="#howtosolve"></A>
The table of contents should work again if you remove the '#'s from
those anchors. You can use an html validator to find such problems:
http://validator.w3.org/check?uri=http%3A%2F%2Fwww.osl.iu.edu%2Fresearch%2Fmtl%2Ffaq.php3&charset=(detect+automatically)&doctype=Inline&group=0
but you'll have a bunch to tackle. You can always run it through tidy,
perhaps with -asxml if you'd like xhtml output.
--
-eric

office: +1.617.258.5741 32-G528, MIT, Cambridge, MA 02144 USA
mobile: +1.617.599.3509

(eric@...)
Feel free to forward this message to any list for any purpose other than
email address distribution.

There are subtle nuances encoded in font variation and clever layout
which can only be seen by printing this message on high-clay paper.

_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/

signature.asc (492 bytes) Download Attachment

MTL: Examples work correctly with Anthony's defines under MS Visual Studio 2005 Express Edition

2009-01-11T22:39:03Z

Problem (which was described in my message on 2008-11-03) can be really resolved by adding two preprocessor definitions:

#define _HAS_ITERATOR_DEBUGGING 0
#define _SECURE_SCL 0

before including of any STL, MTL or ITL files.
After this examples array2D.cc, sparse_copy.cc, sparse_matrix.cc, sparse_mult.cc work correctly.

Thanks a lot to Anthony Lock. Now I can continue my research work with MTL + ITL.

_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/

MTL: 'debug assertion failed' -- work-around

2009-01-10T13:10:01Z

One way to get around the 'debug assertion failed' error I encountered when running ITL on MS Visual Studio 2008 is given below. [See my posting of 2009-01-05 and that of Sergey Avdeev on 2008-11-03.]

In order to get ITL to work, I have included the following preprocessor macros:

1. _HAS_ITERATOR_DEBUGGING=0

2. _SECURE_SCL=0

In the Properties page for the project, add these to: 'C/C++ -> Preprocessor -> Preprocessor definitions'

This isn't ideal, but it works.

_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/

MTL: QMR iterative method using Visual Studio 2008 Express Edition

2009-01-05T12:15:28Z

I am having difficulty trying to run the QMR iterative method, in ITL, using Visual Studio 2008 Express Edition. The method is the first example at:

http://www.osl.iu.edu/research/itl/examples.php3

When I build my 'solution', it compiles without any problems, but when I try to run the code, I get a message saying "debug assertion failed". In the same pop-up message it says "Expression: vector iterator is not dereferencable". The error seems to come from:

C:\program files\microsoft visual studio 9.0\vc\include\vector

When I click 'retry', to debug, a break point is introduced in the call stack window at:

> ITL_Test.exe!std::_Vector_const_iterator<int,std::allocator<int> >::operator*() Line 98 + 0x14 bytes C++

Any help would be appreciated,

Anthony.

_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/

MTL: MTL and Visual Studio 2008 (Express Edition)

2009-01-03T03:05:26Z

I am currently trying to move some C++ code from Linux to Windows; the code uses MTL for all of the linear algebra. As a first test, I have tried to use MS Visual Studio 2008 (express edition) to run the simple example given at:

http://www.osl.iu.edu/research/mtl/examples/euclid_norm.cc

However, I end up with 54 error messages and 11 warnings, all of which seem to stem from the header file 'light_matrix.h'. Please find below the test, together with the error messages.

Any help would be appreciated,

Anthony.

****************************

test

****************************

// 02Jan09

// Taken from: http://www.osl.iu.edu/research/mtl/examples/euclid_norm.cc

// -*- c++ -*-

// $COPYRIGHT$

//===========================================================================

#include "mtl/mtl.h"

#include "mtl/utils.h"

#include "mtl/linalg_vec.h"

#include <mtl/mtl.h>

#include <mtl/utils.h>

#include <mtl/linalg_vec.h>

Sample output;

X: [1,5,3,1,]

The L-2 norm of X is 6

using namespace mtl;

using namespace std;

typedef external_vec<double> Vec;

int

main()

{

double data[] = {1,5,3,1};

Vec x(data, 4);

double s = two_norm(x);

cout << "X: ";

print_vector(x);

cout << "The L-2 norm of X is " << s << endl;

return 0;

}

****************************

error messages

****************************

Compiling...

euclid_norm.cpp

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(128) : warning C4346: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::value_type' : dependent name is not a type

prefix with 'typename' to indicate a type

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(177) : see reference to class template instantiation 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>' being compiled

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(106) : see reference to class template instantiation 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned' being compiled

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(414) : see reference to class template instantiation 'mtl::light_matrix<T,SizeType,Orien,Strided>' being compiled

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(128) : error C2144: syntax error : 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::value_type' should be preceded by ';'

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(128) : error C2208: 'T' : no members defined using this type

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(129) : warning C4346: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::pointer' : dependent name is not a type

prefix with 'typename' to indicate a type

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(129) : error C2144: syntax error : 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::pointer' should be preceded by ';'

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(129) : error C2602: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::pointer' is not a member of a base class of 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>'

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(111) : see declaration of 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::pointer'

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(129) : error C2868: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::pointer' : illegal syntax for using-declaration; expected qualified-name

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(130) : warning C4346: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::size_type' : dependent name is not a type

prefix with 'typename' to indicate a type

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(130) : error C2144: syntax error : 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::size_type' should be preceded by ';'

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(130) : error C2208: 'SizeType' : no members defined using this type

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(131) : warning C4346: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::difference_type' : dependent name is not a type

prefix with 'typename' to indicate a type

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(131) : error C2144: syntax error : 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::difference_type' should be preceded by ';'

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(131) : error C2602: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::difference_type' is not a member of a base class of 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>'

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(113) : see declaration of 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::difference_type'

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(131) : error C2868: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::difference_type' : illegal syntax for using-declaration; expected qualified-name

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(133) : warning C4346: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::const_reference' : dependent name is not a type

prefix with 'typename' to indicate a type

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(133) : error C2923: 'mtl::IF' : 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::const_reference' is not a valid template type argument for parameter 'A'

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(133) : warning C4346: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::reference' : dependent name is not a type

prefix with 'typename' to indicate a type

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(133) : error C2923: 'mtl::IF' : 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::reference' is not a valid template type argument for parameter 'B'

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(133) : error C2146: syntax error : missing ';' before identifier 'reference'

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(133) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(133) : error C2602: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::reference' is not a member of a base class of 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>'

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(108) : see declaration of 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::reference'

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(133) : error C2868: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::reference' : illegal syntax for using-declaration; expected qualified-name

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(137) : error C2061: syntax error : identifier 'size_type'

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(154) : error C2061: syntax error : identifier 'size_type'

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(154) : error C2805: binary 'operator +=' has too few parameters

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(154) : error C2333: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>::operator +=' : error in function declaration; skipping function body

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(159) : error C2061: syntax error : identifier 'size_type'

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(159) : error C2805: binary 'operator -=' has too few parameters

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(159) : error C2333: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>::operator -=' : error in function declaration; skipping function body

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(165) : error C2146: syntax error : missing ';' before identifier 'index'

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(165) : error C2433: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>::size_type' : 'inline' not permitted on data declarations

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(165) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(165) : warning C4183: 'index': missing return type; assumed to be a member function returning 'int'

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(167) : error C2143: syntax error : missing ';' before '&'

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(167) : error C2433: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>::size_type' : 'inline' not permitted on data declarations

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(167) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(167) : warning C4183: 'pos': missing return type; assumed to be a member function returning 'int'

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(168) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(168) : error C2143: syntax error : missing ';' before '&'

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(168) : error C2433: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>::size_type' : 'inline' not permitted on data declarations

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(168) : error C2086: 'oned::oned::size_type size_type' : redefinition

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(130) : see declaration of 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>::operator size_type'

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(168) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(168) : warning C4183: 'pos': missing return type; assumed to be a member function returning 'int'

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(170) : error C2146: syntax error : missing ';' before identifier 'row'

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(170) : error C2433: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>::size_type' : 'inline' not permitted on data declarations

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(170) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(170) : warning C4183: 'row': missing return type; assumed to be a member function returning 'int'

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(171) : error C2146: syntax error : missing ';' before identifier 'column'

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(171) : error C2433: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>::size_type' : 'inline' not permitted on data declarations

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(171) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(171) : warning C4183: 'column': missing return type; assumed to be a member function returning 'int'

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(174) : error C2146: syntax error : missing ';' before identifier 'i'

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(174) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(175) : error C2146: syntax error : missing ';' before identifier 'offset'

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(175) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(176) : error C2146: syntax error : missing ';' before identifier 'stride'

c:\program files\matrix template library\mtl-2.1.2-21\mtl\light_matrix.h(176) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int

Build log was saved at "file://c:\Documents and Settings\Anthony\My Documents\Visual Studio 2008\Projects\MTL_Test\MTL_Test\Debug\BuildLog.htm"

MTL_Test - 54 error(s), 11 warning(s)

========== Rebuild All: 0 succeeded, 1 failed, 0 skipped ==========

_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/

MTL: missing files

2008-11-15T02:05:52Z

Good Morning

I have just installed your library and noticed that

#include <boost/static_assert.hpp>

#include <boost/type_traits.hpp>

#include <boost/utility/enable_if.hpp>

Are missing from the distribution!

I am on a DELL machine running Windows XP.

Regards

-----------------------------------------------------------------------

Dr Paolo Remagnino

Faculty of Computing, Information Systems and Mathematics

Kingston University

voice: +44 (0)20 8547 7930

facsimile: +44 (0)20 8547 7824

e-mail: p.remagnino@...

URL: DrPaoloRemagnino.googlepages.com/home

URL: www.kingston.ac.uk/dirc

_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/

MTL: Problem: Debug assertion failed with sparse types under Microsoft Visual Studio 2005

2008-11-03T00:00:21Z

   I try to execute an example (sparse_matrix.cc) from mtl-2.1.2-22 and get "debug assertion failed" (attached). The error occur in line 85 of the example in expression B(0,2) = 2; . Then I try to execute others example. So, there is a list of break down examples: array2D.cc , sparse_copy.cc, sparse_matrix.cc, sparse_mult.cc . All this examples contain routine with sparse storage types. Error occur deep in library (5 levels of nesting), where MTL works with STL (usually call std::vector methods). And error usually occur in STL code.
   I don't know how debug this, because I must not change STL code and I really don't know how to change MTL. Only one idea in this direction: remove working vector class with another one (may be another std::vector<> implementation or boost::vector implementation). But I have not started realize yet.
   If you are familiar with this problem or suppose anything about reasons, please, tell me.

_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/

MTL_sparse_error.jpg (38K) Download Attachment

MTL: How to compile MTL under Microsoft Visual Studio 2005

2008-11-02T23:14:30Z

My previous post in the mailing list dated 2008-09-20 was about how to compile MTL for Visual C++ 7(MS VS 2003) (mtl-2.1.2-21.zip) under Visual C++ 8 (MS VS 2005). But I did not realize that it can be made easier. For easier compilation you must download MTL for Unix (mtl-2.1.2-22.tar.gz), then patch it for MS VS 2005 (mtl-2.1.2-22ForVS2005.zip). Why I didn't do like this firstly? First time when I try to compile MTL under MS VS 2005 it wasn't obvious for me that I should patch version for Unix. And I downloaded version for VS 2003 and tryed patch this version. That's why my first way of solving this problem was much more complicated.
So, I tell about two ways to compile MTL under MS VS 2005. If you find difficulties in compiling how I was, probably you should make use the easiest compile method.

_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/

MTL: Hello dear MTL-developers

2008-09-20T00:56:52Z

Good day.

My name Sergey Avdeev. And I'm a student of Novosibirsk State University (Russia). In my research work I need to manage with big matrices and vectors and I wanted it faster, then my own BLAS did this. That is why I decided to download and use your library. According to the site (http://www.osl.iu.edu/research/mtl/) MTL have many matrix/vector types and basic operators implementation. Only one problem I had, this is Microsoft Visual Studio 2005 Express Edition (I develop my project with it). I downloaded and installed mtl-2.1.2-21 (For Visual C++ 7 .NET, I know that I have VC8). So, when I try to compile the simplest code like:

#include "mtl/mtl.h"
void main()
{}

I have got a long list of compile errors in MTL's files (Error list 1.txt file in the attachment).
The first and probably the most important error was
mtl-2.1.2-21\mtl\light1d.h(87) : error C2039: '_Ptrit' : is not a member of 'std'

Then I download mtl-2.1.2-22ForVS2005.zip and patch the library with this files.
Try the same code again. And I got another long list of errors (Error list 2.txt)

I was very upset, because I did not know what to do with this error. Did not find the answer, I try to search it in the Internet. And I found something:
http://osdir.com/ml/lib.mtl.devel/2004-06/msg00004.html on this web-page one user asked, and one another answered him by posting some files, which was corrected by him to compile under .NET 2003.
List of files: mtl.h , dense2D.h , linalg_vec.h , dense_iterator.h , light_matrix.h , ligth1D.h .
I patched the library by this files and try to compile again.
That time I got much shorter error list than all previous (Error list 3.txt). I was inspired and try to eliminate errors (by my self). After light_matrix.h was a little bit changed (ligth_matrix.h in the attachment, new lines have now comment // Changed by SergeyAvd ...), then the simplest code was compiled successfully(!!!).
Next I try to compile ITL(4.0.0-1) examples with MTL. Got some errors. Change one file (gmres.h) from ITL. Only then I was able to use libraries.

I write this letter for probably improve MTL and ITL for MS VS 2005 . If my information are urgent, post it on your site to help new MS VS 2005/2008 users to start with MTL, ITL..

In the end of my letter I want to note that your library is not on top of perfection in field of documentation (manual) and have mistakes in pseudo-code (concept files). So, I can help project to improve this things.

If you are interested in my proposal answer me.

Good bye.

*****************************************************************************
**                                                                         **
** WARNING:  This email contains an attachment of a very suspicious type.  **
** You are urged NOT to open this attachment unless you are absolutely     **
** sure it is legitimate.  Opening this attachment may cause irreparable   **
** damage to your computer and your files.  If you have any questions      **
** about the validity of this message, PLEASE SEEK HELP BEFORE OPENING IT. **
**                                                                         **
** This warning was added by the IU Computer Science Dept. mail scanner.   **
*****************************************************************************

1>------ Build started: Project: TryITL2, Configuration: Debug Win32 ------
1>Compiling...
1>Main.cpp
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light1d.h(87) : error C2039: '_Ptrit' : is not a member of 'std'
1> f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light1d.h(200) : see reference to class template instantiation 'mtl::light1D<T,NN,IND_OFFSET>' being compiled
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light1d.h(87) : error C2143: syntax error : missing ';' before '<'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light1d.h(87) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light1d.h(87) : error C2238: unexpected token(s) preceding ';'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light1d.h(88) : error C2039: '_Ptrit' : is not a member of 'std'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light1d.h(88) : error C2059: syntax error : '<'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light1d.h(88) : error C2238: unexpected token(s) preceding ';'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light1d.h(90) : error C2065: 'ptr_iterator' : undeclared identifier
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light1d.h(92) : error C2065: 'ptr_const_iterator' : undeclared identifier
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light1d.h(102) : error C3203: 'dense_iterator' : unspecialized class template can't be used as a template argument for template parameter 'Iter', expected a real type
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light1d.h(104) : error C3203: 'dense_iterator' : unspecialized class template can't be used as a template argument for template parameter 'Iter', expected a real type
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\linalg_vec.h(415) : error C2039: '_Ptrit' : is not a member of 'std'
1> f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\linalg_vec.h(570) : see reference to class template instantiation 'mtl::external_vec<T,NN,SizeType>' being compiled
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\linalg_vec.h(415) : error C2143: syntax error : missing ';' before '<'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\linalg_vec.h(415) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\linalg_vec.h(415) : error C2238: unexpected token(s) preceding ';'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\linalg_vec.h(416) : error C2039: '_Ptrit' : is not a member of 'std'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\linalg_vec.h(416) : error C2059: syntax error : '<'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\linalg_vec.h(416) : error C2238: unexpected token(s) preceding ';'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\linalg_vec.h(428) : error C3203: 'dense_iterator' : unspecialized class template can't be used as a template argument for template parameter 'Iter', expected a real type
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\linalg_vec.h(430) : error C3203: 'dense_iterator' : unspecialized class template can't be used as a template argument for template parameter 'Iter', expected a real type
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(128) : warning C4346: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::value_type' : dependent name is not a type
1> prefix with 'typename' to indicate a type
1> f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(177) : see reference to class template instantiation 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>' being compiled
1> f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(106) : see reference to class template instantiation 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned' being compiled
1> f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(414) : see reference to class template instantiation 'mtl::light_matrix<T,SizeType,Orien,Strided>' being compiled
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(128) : error C2144: syntax error : 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::value_type' should be preceded by ';'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(128) : error C2208: 'T' : no members defined using this type
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(129) : warning C4346: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::pointer' : dependent name is not a type
1> prefix with 'typename' to indicate a type
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(129) : error C2144: syntax error : 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::pointer' should be preceded by ';'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(129) : error C2602: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::pointer' is not a member of a base class of 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>'
1> f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(111) : see declaration of 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::pointer'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(129) : error C2868: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::pointer' : illegal syntax for using-declaration; expected qualified-name
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(130) : warning C4346: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::size_type' : dependent name is not a type
1> prefix with 'typename' to indicate a type
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(130) : error C2144: syntax error : 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::size_type' should be preceded by ';'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(130) : error C2208: 'SizeType' : no members defined using this type
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(131) : warning C4346: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::difference_type' : dependent name is not a type
1> prefix with 'typename' to indicate a type
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(131) : error C2144: syntax error : 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::difference_type' should be preceded by ';'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(131) : error C2602: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::difference_type' is not a member of a base class of 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>'
1> f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(113) : see declaration of 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::difference_type'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(131) : error C2868: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::difference_type' : illegal syntax for using-declaration; expected qualified-name
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(133) : warning C4346: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::const_reference' : dependent name is not a type
1> prefix with 'typename' to indicate a type
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(133) : error C2923: 'mtl::IF' : 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::const_reference' is not a valid template type argument for parameter 'A'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(133) : warning C4346: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::reference' : dependent name is not a type
1> prefix with 'typename' to indicate a type
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(133) : error C2923: 'mtl::IF' : 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::reference' is not a valid template type argument for parameter 'B'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(137) : error C2061: syntax error : identifier 'size_type'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(154) : error C2061: syntax error : identifier 'size_type'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(154) : error C2805: binary 'operator +=' has too few parameters
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(154) : error C2333: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>::operator +=' : error in function declaration; skipping function body
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(159) : error C2061: syntax error : identifier 'size_type'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(159) : error C2805: binary 'operator -=' has too few parameters
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(159) : error C2333: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>::operator -=' : error in function declaration; skipping function body
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(165) : error C2146: syntax error : missing ';' before identifier 'index'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(165) : error C2433: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>::size_type' : 'inline' not permitted on data declarations
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(165) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(165) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(165) : warning C4183: 'index': missing return type; assumed to be a member function returning 'int'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(167) : error C2143: syntax error : missing ';' before '&'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(167) : error C2433: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>::size_type' : 'inline' not permitted on data declarations
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(167) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(167) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(167) : warning C4183: 'pos': missing return type; assumed to be a member function returning 'int'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(168) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(168) : error C2143: syntax error : missing ';' before '&'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(168) : error C2433: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>::size_type' : 'inline' not permitted on data declarations
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(168) : error C2086: 'oned::size_type size_type' : redefinition
1> f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(130) : see declaration of 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>::operator size_type'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(168) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(168) : warning C4183: 'pos': missing return type; assumed to be a member function returning 'int'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(170) : error C2146: syntax error : missing ';' before identifier 'row'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(170) : error C2433: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>::size_type' : 'inline' not permitted on data declarations
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(170) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(170) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(170) : warning C4183: 'row': missing return type; assumed to be a member function returning 'int'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(171) : error C2146: syntax error : missing ';' before identifier 'column'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(171) : error C2433: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>::size_type' : 'inline' not permitted on data declarations
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(171) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(171) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(171) : warning C4183: 'column': missing return type; assumed to be a member function returning 'int'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(174) : error C2146: syntax error : missing ';' before identifier 'i'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(174) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(174) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(174) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(175) : error C2146: syntax error : missing ';' before identifier 'offset'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(175) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(175) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(176) : error C2146: syntax error : missing ';' before identifier 'stride'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(176) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(176) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\dense2d.h(200) : error C2065: 'x' : undeclared identifier
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\dense2d.h(201) : error C2433: 'mtl::rect_offset<size_t,MM,NN>::{ctor}' : 'inline' not permitted on data declarations
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\dense2d.h(201) : warning C4346: 'mtl::rect_offset<size_t,MM,NN>::transpose_type' : dependent name is not a type
1> prefix with 'typename' to indicate a type
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\dense2d.h(201) : error C2350: 'mtl::rect_offset<size_t,MM,NN>::{ctor}' is not a static member
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\utils.h(477) : warning C4244: 'return' : conversion from 'double' to 'float', possible loss of data
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\mtl.h(1046) : error C2906: 'mtl::givens_rotation<std::complex<double>>' : explicit specialization requires 'template <>'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\mtl.h(1155) : error C2906: 'mtl::givens_rotation<std::complex<float>>' : explicit specialization requires 'template <>'
1>Build log was saved at "file://f:\Sergey\Work\Current projects\Try BLAS\TryITL2\TryITL2\Debug\BuildLog.htm"
1>TryITL2 - 75 error(s), 13 warning(s)
========== Build: 0 succeeded, 1 failed, 0 up-to-date, 0 skipped ==========

1>------ Build started: Project: TryITL2, Configuration: Debug Win32 ------
1>Compiling...
1>Main.cpp
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(128) : warning C4346: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::value_type' : dependent name is not a type
1> prefix with 'typename' to indicate a type
1> f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(177) : see reference to class template instantiation 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>' being compiled
1> f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(106) : see reference to class template instantiation 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned' being compiled
1> f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(414) : see reference to class template instantiation 'mtl::light_matrix<T,SizeType,Orien,Strided>' being compiled
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(128) : error C2144: syntax error : 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::value_type' should be preceded by ';'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(128) : error C2208: 'T' : no members defined using this type
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(129) : warning C4346: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::pointer' : dependent name is not a type
1> prefix with 'typename' to indicate a type
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(129) : error C2144: syntax error : 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::pointer' should be preceded by ';'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(129) : error C2602: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::pointer' is not a member of a base class of 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>'
1> f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(111) : see declaration of 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::pointer'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(129) : error C2868: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::pointer' : illegal syntax for using-declaration; expected qualified-name
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(130) : warning C4346: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::size_type' : dependent name is not a type
1> prefix with 'typename' to indicate a type
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(130) : error C2144: syntax error : 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::size_type' should be preceded by ';'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(130) : error C2208: 'SizeType' : no members defined using this type
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(131) : warning C4346: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::difference_type' : dependent name is not a type
1> prefix with 'typename' to indicate a type
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(131) : error C2144: syntax error : 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::difference_type' should be preceded by ';'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(131) : error C2602: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::difference_type' is not a member of a base class of 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>'
1> f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(113) : see declaration of 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::difference_type'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(131) : error C2868: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::difference_type' : illegal syntax for using-declaration; expected qualified-name
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(133) : warning C4346: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::const_reference' : dependent name is not a type
1> prefix with 'typename' to indicate a type
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(133) : error C2923: 'mtl::IF' : 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::const_reference' is not a valid template type argument for parameter 'A'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(133) : warning C4346: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::reference' : dependent name is not a type
1> prefix with 'typename' to indicate a type
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(133) : error C2923: 'mtl::IF' : 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::reference' is not a valid template type argument for parameter 'B'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(133) : error C2039: 'RET' : is not a member of 'mtl::IF'
1> f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\meta_if.h(31) : see declaration of 'mtl::IF'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(133) : error C2146: syntax error : missing ';' before identifier 'reference'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(133) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(133) : error C2602: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::reference' is not a member of a base class of 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>'
1> f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(108) : see declaration of 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::reference'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(133) : error C2868: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::reference' : illegal syntax for using-declaration; expected qualified-name
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(137) : error C2061: syntax error : identifier 'size_type'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(154) : error C2061: syntax error : identifier 'size_type'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(154) : error C2805: binary 'operator +=' has too few parameters
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(154) : error C2333: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>::operator +=' : error in function declaration; skipping function body
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(159) : error C2061: syntax error : identifier 'size_type'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(159) : error C2805: binary 'operator -=' has too few parameters
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(159) : error C2333: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>::operator -=' : error in function declaration; skipping function body
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(165) : error C2146: syntax error : missing ';' before identifier 'index'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(165) : error C2433: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>::size_type' : 'inline' not permitted on data declarations
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(165) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(165) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(165) : warning C4183: 'index': missing return type; assumed to be a member function returning 'int'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(167) : error C2143: syntax error : missing ';' before '&'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(167) : error C2433: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>::size_type' : 'inline' not permitted on data declarations
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(167) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(167) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(167) : warning C4183: 'pos': missing return type; assumed to be a member function returning 'int'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(168) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(168) : error C2143: syntax error : missing ';' before '&'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(168) : error C2433: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>::size_type' : 'inline' not permitted on data declarations
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(168) : error C2086: 'oned::size_type size_type' : redefinition
1> f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(130) : see declaration of 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>::operator size_type'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(168) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(168) : warning C4183: 'pos': missing return type; assumed to be a member function returning 'int'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(170) : error C2146: syntax error : missing ';' before identifier 'row'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(170) : error C2433: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>::size_type' : 'inline' not permitted on data declarations
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(170) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(170) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(170) : warning C4183: 'row': missing return type; assumed to be a member function returning 'int'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(171) : error C2146: syntax error : missing ';' before identifier 'column'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(171) : error C2433: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>::size_type' : 'inline' not permitted on data declarations
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(171) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(171) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(171) : warning C4183: 'column': missing return type; assumed to be a member function returning 'int'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(174) : error C2146: syntax error : missing ';' before identifier 'i'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(174) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(174) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(174) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(175) : error C2146: syntax error : missing ';' before identifier 'offset'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(175) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(175) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(176) : error C2146: syntax error : missing ';' before identifier 'stride'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(176) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(176) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int
1>Build log was saved at "file://f:\Sergey\Work\Current projects\Try BLAS\TryITL2\TryITL2\Debug\BuildLog.htm"
1>TryITL2 - 55 error(s), 11 warning(s)
========== Build: 0 succeeded, 1 failed, 0 up-to-date, 0 skipped ==========

1>------ Build started: Project: TryITL2, Configuration: Debug Win32 ------
1>Compiling...
1>Main.cpp
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(129) : warning C4346: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::pointer' : dependent name is not a type
1> prefix with 'typename' to indicate a type
1> f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(177) : see reference to class template instantiation 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>' being compiled
1> f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(106) : see reference to class template instantiation 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned' being compiled
1> f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(414) : see reference to class template instantiation 'mtl::light_matrix<T,SizeType,Orien,Strided>' being compiled
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(129) : error C2144: syntax error : 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::pointer' should be preceded by ';'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(129) : error C2602: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::pointer' is not a member of a base class of 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>'
1> f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(111) : see declaration of 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::pointer'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(129) : error C2868: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::pointer' : illegal syntax for using-declaration; expected qualified-name
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(131) : warning C4346: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::difference_type' : dependent name is not a type
1> prefix with 'typename' to indicate a type
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(131) : error C2144: syntax error : 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::difference_type' should be preceded by ';'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(131) : error C2602: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::difference_type' is not a member of a base class of 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>'
1> f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(113) : see declaration of 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::difference_type'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(131) : error C2868: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::difference_type' : illegal syntax for using-declaration; expected qualified-name
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(133) : warning C4346: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::const_reference' : dependent name is not a type
1> prefix with 'typename' to indicate a type
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(133) : error C2923: 'mtl::IF' : 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::const_reference' is not a valid template type argument for parameter 'A'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(133) : warning C4346: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::reference' : dependent name is not a type
1> prefix with 'typename' to indicate a type
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(133) : error C2923: 'mtl::IF' : 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::reference' is not a valid template type argument for parameter 'B'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(133) : error C2039: 'RET' : is not a member of 'mtl::IF'
1> f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\meta_if.h(31) : see declaration of 'mtl::IF'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(133) : error C2146: syntax error : missing ';' before identifier 'reference'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(133) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(133) : error C2602: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::reference' is not a member of a base class of 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::__iterator<isConst>'
1> f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(108) : see declaration of 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::reference'
1>f:\sergey\work\current projects\try blas\iterative template library\mtl-2.1.2-21\mtl\light_matrix.h(133) : error C2868: 'mtl::light_matrix<T,SizeType,Orien,Strided>::oned::reference' : illegal syntax for using-declaration; expected qualified-name
1>Build log was saved at "file://f:\Sergey\Work\Current projects\Try BLAS\TryITL2\TryITL2\Debug\BuildLog.htm"
1>TryITL2 - 13 error(s), 4 warning(s)
========== Build: 0 succeeded, 1 failed, 0 up-to-date, 0 skipped ==========

#ifndef MTL_LIGHT_MATRIX_H
#define MTL_LIGHT_MATRIX_H

#include "mtl/matrix_traits.h"
#include "mtl/dimension.h"
#include "mtl/meta_if.h"
#include "mtl/meta_equal.h"

namespace mtl {

template <int Orien>
struct TRANS {
enum { RET = 0 };
};

template<>
struct TRANS<ROW_MAJOR> {
enum { RET = COL_MAJOR };
};

template<>
struct TRANS<COL_MAJOR> {
enum { RET = ROW_MAJOR };
};

template <class T, class SizeType, int Orien, int Strided>
class light_matrix {
public:
typedef light_matrix self;
typedef light_matrix light_matrix_t; // VC++ workaround
typedef T* DataPtr;

typedef rectangle_tag shape;
typedef typename IF< EQUAL<Orien,ROW_MAJOR>::RET,
row_tag, column_tag>::RET orientation; // mostly wrong

typedef typename IF< EQUAL<Orien,ROW_MAJOR>::RET,
row_orien, column_orien>::RET orien;

typedef light_matrix<T, SizeType, TRANS<Orien>::RET, Strided> transpose_type;
typedef light_matrix<T, SizeType, Orien, !Strided> strided_type;
typedef light_matrix<T, SizeType, Orien, Strided> scaled_type;// wrong

typedef light_matrix<T, SizeType, Orien, Strided> submatrix_type;

typedef int DiffType;

//: The size type
typedef SizeType size_type;
//: The type for differences between iterators
typedef DiffType difference_type;

typedef T value_type;
typedef value_type& reference;
typedef const value_type& const_reference;
typedef value_type* pointer;

enum { M = 0, N = 0 };

protected:

static inline size_type& twod_pos(size_type& i, size_type& j) {
if (Orien == ROW_MAJOR)
return i;
else
return j;
}

static inline const size_type& twod_pos(const size_type& i,
const size_type& j) {
if (Orien == ROW_MAJOR)
return i;
else
return j;
}

static inline size_type& oned_pos(size_type& i, size_type& j) {
if (Orien == ROW_MAJOR)
return j;
else
return i;
}

static inline const size_type& oned_pos(const size_type& i,
const size_type& j) {
if (Orien == ROW_MAJOR)
return j;
else
return i;
}

// idea: completely separate stride/offset/positioning from indexing
// but encapsulate both somehow

public:

//: This is a dense 2D container
typedef dense_tag sparsity;
//: This has external storage
typedef external_tag storage_loc;
//: This is strideable
typedef strideable strideability;

class oned {
public:
typedef T& reference;
typedef const T& const_reference;
typedef T value_type;
typedef T* pointer;
typedef SizeType size_type;
typedef int difference_type;

enum { M = 0, N = 0 };

typedef oned subrange_type;
typedef dense_tag sparsity;
typedef oned IndexArray; /* bogus */
typedef oned IndexArrayRef; /* bogus */

typedef oned_tag dimension; /* bogus */

template <int isConst>
class __iterator {
typedef __iterator self;
public:
typedef typename oned::value_type value_type;
typedef typename oned::pointer pointer; // changed by SergeyAvd, old version: typedef oned::pointer pointer;
typedef typename oned::size_type size_type;
typedef typename oned::difference_type difference_type; // changed by SergeyAvd, old version: typedef oned::difference_type difference_type;

typedef typename IF<isConst, typename oned::const_reference, typename oned::reference>::RET reference; // changed by SergeyAvd, old version: typedef typename IF<isConst, oned::const_reference, oned::reference>::RET reference;

typedef std::random_access_iterator_tag iterator_category;

inline __iterator(DataPtr d,
size_type ii, size_type jj,
size_type os, size_type s)
: data(d), i(ii), j(jj), offset(os), stride(s) { }

inline __iterator(const self& x)
: data(x.data), i(x.i), j(x.j), offset(x.offset), stride(x.stride) { }

inline self& operator=(const self& x) {
data = x.data; i = x.i; j = x.j; offset = x.offset; stride = x.stride;
return *this;
}

inline __iterator() : data(0), i(0), j(0), offset(0), stride(0) { }

inline reference operator*() const { return data[offset]; }
inline self& operator++() { ++pos(); offset += stride; return *this; }
inline self& operator+=(size_type n) {
pos() += n; offset += stride*n; return *this;
}
inline self operator++(int) { self t = *this; ++(*this); return t; }
inline self& operator--() { --pos(); offset -= stride; return *this; }
inline self& operator-=(size_type n) {
pos() -= n; offset -= stride*n; return *this; }
inline self operator--(int) { self t = *this; --(*this); return t; }
inline bool operator!=(const self& x) const { return pos() != x.pos(); }
inline bool operator==(const self& x) const { return pos() == x.pos(); }
inline bool operator<(const self& x) const { return pos() < x.pos(); }
inline size_type index() const { return pos(); }

inline size_type& pos() { return oned_pos(i,j); }
inline const size_type& pos() const { return oned_pos(i,j); }

inline size_type row() const { return i; }
inline size_type column() const { return j; }
protected:
DataPtr data;
size_type i, j;
size_type offset;
size_type stride;
};

typedef __iterator<0> iterator;
typedef __iterator<1> const_iterator;

inline oned(DataPtr d, size_type ii, size_type jj,
size_type ie, size_type je,
size_type os, size_type ld)
: data(d), i(ii), j(jj), iend(ie), jend(je),
offset(os), ldim(ld) { }

inline oned(const oned& x)
: data(x.data), i(x.i), j(x.j),
iend(x.iend), jend(x.jend),
offset(x.offset), ldim(x.ldim) { }

inline oned& operator=(const oned& x) {
data = x.data; i = x.i; j = x.j;
iend = x.iend; jend = x.jend;
offset = x.offset; ldim = x.ldim;
return *this;
}
inline oned()
: data(0), i(0), j(0), iend(0), jend(0), offset(0), ldim(0) { }

inline ~oned() { }

inline reference operator[](size_type n) {
return data[ Strided ? offset + n * ldim : offset + n];
}

inline const_reference operator[](size_type n) const {
return data[ Strided ? offset + n * ldim : offset + n];
}

inline iterator begin() {
return iterator(data, i, j, offset, Strided ? ldim : 1);
}
inline iterator end() {
size_type iiend, jjend;
if (Orien == ROW_MAJOR) { iiend = i; jjend = jend; }
else { iiend = iend; jjend = j; }

return iterator(data, iiend, jjend, offset, Strided ? ldim: 1);
}

inline const_iterator begin() const {
return const_iterator(data, i, j, offset, Strided ? ldim : 1);
}
inline const_iterator end() const {
size_type iiend, jjend;
if (Orien == ROW_MAJOR) { iiend = i; jjend = jend; }
else { iiend = iend; jjend = j; }

return const_iterator(data, iiend, jjend, offset, Strided ? ldim : 1);
}

protected:
DataPtr data;
size_type i, j;
size_type iend, jend;
size_type offset;
size_type ldim;
};

typedef oned OneD;
typedef OneD OneDRef;
typedef OneD Row;
typedef OneD RowRef;
typedef OneD Column;
typedef OneD ColumnRef;

//: The iterator type
template <int Const>
class __iterator {
typedef __iterator self;
public:
typedef std::random_access_iterator_tag iterator_category;
typedef oned value_type;
typedef value_type* pointer;

#if defined(_MSVCPP_)
typedef typename light_matrix_t::size_type size_type;
typedef typename light_matrix_t::difference_type difference_type;
#else
typedef SizeType size_type;
typedef DiffType difference_type;
#endif

typedef typename IF<Const, const oned, oned>::RET reference;

inline __iterator(DataPtr d, size_type ii, size_type jj,
size_type ie, size_type je, size_type ld)
: data(d), i(ii), j(jj), iend(ie), jend(je), offset(0), ldim(ld) {
if (Strided) stride = 1; else stride = ldim;
}

inline __iterator() : data(0), i(0), j(0), iend(0), jend(0),
offset(0), stride(0), ldim(0) { }

inline __iterator(const self& x)
: data(x.data), i(x.i), j(x.j),
iend(x.iend), jend(x.jend), offset(x.offset),
stride(x.stride), ldim(x.ldim) { }

inline self& operator=(const self& x) {
data = x.data; i = x.i; j = x.j;
iend = x.iend; jend = x.jend;
offset = x.offset; stride = x.stride; ldim = x.ldim;
return *this;
}
inline reference operator*() const {
return oned(data, i, j, iend, jend, offset, ldim);
}

inline self& operator++() { ++pos(); offset += stride; return *this; }
inline self& operator+=(size_type n) {
pos() += n; offset += stride*n; return *this; }
inline self operator++(int) { self t = *this; ++(*this); return t; }
inline self& operator--() { --pos(); offset -= stride; return *this; }
inline self& operator-=(size_type n) {
pos() -= n; offset -= stride*n; return *this; }
inline self operator--(int) { self t = *this; --(*this); return t; }
inline bool operator!=(const self& x) const { return pos() != x.pos(); }
inline bool operator==(const self& x) const { return pos() == x.pos(); }
inline bool operator<(const self& x) const { return pos() < x.pos(); }
inline size_type index() const { return pos(); }

inline size_type& pos() { return twod_pos(i,j); }
inline const size_type& pos() const { return twod_pos(i,j); }

inline size_type row() const { return i; }
inline size_type column() const { return j; }

protected:
DataPtr data;
size_type i, j;
size_type iend, jend;
size_type offset;
size_type stride;
size_type ldim;
};

typedef __iterator<0> iterator;
typedef __iterator<1> const_iterator;

//: Standard Constructor
inline light_matrix(DataPtr d, size_type m, size_type n, size_type ld)
: data_(d), nrows_(m), ncols_(n), ldim(ld) { }

inline light_matrix(DataPtr d, size_type m, size_type n)
: data_(d), nrows_(m), ncols_(n), ldim(Orien == ROW_MAJOR ? n : m) { }

//: Copy Constructor
inline light_matrix(const light_matrix& x)
: data_(x.data_), nrows_(x.nrows_), ncols_(x.ncols_), ldim(x.ldim) { }

//: Assignment Operator
inline const light_matrix& operator=(const light_matrix& x) {
data_ = x.data_; nrows_ = x.nrows_; ncols_ = x.ncols_; ldim = x.ldim;
return *this;
}
//: Default Constructor
inline light_matrix() : data_(0), nrows_(0), ncols_(0), ldim(0) { }

inline light_matrix(const strided_type& x, do_strided s)
: data_(x.data_), nrows_(x.nrows_), ncols_(x.ncols_), ldim(x.ldim) { }

template <class StridedType>
inline light_matrix(const StridedType& x, do_strided s)
: data_(x.data_), nrows_(x.nrows_), ncols_(x.ncols_), ldim(x.ldim) { }

//: Destructor
inline ~light_matrix() { }

//: Return an iterator pointing to the first 1D container
inline iterator begin() {
return iterator(data_, 0, 0, nrows_, ncols_, ldim);
}
//: Return an iterator pointing past the end of the 2D container
inline iterator end() {
return iterator(data_, nrows_, ncols_, nrows_, ncols_, ldim);
}

//: Return a const iterator pointing to the first 1D container
inline const_iterator begin() const {
return const_iterator(data_, 0, 0, nrows_, ncols_, ldim);
}
//: Return a const iterator pointing past the end of the 2D container
inline const_iterator end() const {
return const_iterator(data_, nrows_, ncols_, nrows_, ncols_, ldim);
}

//: Return a reference to the ith 1D container
inline oned operator[](size_type n) {
if (Orien == ROW_MAJOR)
return oned(data_, n, 0, nrows_, ncols_, Strided ? n : ldim * n, ldim);
else
return oned(data_, 0, n, nrows_, ncols_, Strided ? n : ldim * n, ldim);
}

inline const oned operator[](size_type n) const {
if (Orien == ROW_MAJOR)
return oned(data_, n, 0, nrows_, ncols_, Strided ? n : ldim * n, ldim);
else
return oned(data_, 0, n, nrows_, ncols_, Strided ? n : ldim * n, ldim);
}

//: Return a reference to the (i,j) element, where (i,j) is in the 2D coordinate system
inline reference operator()(size_type i, size_type j) {
return Orien == ROW_MAJOR ? operator[](i)[j] : operator[](j)[i];
}

//: Return a const reference to the (i,j) element, where (i,j) is in the 2D coordinate system
inline const_reference operator()(size_type i, size_type j) const {
return Orien == ROW_MAJOR ? operator[](i)[j] : operator[](j)[i];
}

inline size_type nrows() const { return nrows_; }
inline size_type ncols() const { return ncols_; }

inline submatrix_type sub_matrix(size_type i, size_type iend,
size_type j, size_type jend) const
{
if (Strided)
return submatrix_type(data_ + oned_pos(i,j) * ldim + twod_pos(iend,jend),
iend - i, jend - j, ldim);
else
return submatrix_type(data_ + twod_pos(i,j) * ldim + oned_pos(iend,jend),
iend - i, jend - j, ldim);
}

DataPtr data_;
size_type nrows_, ncols_;
size_type ldim;
};

} /* namespace mtl */

#endif /* MTL_LIGHT_MATRIX_H */

// -*- c++ -*-
// $COPYRIGHT$
//

#ifndef ITL_KRYLOV_GMRES_H
#define ITL_KRYLOV_GMRES_H

#include <vector>
#include <algorithm>
#include "itl/itl.h"
#include "itl/givens_rotation.h"
#include "itl/number_traits.h"
#include <itl/modified_gram_schmidt.h>

namespace itl {

// Generalized Minimum Residual
//
// This solve the unsymmetric linear system Ax = b using restarted GMRES.
//
// A return value of 0 indicates convergence within the
// maximum number of iterations (determined by the iter object).
// A return value of 1 indicates a failure to converge.
//
// On instantiating Iteration object outer, the first parameter
// Vector w must be the precondtioned one, i.e., solve(M, b, w) where
// b is right side of linear system. See test_gmres.cc for example.
//
// See: Y. Saad and M. Schulter. GMRES: A generalized minimum residual
// algorithm for solving nonsysmmetric linear systems, SIAM
// J. Sci. Statist. Comp. 7(1986), pp, 856-869
//
/* required operations: mult,copy,dot_conj,add,scaled,two_norm,tri_solve */
template < class Matrix, class Vector, class VectorB, class Preconditioner,
class Iter, class Basis >
int
gmres(const Matrix &A, Vector &x, const VectorB &b,
const Preconditioner &M, int restart, Iter& outer, Basis& KS)
{
typedef typename itl_traits<Vector>::value_type T;
typedef typename itl_traits<Vector>::size_type size_type;

typedef mtl::dense1D<T> TmpVec; // Changed by SergeyAvd, old version: typedef std::vector<T> TmpVec;
typedef Vector internal_vector;

//These must be arithmetically compatible with x
internal_vector w(size(x)), r(size(x)), u(size(x));

typedef typename number_traits<T>::magnitude_type Real;
typedef typename internal_matrix_traits<T>::Matrix HMat;
HMat H(restart+1, restart); //Elements in H must be real type

TmpVec s(restart+1);

std::vector< itl::givens_rotation<T> > rotations(restart+1);

itl::mult(A, itl::scaled(x, -1.0), b, w);

itl::solve(M, w, r);
Real beta = std::abs(itl::two_norm(r));

while (! outer.finished(beta)) {

itl::copy(itl::scaled(r, 1./beta), KS[0]);
std::fill(s.begin(), s.end(), 0.0);
s[0] = beta;

size_type i = 0;
Iter inner(outer.normb(), restart, outer.tol(), outer.atol());

do {
size_type k;
itl::mult(A, KS[i], u);
itl::solve(M, u, KS[i+1]);

itl::orthogonalize(KS, H[i], i);

Real H_ip1_i = itl::two_norm(KS[i+1]);
H(i+1, i) = H_ip1_i;
itl::scale(KS[i+1], 1./H_ip1_i);

for (k = 0; k < i; k++)
rotations[k].scalar_apply(H(k,i), H(k+1,i));

rotations[i] = itl::givens_rotation<T>(H(i,i), H(i+1,i));
rotations[i].scalar_apply(H(i,i), H(i+1,i));
rotations[i].scalar_apply(s[i], s[i+1]);

++inner, ++outer, ++i;

} while (! inner.finished(std::abs(s[i])));

itl::upper_tri_solve(H, s, i);

itl::combine(KS, s, x, i);

itl::mult(A, itl::scaled(x, -1.0), b, w);
itl::solve(M, w, r);
beta = std::abs(itl::two_norm(r));
}

return outer.error_code();
}

template < class Matrix, class Vector, class VectorB, class Preconditioner,
class Iter >
int
gmres(const Matrix &A, Vector &x, const VectorB &b,
const Preconditioner &M, int restart, Iter& outer)
{
itl::modified_gram_schmidt<Vector> orth(restart, size(x));
return itl::gmres(A, x, b, M, restart, outer, orth);
}

}

#endif

_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/

Aaron patch.zip (49K) Download Attachment

MTL: installing and using MTL failed.

2008-08-10T03:22:55Z

Dear,

I am trying to install MTL from tar ball mtl-2.1.2-22.tar.gz <http://www.osl.iu.edu/download/research/mtl//mtl-2.1.2-22.tar.gz> given on the webpage. However, I recently had some problems on autoconf and automake and after some search found out that, some version inconsistency was the reason for this. I change to autoconf 2.13 and automake 1.4-p6. I applied also the patch given on the mtl installation and download page, however that patch is already for gcc version 4.02 and I have been using gcc 4.2.3 on ubuntu linux. Patch application succeded however(I think this was the 1st time I applied a patch. I am glad I learnt how to do that besides.) And following the order given on the download and installation page of the mtl, I was at the configure step, I read from the INSTALL file, some extra options have to supplied to work with LAPACK and so on.

I used ./configure --with-blais --with-lapack="-llapack -lblas -lg2c"

which went ok except
.
.
.
checking whether to use interface to LAPACK... yes
test: 2807: -lblas: unexpected operator
checking whether the Fortran 77 compiler (g77 ) works... yes
.
.
.
The I tried make and I got this errors

Making all in test
make[1]: Entering directory `/home/utabak/Desktop/mtl-2.1.2-22/test'
Making all in src
make[2]: Entering directory `/home/utabak/Desktop/mtl-2.1.2-22/test/src'
g++ -DHAVE_CONFIG_H -I. -I. -I../../mtl -I../.. -ftemplate-depth-30 -Wall -c matrix_iterator_test.cc
In file included from ../../mtl/matrix_implementation.h:32,
from ../../mtl/matrix.h:33,
from matrix_iterator_test.cc:15:
../../mtl/initialize.h: In function 'void mtl::__initialize(Matrix&, mtl::matrix_market_stream<T>&, mtl::symmetric_tag)':
../../mtl/initialize.h:49: error: there are no arguments to 'assert' that depend on a template parameter, so a declaration of 'assert' must be available
../../mtl/initialize.h:49: error: (if you use '-fpermissive', G++ will accept your code, but allowing the use of an undeclared name is deprecated)
make[2]: *** [matrix_iterator_test.o] Error 1
make[2]: Leaving directory `/home/utabak/Desktop/mtl-2.1.2-22/test/src'
make[1]: *** [all-recursive] Error 1
make[1]: Leaving directory `/home/utabak/Desktop/mtl-2.1.2-22/test'
make: *** [all-recursive] Error 1

Could somebody direct me from this point on to find the point where I am mistaken. Do I strictly need gcc 4.02? Or there is another workaround to this problem...

Thanks for your help and time in advance.

Umut

_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/

MTL: diag_matrix bug

2008-05-23T11:42:59Z

Hello,

I think I've found a bug when multiplying a diagonal Matrix times a Vector.

It appears that mult_generic__ was being called, picking up the banded_tag in the diagonal matrix.

Attached is a test program which produces the following output:

V=
[6,7,8,9,10,]
Diag A=
5x5
[
[1,0,0,0,0],
[0,2,0,0,0],
[0,0,3,0,0],
[0,0,0,4,0],
[0,0,0,0,5]
]
Full B=
5x5
[
[1,0,0,0,0],
[0,2,0,0,0],
[0,0,3,0,0],
[0,0,0,4,0],
[0,0,0,0,5]
]
Output AxV=
[6,12,18,24,30,]
Output BxV=
[6,14,24,36,50,]

Also attached is a patch that adds a diag_mult__ method to mtl.h. I have to admit: templates as well as matrices are a little out of my league, but the patch solved my problem!

Great work on the project! Hope this helps.

Thanks,

Aaron

_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/

mtl_bug.cpp (1K) Download Attachment
mtl_diag.patch (1K) Download Attachment

MTL: Matlab equivalent functions: rref and null

2008-03-15T12:32:42Z

Hi,
I would like to convert some matlab code in C++. It uses rref and null functions all over the place. What are the equivalent functions for these in mtl?

--
S M Mahbub Murshed
_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/

MTL: BIT Mask and change of Matrix order

2008-01-10T04:56:44Z

Hello all,
I am resending my message to this list after I have send it to the
MTL4 list in case someone can help here.
I am working on Peano order matrix recursive multiplication traversal.
In the paper "representation-transparent Matrix Algorithms with
scalable performance" by Gottschling et al., it is mentioned that
using MTL4 it is easy to choose the order by changing the Bit-Mask
representation of the matrix.
I would like to ask if anyone knows or have experienced on changing
the order using the bit-masks and whether this can be done through
passing some parameters to the matrix type, or it requires some work
on the predefined bit-masks in the library header files
(predefinedmasks.hpp, recursion.hpp, bit-masking.hpp(generatemask
class) ).

thanks.

----------------------------------------------------------------

_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/

Re: MTL: Loop-free submatrix extraction

2007-12-27T04:18:36Z

Hi Michael,

Thank you for reminding me on that. I thought I already fixed that. Early January I hope finding the time for that.

Cheers,

Peter

On Dec 26, 2007, at 8:44 PM, Michael Smolsky wrote:

Peter,

I seem to be unable to say

result += op1 * transposed_view <const M> (op2),

where op1, op2, result are of the type dense2D<double> and M is a typedef for dense2D<double>.

The compiler barks at the absense of a type named type in struct mtl::traits::mult_result<compressed2D<...>, mtl::transposed_view<compressed2D<..> > >.

The code, that tests a product of a dense and a transposed matrix seems to be commented out from matrix_product_test.cpp.

Does MTL support calculation of the product of a transposed and dense matrix?

Thank you,

Michael

--
Got No Time? Shop Online for Great Gift Ideas!
mail.com Shopping
_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/

------------

Peter Gottschling

Research Associate

Open Systems Laboratory

Indiana University

135 Lindley Hall

Bloomington, IN 47405

Tel.: +1-812-855-3608 Fax: +1-812-856-0853

http://www.osl.iu.edu/~pgottsch

_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/

Re: MTL: Loop-free submatrix extraction

2007-12-26T11:44:53Z

Peter,

I seem to be unable to say

result += op1 * transposed_view <const M> (op2),

where op1, op2, result are of the type dense2D<double> and M is a typedef for dense2D<double>.

The compiler barks at the absense of a type named type in struct mtl::traits::mult_result<compressed2D<...>, mtl::transposed_view<compressed2D<..> > >.

The code, that tests a product of a dense and a transposed matrix seems to be commented out from matrix_product_test.cpp.

Does MTL support calculation of the product of a transposed and dense matrix?

Thank you,

Michael

--

Got No Time? Shop Online for Great Gift Ideas!
mail.com Shopping

_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/

Re: MTL: creation of matrix elements

2007-12-23T04:01:23Z

On Sunday 23 December 2007, Michael Smolsky wrote:
> Thank you for your insight, Peter.
>
> By looking at your code (contiguous_memory_block.hpp,
> generic_array<..>::alloc(..)) I've noticed, that you don't initialize
> matrix elements on allocation. That method seems to be called from some
> methods of dense2D, such as resize(..). If that is the case, the library
> won't work with any class, that has a non-trivial default constructor.
>

If you look at the current boost::ublas, you'll find e.g.:

unbounded_array (size_type size, const ALLOC &a = ALLOC()):
alloc_(a), size_ (size) {
if (size_) {
data_ = alloc_.allocate (size_);
if (! detail::has_trivial_constructor<T>::value) {
for (pointer d = data_; d != data_ + size_; ++d)
alloc_.construct(d, value_type());

also:
namespace detail {

// specialisation which define whether a type has a trivial
constructor
// or not. This is used by array types.
template<typename T>
struct has_trivial_constructor : public
boost::has_trivial_constructor<T> {};

template<typename T>
struct has_trivial_destructor : public
boost::has_trivial_destructor<T> {};

template<typename FLT>
struct has_trivial_constructor<std::complex<FLT> > : public
boost::true_type {};

template<typename FLT>
struct has_trivial_destructor<std::complex<FLT> > : public
boost::true_type {};

}

I believe this is a good approach
_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/

MTL: creation of matrix elements

2007-12-22T22:36:19Z

Thank you for your insight, Peter.

By looking at your code (contiguous_memory_block.hpp, generic_array<..>::alloc(..)) I've noticed, that you don't initialize matrix elements on allocation. That method seems to be called from some methods of dense2D, such as resize(..). If that is the case, the library won't work with any class, that has a non-trivial default constructor.

Is this by design?

In order to cover classes with non-trivial default constructors (STL has no problem with this), you might want to invoke in-place operator new on the newly allocated memory in generic_array<..>::alloc(..) and perhaps other places as well. If the default constructor of the element type is trivial, the compiler will probably optimize this call out and you won't loose in performance.

Michael.

----- Original Message -----
From: "Peter Gottschling"
To: "General Matrix Template Library $MTL$ list"
Subject: Re: MTL: Loop-free submatrix extraction
Date: Sat, 22 Dec 2007 17:24:03 +0100

Hi again, Michael,

For dense2D and morton_dense there is a function sub_matrix (shame on me that there is no documentation yet). You can use it like:
b= sub_matrix(a, begin_row, end_row, begin_column, end_column);

As you see it only supports intervals of rows and columns not arbitrary sets. The reason is that this can be realized efficiently with a view.

If you are going to write a function I would do something like (only sketch not tested code):

template <typename Matrix, typename RowVector, typename ColVector>
Matrix list_sub_matrix(const Matrix& matrix, const RowVector& row_vector, const ColVector& col_vector)
{
  Matrix s(size(row_vector), size(col_vector));
  {
   matrix::inserter<Matrix> ins(s);
   for (size_t r= 0; r < size(row_vector); r++)
   for (size_t c= 0; c < size(col_vector); c++)
   ins(r, c) << matrix[row_vector[r]][col_vector[c]];
  }
  return s;
}

Returning a matrix is quite expensive at the moment but will be okay once I've implemented move semantics. I hope I haven't overseen a detail otherwise you can asked again. Once you have it done and tested I'd like to put it in the repository.

Good luck and Merry Christmas,
Peter

Want an e-mail address like mine?
Get a free e-mail account today at www.mail.com!

_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/

Re: MTL: Loop-free submatrix extraction

2007-12-22T08:24:03Z

Hi again, Michael,

For dense2D and morton_dense there is a function sub_matrix (shame on me that there is no documentation yet). You can use it like:

b= sub_matrix(a, begin_row, end_row, begin_column, end_column);

As you see it only supports intervals of rows and columns not arbitrary sets. The reason is that this can be realized efficiently with a view.

If you are going to write a function I would do something like (only sketch not tested code):

template <typename Matrix, typename RowVector, typename ColVector>

Matrix list_sub_matrix(const Matrix& matrix, const RowVector& row_vector, const ColVector& col_vector)

{

Matrix s(size(row_vector), size(col_vector));

{

matrix::inserter<Matrix> ins(s);

for (size_t r= 0; r < size(row_vector); r++)

for (size_t c= 0; c < size(col_vector); c++)

ins(r, c) << matrix[row_vector[r]][col_vector[c]];

}

return s;

}

Returning a matrix is quite expensive at the moment but will be okay once I've implemented move semantics. I hope I haven't overseen a detail otherwise you can asked again. Once you have it done and tested I'd like to put it in the repository.

Good luck and Merry Christmas,

Peter

On Dec 22, 2007, at 8:24 AM, Michael Smolsky wrote:

Peter,

Thank you for your detailed response.

I'm exploring MTL further. I hope you don't mind if I ask you some questions, the answers to which I couldn't find in the documentation.

I need to extract a submatrix of a matrix, in a way similar to:

octave:4> a=rand(4),b=a([1,2],[2,4])
a =

   0.622575 0.844744 0.424293 0.207615
   0.642011 0.102342 0.011919 0.168107
   0.774028 0.446422 0.076552 0.759136
   0.633285 0.792054 0.738877 0.982645

b =

   0.84474 0.20762
   0.10234 0.16811

Is there an elegant way of doing this, or I need to write loops? I feel optimistic having read about inserters in your documentation...

My matrix is dense for now, but I will probably need to apply this to a sparse matrix as well.

Thanks you in advance,

Michael.

----- Original Message -----
From: "Peter Gottschling" <pgottsch@...>
To: "General Matrix Template Library (MTL) list" <mtl-devel@...>, "Matrix Template Library 4" <mtl4@...>
Subject: Re: MTL: MTL development road
Date: Mon, 03 Dec 2007 10:53:38 +0100

--
Want an e-mail address like mine?
Get a free e-mail account today at www.mail.com!

_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/

------------

Peter Gottschling

Research Associate

Open Systems Laboratory

Indiana University

135 Lindley Hall

Bloomington, IN 47405

Tel.: +1-812-855-3608 Fax: +1-812-856-0853

http://www.osl.iu.edu/~pgottsch

_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/

Re: MTL: Patch for mtl4-alpha-1-r6398, gcc 4.1.2

2007-12-22T07:58:45Z

Hi Michael,

Thank you for the patch. I'll apply it eventually.

I'm kind of surprised that a C++ compiler complains about this, esp. as I already tried gcc 4.1.2 on. It's also surprising that I put semicolons after function definitions so often.

Cheers,

Peter

On Dec 22, 2007, at 6:13 AM, Michael Smolsky wrote:

Hi,

I'm trying to compile mtl4-alpha-1-r6398.tar.gz using gcc 4.1.2
20070925 (Red Hat 4.1.2-27), and it fails on a bunch of files. Essentially, there're extra semicolons after definitions of some functions, which gcc doesn't like any more.

I've created a patch file, that I'm attaching. To apply the patch, go to mtl4 directory and execute

patch -p1 < /path/to/patch/file.txt

Thank you for supporting MTL,

Michael

--
Want an e-mail address like mine?
Get a free e-mail account today at www.mail.com!

<patch-alpha-1-r6398-remove-semicolons.txt>_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/

------------

Peter Gottschling

Research Associate

Open Systems Laboratory

Indiana University

135 Lindley Hall

Bloomington, IN 47405

Tel.: +1-812-855-3608 Fax: +1-812-856-0853

http://www.osl.iu.edu/~pgottsch

_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/

MTL: Loop-free submatrix extraction

2007-12-21T23:24:05Z

Peter,

Thank you for your detailed response.

I'm exploring MTL further. I hope you don't mind if I ask you some questions, the answers to which I couldn't find in the documentation.

I need to extract a submatrix of a matrix, in a way similar to:

octave:4> a=rand(4),b=a([1,2],[2,4])
a =

0.622575 0.844744 0.424293 0.207615
0.642011 0.102342 0.011919 0.168107
0.774028 0.446422 0.076552 0.759136
0.633285 0.792054 0.738877 0.982645

b =

0.84474 0.20762
0.10234 0.16811

Is there an elegant way of doing this, or I need to write loops? I feel optimistic having read about inserters in your documentation...

My matrix is dense for now, but I will probably need to apply this to a sparse matrix as well.

Thanks you in advance,

Michael.

> ----- Original Message -----
> From: "Peter Gottschling" <pgottsch@...>
> To: "General Matrix Template Library (MTL) list" <mtl-devel@...>, "Matrix Template Library 4" <mtl4@...>
> Subject: Re: MTL: MTL development road
> Date: Mon, 03 Dec 2007 10:53:38 +0100
>
>

--
Want an e-mail address like mine?
Get a free e-mail account today at www.mail.com!

_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/

MTL: Patch for mtl4-alpha-1-r6398, gcc 4.1.2

2007-12-21T21:13:54Z

Hi,

I'm trying to compile mtl4-alpha-1-r6398.tar.gz using gcc 4.1.2
20070925 (Red Hat 4.1.2-27), and it fails on a bunch of files. Essentially, there're extra semicolons after definitions of some functions, which gcc doesn't like any more.

I've created a patch file, that I'm attaching. To apply the patch, go to mtl4 directory and execute

patch -p1 < /path/to/patch/file.txt

Thank you for supporting MTL,

Michael

--
Want an e-mail address like mine?
Get a free e-mail account today at www.mail.com!

diff -ru mtl4/boost/numeric/mtl/operation/conj.hpp mtl4-fixed/boost/numeric/mtl/operation/conj.hpp
--- mtl4/boost/numeric/mtl/operation/conj.hpp 2007-10-26 05:33:04.000000000 -0400
+++ mtl4-fixed/boost/numeric/mtl/operation/conj.hpp 2007-12-21 23:50:35.000000000 -0500
@@ -97,7 +97,7 @@
inline conj(const Value& v)
{
return sfunctor::conj<Value, typename traits::algebraic_category<Value>::type>::apply(v);
-};
+}

namespace sfunctor {
@@ -129,7 +129,7 @@
inline typename sfunctor::real<Value>::result_type real(const Value& v)
{
return sfunctor::real<Value>::apply(v);
-};
+}

namespace sfunctor {
@@ -163,7 +163,7 @@
inline typename sfunctor::imag<Value>::result_type imag(const Value& v)
{
return sfunctor::imag<Value>::apply(v);
-};
+}

} // namespace mtl
diff -ru mtl4/boost/numeric/mtl/operation/trans.hpp mtl4-fixed/boost/numeric/mtl/operation/trans.hpp
--- mtl4/boost/numeric/mtl/operation/trans.hpp 2007-10-26 05:33:04.000000000 -0400
+++ mtl4-fixed/boost/numeric/mtl/operation/trans.hpp 2007-12-21 23:50:51.000000000 -0500
@@ -49,7 +49,7 @@
inline trans(const Value& v)
{
return sfunctor::trans<Value, typename traits::algebraic_category<Value>::type>::apply(v);
-};
+}

} // namespace mtl

_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/

Re: MTL: MTL development road

2007-12-03T02:27:47Z

On Mon, 2007-12-03 at 10:53 +0100, Peter Gottschling wrote:
> Hi Michael,
>
> Michael Smolsky schrieb:
> > Hi,
> >
> > MTL looks like a very promising development. I am considering
> > switching my project to MTL.

Snip

> That's an excellent idea!
> > * Is there a reason, why standard matrix-algebraic operations, such as
> > matrix products, addition, etc are not expressed in terms of C++
> > operators (*, +, etc)?
> >
> Because it is not trivial to dispatch between all the possible
> multiplication operations. If one says alpha * A it should return a
> scaled view on a matrix, A * x it should perform a matrix vector
> product, ...
>
> Furthermore, if one has three matrices that approach memory size, one
> don't want to set up another as temporary. Another reason is
> genericity, mult(A, B, C) allows you to store the temporary results
> directly into C whether it is the same type as A and B or not. If you
> have an operator A * B you don't know what type C has and if A and B
> have different types what type should have the temporary matrix (that
> might or might not shallow copied to C (but this is already another
> non-trivial problem)).
>
> However, MTL4 supports C= A * B;
> - without creating temporaries;
> - while supporting arbitrary type combinations; and
> - performing dgemm, sgemm, ... if the type combination allows.

We use MTL in real-time applications. There it is extremely important to
avoid temporaries, since allocation / freeing of memory will quickly
lead to unpredictable timing. A "fancier" matrix library like newmat is
not compatible with real-time software.

--
M.M. (René) van Paassen | ______o____/_| M.M.vanPaassen@...
Associate Professor <[___\_\_-----< t/f: +31 15 27 85370/86480
Aerospace Engineering | o' http://www.cs.lr.tudelft.nl
Delft University of Technology Kluijverweg 1, NL-2629 HS Delft
###########################################

This message has been scanned by F-Secure Anti-Virus for Microsoft Exchange.
For more information, connect to http://www.f-secure.com/

_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/

Re: MTL: MTL development road

2007-12-03T01:53:38Z

Hi Michael,

Michael Smolsky schrieb:
> Hi,
>
> MTL looks like a very promising development. I am considering
> switching my project to MTL.
>
Thank you for the compliment and we are glad about your intentions.
> I would like to know MTL maintainers' opinion about the immediate
> future of MTL. The first set of questions, that comes to mind is:
>
> * Is there a document, that describes the roadmap for MTL?
>
Currently not. We have to clarify the funding situation for the project
at the moment. Otherwise the roadmap would be pure fiction.

> * I noticed, that MTL currently doesn't do bound checking: my test
> code, in which I deliberately assigned a value to an out-of-bound
> element of a dense matrix didn't crash in any way. I find this pretty
> inconvenient. Having bound checking (when compiled under -DDEBUG) can
> save a lot of development time. From talking to Java folks, I noticed,
> that presence such checks in Java libraries, and lack of those in STL
> implementations was one of the reasons, that drove them away from C++
> and towards Java. Is there a plan to add bound checking? A place, such
> as dense2D.h:rect_offset<..>::elt(..) would be my guess?
>

Are you talking about MTL2 or MTL4? For the latter this is mostly
prepared and was later disabled for simple annoying reason. For STL
like implementations one needs end-iterators for half-open intervals.
By nature they are past the end and unfortunately I implemented them by
&A[i][end_column]. Therefore, the library provokes the exception
itself. This annoying detail will be fixed soon.

Thank you for pointing out the importance of it.
> * It's great, that MTL supports some LAPACK/BLAS. With all due respect
> to optimizations, that are built into MTL, some vendors optimize these
> libraries to their hardware, and that perhaps can beat sophisticated
> template unwrapping techniques. Hence my questions:
>
> ** Is there a plan to support sparse BLAS routines? Maybe MTL supports
> those now, and I just missed it?
>
Is there a sparse BLAS implementation that is worth the effort? AFAIK
there is only a prototype in C++ and no hand-tuned implementation. Even
if there was I doubt that sparse operations could be implemented much
faster than in C++. For dense matrix multiplication the run-time is
practically independent from memory bandwidth if you use caches well and
the performance depends on a filigree register choreography that
compilers these days do quite well but still not with the same
perfection as hand-tuned code. For sparse matrix operations the run time
is mainly determined by loading the matrix/matrices and there is not
much one can do with hardware optimization, eventually some prefetching
but even this is questionable since there would be not spare bandwidth
for the it.

I must admit that I only skimmed the sparse BLAS proposal but there are
several operations supported by MTL that are not part of sparse BLAS,
e.g., sparse matrix multiplication (which one could use in algebraic
multigrid methods), sparse dense matrix product (which could be used to
multiply a sparse matrix simultaneously with multiple vectors). I
neither recall that the standard addressed the fundamental question how
the matrix should be filled. MTL4 has a very efficient implementation
for this and the matrix insertion can be done in generic fashion, i.e.
sparse and dense matrices can be filled in the same way.
> ** It seems, that matrix/vector products are handled by MTL itself,
> rather than by BLAS. Is this the case? Is there a doc, that describes
> what functionalities of BLAS/LAPACK are utilized, rather than being
> implemented within MTL?
>
In MTL2 are 2 files you can just look it up and in MTL4 there is
currently only matrix multiplication, but this can be done in
user-transparent manner (a feature I haven't found in another generic
library yet).

> AMD promises to start selling FireStream boards fairly soon
> (http://ati.amd.com/products/streamprocessor/specs.html), and those
> might turn out to be ideal platform for running applications, built on
> top of MTL. Has anybody tried linking an MTL-based application to
> ACML, which is AMD's implementation of BLAS/LAPACK?
>
> Regarding the previous questions, it looks like MTL does not use
> utilize exceptions for reporting BLAS/LAPACK error codes. Is there a
> plan to wrap LAPACK's INFO variable into exceptions?
>

That's an excellent idea!
> * Is there a reason, why standard matrix-algebraic operations, such as
> matrix products, addition, etc are not expressed in terms of C++
> operators (*, +, etc)?
>
Because it is not trivial to dispatch between all the possible
multiplication operations. If one says alpha * A it should return a
scaled view on a matrix, A * x it should perform a matrix vector
product, ...

Furthermore, if one has three matrices that approach memory size, one
don't want to set up another as temporary. Another reason is
genericity, mult(A, B, C) allows you to store the temporary results
directly into C whether it is the same type as A and B or not. If you
have an operator A * B you don't know what type C has and if A and B
have different types what type should have the temporary matrix (that
might or might not shallow copied to C (but this is already another
non-trivial problem)).

However, MTL4 supports C= A * B;
- without creating temporaries;
- while supporting arbitrary type combinations; and
- performing dgemm, sgemm, ... if the type combination allows.

Nevertheless there is still a lot to do. This said, if somebody
volunteers to extend the interface to BLAS and LAPACK it would be more
then welcome. This could also be done as a master thesis, the principle
is clear and I would be glad to introduce it to the student and if
desired talk to his/her professor.

Best Regards,

Peter
_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/

MTL: MTL development road

2007-12-01T12:19:44Z

Hi,

MTL looks like a very promising development. I am considering switching my project to MTL.

I would like to know MTL maintainers' opinion about the immediate future of MTL. The first set of questions, that comes to mind is:

* Is there a document, that describes the roadmap for MTL?

* I noticed, that MTL currently doesn't do bound checking: my test code, in which I deliberately assigned a value to an out-of-bound element of a dense matrix didn't crash in any way. I find this pretty inconvenient. Having bound checking (when compiled under -DDEBUG) can save a lot of development time. From talking to Java folks, I noticed, that presence such checks in Java libraries, and lack of those in STL implementations was one of the reasons, that drove them away from C++ and towards Java. Is there a plan to add bound checking? A place, such as dense2D.h:rect_offset<..>::elt(..) would be my guess?

* It's great, that MTL supports some LAPACK/BLAS. With all due respect to optimizations, that are built into MTL, some vendors optimize these libraries to their hardware, and that perhaps can beat sophisticated template unwrapping techniques. Hence my questions:

** Is there a plan to support sparse BLAS routines? Maybe MTL supports those now, and I just missed it?

** It seems, that matrix/vector products are handled by MTL itself, rather than by BLAS. Is this the case? Is there a doc, that describes what functionalities of BLAS/LAPACK are utilized, rather than being implemented within MTL?

AMD promises to start selling FireStream boards fairly soon (http://ati.amd.com/products/streamprocessor/specs.html), and those might turn out to be ideal platform for running applications, built on top of MTL. Has anybody tried linking an MTL-based application to ACML, which is AMD's implementation of BLAS/LAPACK?

Regarding the previous questions, it looks like MTL does not use utilize exceptions for reporting BLAS/LAPACK error codes. Is there a plan to wrap LAPACK's INFO variable into exceptions?

* Is there a reason, why standard matrix-algebraic operations, such as matrix products, addition, etc are not expressed in terms of C++ operators (*, +, etc)?

Thank you for programming and supporting something, that looks like a very useful library.

Regards,

Michael.

--

Over 2 Million Holiday Gift Ideas - Take a Look!
mail.com Shopping

_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/

Re: MTL: Symmetric Matrix Resize

2007-11-04T19:44:06Z

PS = Pete Stieber
PS>> MTL version: mtl-2.1.2-23.tenative
PS>> Compiler: g++ (GCC) 4.1.2 20070925 (Red Hat 4.1.2-27)
PS>> Platform: Fedora 7 x86_64
PS>>
PS>> I'm trying to resize a symmetric matrix using a simple
PS>> test jig called SymmetricPackedResize.cpp ...
PS>>
PS>> #include <mtl/mtl.h>
PS>>
PS>> int main()
PS>> {
PS>> typedef mtl::matrix<
PS>> double,
PS>> mtl::symmetric<mtl::lower>,
PS>> mtl::packed<> >::type TDSymmetricMatrix;
PS>>
PS>> TDSymmetricMatrix Symmetric(6, 6);
PS>>
PS>> Symmetric.resize(16, 16);
PS>>
PS>> return 0;
PS>> }
PS>>
PS>> I'm compiling with
PS>>
PS>> g++ -IPathToMtlHeaders SymmetricPackedResize.cpp
PS>>
PS>> The first set of errors I receive follow...
Snip
PS>> So the packed_offset class has the dim data member
PS>> declared private. Looking at the rect_offset class
PS>> for which resize actually works, there is a
PS>> suspicious commenting out for the private: scope
PS>> qualifier ahead of the data members.
PS>>
PS>> I could do the same in packed_offset, but this
PS>> would still leave me with the "packed_offset<...>
PS>> has no member named ld. It looks like resize
PS>> will only work for rect_offset storage.
PS>>
PS>> Does anyone see a workaround for this?
PS>>
PS>> Peter, does version 4 deal with this?
PS>>
PS>> How do I become an alpha tester of version 4?

PS> Instead of explicitly setting the values of the
PS> offset class in the resize function, I changed
PS> the code to call a new resize member function of
PS> the offset class. The following code implements
PS> resize for the rect_offset class by performing
PS> the same functionality that was in the dense2D
PS> class resize function.
PS>
PS> inline void resize(size_type m, size_type n)
PS> {
PS> dim = dim_type(m, n);
PS> ld = n;
PS> }
PS>
PS> Now I need to add a resize function to the
PS> packed_offset class. Here's my first cut...
PS>
PS> inline void resize(size_type m, size_type n)
PS> {
PS> dim = dim_type(m, n);
PS> }
PS>
PS> but the packed_offset class has an addition
PS> data member called bw. Can anyone explain the
PS> purpose of this data member so I can figure out
PS> if I need to alter it in packed_offset::resize?

OK. I believe bw stands for band width, as in which bands are required
with respect to the storage scheme. My new packed_offset resize now
looks like the following:

inline void resize(size_type m, size_type n)
{
dim = dim_type(m, n);
bw = band_type(0, n);
}

since my symmetric matrix is full. I also had to modify the storage
allocation call in dense2d::resize from

rep_ptr newdata = new reptype(Offset::size(m, n, 0, 0));

to

rep_ptr newdata = new reptype(Offset::size(m, n, 0, n));

to generate the proper amount of storage for a full symmetric matrix. A
general dense matrix using rect_offset doesn't use the last two
parameters, so it still works.

Unfortunately, I couldn't figure out how to copy the previously existing
existing matrix elements to the new storage. The (size_type, size_type)
operator doesn't work properly at this level of scope. If anyone can
help me with that it would be greatly appreciated.

I also noticed that the following code to zero out any extra matrix
elements if the new matrix is bigger

for (; i < m; ++i)
for (; j < n; ++j)
(*newdata)[i * n + j] = T();

was incorrect and seemed to be unneeded because the allocation code
seems to zero out the data.

If anyone is interested in a patch to the dense2d.h file that will
resize both general and symmetric matrices, but will not copy old
elements, I believe I have one.

If anyone can help me figure out how to copy the old elements, I would
greatly appreciate it.

Pete

_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/

Re: MTL: Test fails in MTL4 alpha 1

2007-10-30T10:37:35Z

Hi Jari,

You are right, my hand-written CCS matrix was wrong. Which shows how
useful the inserter is. ;-) The version on the web page is fixed.

I also threw out the test for Fortran-indexed matrices since I stopped
supporting them some day.

Thank you for helping us to make the lib bug-free.

Cheers,

Peter

Jari Häkkinen schrieb:

> Hi,
>
> One of the tests fails on my mac. The problem seems to be related to the
> function
>
> template <typename Matrix>
> void raw_copy_test(Matrix& matrix, col_major)
>
> in compressed2D_inserter.cpp. The index array created in the problematic
> function looks strange. I get the output attached below.
>
> I am using:
>
> boost 1.34.1
> i686-apple-darwin8-g++-4.0.1 (GCC) 4.0.1 (Apple Computer, Inc. build 5363)
>
>
> The other tests works fine.
>
>
> Cheers,
>
> Jari
>
>
>
> # scons -D . check=1 with-blas=1
>
> scons: Entering directory `/projects/osd/mtl4'
> scons: Reading SConscript files ...
> -DMTL_ASSERT_FOR_THROW
> scons: done reading SConscript files.
> scons: Building targets ...
> g++ -o libs/numeric/mtl/test/compressed2D_inserter.o -c -Wno-long-double
> -g -DMTL_ASSERT_FOR_THROW -DNDEBUG -Wno-long-double -I.
> libs/numeric/mtl/test/compressed2D_inserter.cpp
> g++ -o libs/numeric/mtl/test/compressed2D_inserter
> libs/numeric/mtl/test/compressed2D_inserter.o
> /projects/osd/mtl4/libs/numeric/mtl/test/compressed2D_inserter
> [ 0 1 0 0 2 0]
> [ 0 0 3 0 0 0]
> [ 4 5 6 7 0 0]
> [ 0 0 0 8 0 9]
> [ 0 10 0 11 12 0]
> [ 0 13 0 14 15 16]
> [ 0 0 0 0 0 0]
> [ 0 0 17 0 0 18]
>
>
> [ 0 1 0 31 2 0]
> [ 0 0 3 0 0 0]
> [ 4 5 6 7 0 0]
> [ 0 0 0 33 0 9]
> [ 0 10 0 11 37 0]
> [ 0 13 0 14 15 16]
> [ 35 0 0 0 0 0]
> [ 0 0 17 0 0 18]
>
>
> [ 0 1 27 31 2 0]
> [ 0 0 3 0 0 0]
> [ 4 5 27 7 22 26]
> [ 0 28 0 33 0 9]
> [ 0 10 54 11 37 0]
> [ 0 13 0 14 15 16]
> [ 35 23 0 0 0 0]
> [ 0 0 43 0 0 18]
>
>
>
> matrix[2][2] = 28
> [ 0 0 0 0 0 0]
> [ 0 0 1 0 0 2]
> [ 0 0 0 3 0 0]
> [ 0 4 5 6 7 0]
> [ 0 0 0 0 8 0]
> [ 0 0 10 0 11 12]
> [ 0 0 13 0 14 15]
> [ 0 0 0 0 0 0]
>
>
> [ 0 0 0 0 0 0]
> [ 0 0 1 0 31 2]
> [ 0 0 0 3 0 0]
> [ 0 4 5 6 7 0]
> [ 0 0 0 0 33 0]
> [ 0 0 10 0 11 37]
> [ 0 0 13 0 14 15]
> [ 0 35 0 0 0 0]
>
>
> [ 0 0 0 0 0 0]
> [ 0 0 1 27 31 2]
> [ 0 0 0 3 0 0]
> [ 0 4 5 27 7 22]
> [ 0 0 28 0 33 0]
> [ 0 0 10 54 11 37]
> [ 0 0 13 0 14 15]
> [ 0 35 23 0 0 0]
>
>
>
> matrix[2][2] = 28
> [ 0 1 0 0 2 0]
> [ 0 0 3 0 0 0]
> [ 4 5 6 7 0 0]
> [ 0 0 0 8 0 0]
> [ 0 10 0 0 12 0]
> [ 0 13 0 0 15 16]
> [ 0 0 0 0 0 0]
> [ 0 0 17 0 0 18]
>
>
> [ 0 1 0 31 2 0]
> [ 0 0 3 0 0 0]
> [ 4 5 6 7 0 0]
> [ 0 0 0 33 0 0]
> [ 0 10 0 0 37 0]
> [ 0 13 0 0 15 16]
> [ 35 0 0 0 0 0]
> [ 0 0 17 0 0 18]
>
>
> [ 0 1 27 31 2 0]
> [ 0 0 3 0 0 0]
> [ 4 5 27 7 22 26]
> [ 0 28 0 33 0 0]
> [ 0 10 54 0 37 0]
> [ 0 13 0 0 15 16]
> [ 35 23 0 0 0 0]
> [ 0 0 43 0 0 18]
>
>
>
> matrix[2][2] = 28
> /usr/local/include/boost/test/minimal.hpp(122): exception "memory access
> violation" caught in function: 'int main(int, char**)'
>
> **** Testing aborted.
> **** 1 error detected
> [scons: *** [libs/numeric/mtl/test/compressed2D_inserter] Error 201
> scons: building terminated because of errors.
>
> _______________________________________________
> This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/
>
>

--
Dr. rer. nat. Peter Gottschling

Dresdner Grundwasserforschungszentrum e.V. (DGFZ)
Meraner Str. 10
01217 Dresden
Germany
.....................................................................
phone: +49 (0) 351 / 4050651
fax: +49 (0) 351 / 4050679
web: http://www.dgfz.de
.....................................................................
VR 210 Ust.-Id DE140134319
GF: Dr. F. Börner
.....................................................................

Veranstaltungshinweise:

Weiterbildung im DGFZ e.V.

14.-16.11.2007 Modellierung mit PCGEOFIM

Programm und Termine:
http://www.dgfz.de/aktuell

06./07.11.2007

XVIII. Sächsisches Altlastenkolloquium
und regionaler Workshop des BMBF-Förderschwerpunktes KORA

"Natürliche Schadstoffminderungsprozesse bei der Altlastenbearbeitung"

Programm und Anmeldung:
http://www.dgfz.de/index_altlastenkolloquium

_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/

Re: MTL: Symmetric Matrix Resize

2007-10-30T10:23:16Z

Peter J. Stieber wrote:

> MTL version: mtl-2.1.2-23.tenative
> Compiler: g++ (GCC) 4.1.2 20070925 (Red Hat 4.1.2-27)
> Platform: Fedora 7 x86_64
>
> I'm trying to resize a symmetric matrix using a simple test jig called
> SymmetricPackedResize.cpp ...
>
> #include <mtl/mtl.h>
>
> int main()
> {
> typedef mtl::matrix<
> double,
> mtl::symmetric<mtl::lower>,
> mtl::packed<> >::type TDSymmetricMatrix;
>
> TDSymmetricMatrix Symmetric(6, 6);
>
> Symmetric.resize(16, 16);
>
> return 0;
> }
>
> I'm compiling with
>
> g++ -IPathToMtlHeaders SymmetricPackedResize.cpp
>
> The first set of errors I receive follow. I've tried to separate them
> onto different lines for readability:
>
> mtl-2.1.2-23/mtl/dense2D.h: In member function ‘void mtl::dense2D<T,
> OffsetGen, MM, NN>::resize(typename mtl::generic_dense2D<std::vector<T,
> std::allocator<_CharT> >, mtl::refcnt_ptr<std::vector<T,
> std::allocator<_CharT> > >, OffsetGen, MM, NN>::size_type, typename
> mtl::generic_dense2D<std::vector<T, std::allocator<_CharT> >,
> mtl::refcnt_ptr<std::vector<T, std::allocator<_CharT> > >, OffsetGen,
> MM, NN>::size_type) [with T = double, OffsetGen =
> mtl::gen_packed_offset<0, 0>, int MM = 0, int NN = 0]’:
> mtl-2.1.2-23/mtl/matrix_implementation.h:616: instantiated from ‘void
> mtl::row_matrix<TwoDGen, IndexerGen>::resize(typename
> mtl::matrix_implementation<TwoDGen, IndexerGen>::size_type, typename
> mtl::matrix_implementation<TwoDGen, IndexerGen>::size_type) [with
> TwoDGen = mtl::gen_dense2D<double, mtl::gen_packed_offset<0, 0>, 0, 0>,
> IndexerGen = mtl::gen_banded_indexer<mtl::row_orien, 0, 0, long unsigned
> int>]’
>
> SymmetricPackedResize.cpp:12: instantiated from here
>
> mtl-2.1.2-23/mtl/dense2D.h:548: error: ‘mtl::dimension<long unsigned
> int, 0, 0> mtl::packed_offset<long unsigned int, 0, 0>::dim’ is private
>
> mtl-2.1.2-23/mtl/dense2D.h:1235: error: within this context
>
> mtl-2.1.2-23/mtl/dense2D.h:1236: error: ‘class mtl::packed_offset<long
> unsigned int, 0, 0>’ has no member named ‘ld’
>
> So the packed_offset class has the dim data member declared private.
> Looking at the rect_offset class for which resize actually works, there
> is a suspicious commenting out for the private: scope qualifier ahead of
> the data members.
>
> I could do the same in packed_offset, but this would still leave me with
> the "packed_offset<...> has no member named ld. It looks like resize
> will only work for rect_offset storage.
>
> Does anyone see a workaround for this?
>
> Peter, does version 4 deal with this?
>
> How do I become an alpha tester of version 4?

Instead of explicitly setting the values of the offset class in the
resize function, I changed the code to call a new resize member function
of the offset class. The following code implements resize for the
rect_offset class by performing the same functionality that was in the
dense2D class resize function.

inline void resize(size_type m, size_type n)
{
dim = dim_type(m, n);
ld = n;
}

Now I need to add a resize function to the packed_offset class. Here's
my first cut...

inline void resize(size_type m, size_type n)
{
dim = dim_type(m, n);
}

but the packed_offset class has an addition data member called bw. Can
anyone explain the purpose of this data member so I can figure out if I
need to alter it in packed_offset::resize?

Pete

_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/

MTL: MTL 4 released

2007-10-30T10:07:15Z

Hi everybody on the two mailing lists,

We are happy announce that we now release the Alpha version of MTL 4. I
fixed the few bugs that I heard of in the last weeks. The latest
version can be downloaded at
http://www.osl.iu.edu/research/mtl/mtl4/download.php3

I like to thank for all the helpful and encouraging comments I received.

To separate concerns a little I want to ask sending mails concerning MTL
4 to mtl4@.... You can subscribe at:
http://www.osl.iu.edu/mailman/listinfo.cgi/mtl4
The mailing list was fairly silent in the last years but it is now a
good moment to revive it.

Best,
Peter

_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/

MTL: Symmetric Matrix Resize

2007-10-29T15:26:08Z

MTL version: mtl-2.1.2-23.tenative
Compiler: g++ (GCC) 4.1.2 20070925 (Red Hat 4.1.2-27)
Platform: Fedora 7 x86_64

I'm trying to resize a symmetric matrix using a simple test jig called
SymmetricPackedResize.cpp ...

#include <mtl/mtl.h>

int main()
{
typedef mtl::matrix<
double,
mtl::symmetric<mtl::lower>,
mtl::packed<> >::type TDSymmetricMatrix;

TDSymmetricMatrix Symmetric(6, 6);

Symmetric.resize(16, 16);

return 0;
}

I'm compiling with

g++ -IPathToMtlHeaders SymmetricPackedResize.cpp

The first set of errors I receive follow. I've tried to separate them
onto different lines for readability:

mtl-2.1.2-23/mtl/dense2D.h: In member function ‘void mtl::dense2D<T,
OffsetGen, MM, NN>::resize(typename mtl::generic_dense2D<std::vector<T,
std::allocator<_CharT> >, mtl::refcnt_ptr<std::vector<T,
std::allocator<_CharT> > >, OffsetGen, MM, NN>::size_type, typename
mtl::generic_dense2D<std::vector<T, std::allocator<_CharT> >,
mtl::refcnt_ptr<std::vector<T, std::allocator<_CharT> > >, OffsetGen,
MM, NN>::size_type) [with T = double, OffsetGen =
mtl::gen_packed_offset<0, 0>, int MM = 0, int NN = 0]’:
mtl-2.1.2-23/mtl/matrix_implementation.h:616: instantiated from ‘void
mtl::row_matrix<TwoDGen, IndexerGen>::resize(typename
mtl::matrix_implementation<TwoDGen, IndexerGen>::size_type, typename
mtl::matrix_implementation<TwoDGen, IndexerGen>::size_type) [with
TwoDGen = mtl::gen_dense2D<double, mtl::gen_packed_offset<0, 0>, 0, 0>,
IndexerGen = mtl::gen_banded_indexer<mtl::row_orien, 0, 0, long unsigned
int>]’

SymmetricPackedResize.cpp:12: instantiated from here

mtl-2.1.2-23/mtl/dense2D.h:548: error: ‘mtl::dimension<long unsigned
int, 0, 0> mtl::packed_offset<long unsigned int, 0, 0>::dim’ is private

mtl-2.1.2-23/mtl/dense2D.h:1235: error: within this context

mtl-2.1.2-23/mtl/dense2D.h:1236: error: ‘class mtl::packed_offset<long
unsigned int, 0, 0>’ has no member named ‘ld’

So the packed_offset class has the dim data member declared private.
Looking at the rect_offset class for which resize actually works, there
is a suspicious commenting out for the private: scope qualifier ahead of
the data members.

I could do the same in packed_offset, but this would still leave me with
the "packed_offset<...> has no member named ld. It looks like resize
will only work for rect_offset storage.

Does anyone see a workaround for this?

Peter, does version 4 deal with this?

How do I become an alpha tester of version 4?

Pete

_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/

MTL: Test fails in MTL4 alpha 1

2007-10-26T14:21:29Z

Hi,

One of the tests fails on my mac. The problem seems to be related to the
function

template <typename Matrix>
void raw_copy_test(Matrix& matrix, col_major)

in compressed2D_inserter.cpp. The index array created in the problematic
function looks strange. I get the output attached below.

I am using:

boost 1.34.1
i686-apple-darwin8-g++-4.0.1 (GCC) 4.0.1 (Apple Computer, Inc. build 5363)

The other tests works fine.

Cheers,

Jari

# scons -D . check=1 with-blas=1

scons: Entering directory `/projects/osd/mtl4'
scons: Reading SConscript files ...
-DMTL_ASSERT_FOR_THROW
scons: done reading SConscript files.
scons: Building targets ...
g++ -o libs/numeric/mtl/test/compressed2D_inserter.o -c -Wno-long-double
-g -DMTL_ASSERT_FOR_THROW -DNDEBUG -Wno-long-double -I.
libs/numeric/mtl/test/compressed2D_inserter.cpp
g++ -o libs/numeric/mtl/test/compressed2D_inserter
libs/numeric/mtl/test/compressed2D_inserter.o
/projects/osd/mtl4/libs/numeric/mtl/test/compressed2D_inserter
[ 0 1 0 0 2 0]
[ 0 0 3 0 0 0]
[ 4 5 6 7 0 0]
[ 0 0 0 8 0 9]
[ 0 10 0 11 12 0]
[ 0 13 0 14 15 16]
[ 0 0 0 0 0 0]
[ 0 0 17 0 0 18]

[ 0 1 0 31 2 0]
[ 0 0 3 0 0 0]
[ 4 5 6 7 0 0]
[ 0 0 0 33 0 9]
[ 0 10 0 11 37 0]
[ 0 13 0 14 15 16]
[ 35 0 0 0 0 0]
[ 0 0 17 0 0 18]

[ 0 1 27 31 2 0]
[ 0 0 3 0 0 0]
[ 4 5 27 7 22 26]
[ 0 28 0 33 0 9]
[ 0 10 54 11 37 0]
[ 0 13 0 14 15 16]
[ 35 23 0 0 0 0]
[ 0 0 43 0 0 18]

matrix[2][2] = 28
[ 0 0 0 0 0 0]
[ 0 0 1 0 0 2]
[ 0 0 0 3 0 0]
[ 0 4 5 6 7 0]
[ 0 0 0 0 8 0]
[ 0 0 10 0 11 12]
[ 0 0 13 0 14 15]
[ 0 0 0 0 0 0]

[ 0 0 0 0 0 0]
[ 0 0 1 0 31 2]
[ 0 0 0 3 0 0]
[ 0 4 5 6 7 0]
[ 0 0 0 0 33 0]
[ 0 0 10 0 11 37]
[ 0 0 13 0 14 15]
[ 0 35 0 0 0 0]

[ 0 0 0 0 0 0]
[ 0 0 1 27 31 2]
[ 0 0 0 3 0 0]
[ 0 4 5 27 7 22]
[ 0 0 28 0 33 0]
[ 0 0 10 54 11 37]
[ 0 0 13 0 14 15]
[ 0 35 23 0 0 0]

matrix[2][2] = 28
[ 0 1 0 0 2 0]
[ 0 0 3 0 0 0]
[ 4 5 6 7 0 0]
[ 0 0 0 8 0 0]
[ 0 10 0 0 12 0]
[ 0 13 0 0 15 16]
[ 0 0 0 0 0 0]
[ 0 0 17 0 0 18]

[ 0 1 0 31 2 0]
[ 0 0 3 0 0 0]
[ 4 5 6 7 0 0]
[ 0 0 0 33 0 0]
[ 0 10 0 0 37 0]
[ 0 13 0 0 15 16]
[ 35 0 0 0 0 0]
[ 0 0 17 0 0 18]

[ 0 1 27 31 2 0]
[ 0 0 3 0 0 0]
[ 4 5 27 7 22 26]
[ 0 28 0 33 0 0]
[ 0 10 54 0 37 0]
[ 0 13 0 0 15 16]
[ 35 23 0 0 0 0]
[ 0 0 43 0 0 18]

matrix[2][2] = 28
/usr/local/include/boost/test/minimal.hpp(122): exception "memory access
violation" caught in function: 'int main(int, char**)'

**** Testing aborted.
**** 1 error detected
[scons: *** [libs/numeric/mtl/test/compressed2D_inserter] Error 201
scons: building terminated because of errors.

_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/