[Dune] [#417] SuperLU< BCRSMatrix... > produces memory errors in constructor (?)

Dune dune at dune-project.org
Tue Jul 22 16:05:52 CEST 2008


THIS IS AN AUTOMATED MESSAGE, DO NOT REPLY.

A new Flyspray task has been opened.  Details are below. 

User who did this - Uli Sack (usack) 

Attached to Project - Dune
Summary - SuperLU< BCRSMatrix... > produces memory errors in constructor (?)
Task Type - Bug Report
Category - ISTL
Status - Unconfirmed
Assigned To - 
Operating System - Linux
Severity - Low
Priority - Normal
Reported Version - SVN
Due in Version - Undecided
Due Date - Undecided
Details - calling the constructor of Dune::SuperLU< BCRSMatrix< FieldMatrix< T, n, m >, A > > causes memory errors somewhere in my BLAS: (valgrind)

==21436== Use of uninitialised value of size 4
==21436==    at 0x40A127C: (within /usr/lib/atlas/sse2/libblas.so.3.0)
==21436==
==21436== Use of uninitialised value of size 4
==21436==    at 0x40991EC: (within /usr/lib/atlas/sse2/libblas.so.3.0)
==21436==
==21436== Use of uninitialised value of size 4
==21436==    at 0x40991F2: (within /usr/lib/atlas/sse2/libblas.so.3.0)

The (more or less) minimal code example to produce this (on my machine) is:

#ifdef HAVE_CONFIG_H
# include "config.h"     
#endif

#include <iostream>

#include <dune/istl/bcrsmatrix.hh>
#include <dune/istl/matrixindexset.hh>
#include <dune/istl/bvector.hh>
#include <dune/istl/solvers.hh>
#include <dune/istl/superlu.hh>

using namespace Dune;

int main (int argc, char *argv[]) try
{

	typedef Dune::BCRSMatrix< Dune::FieldMatrix<double,1,1> > Matrix;
	typedef Dune::BlockVector< Dune::FieldVector<double,1> > Vector;

	const int N=9;

	Matrix A;
	MatrixIndexSet indices(9,9);

	for (int i=0; i<N; ++i) {
		for (int j=0; j<N; ++j) {
			if ((j<N/4+i and i<N/4+j) or (j==N/2+i or i==N/2+j) or (j==N-2+i or i==N-2+j)) {
				indices.add(i,j);
				std::cout << "

More information can be found at the following URL:
http://www.dune-project.org/flyspray/index.php?do=details&task_id=417

You are receiving this message because you have requested it from the Flyspray bugtracking system.  If you did not expect this message or don't want to receive mails in future, you can change your notification settings at the URL shown above.




More information about the Dune mailing list