Re: [ublas] indirect_array<> and resize()

Subject: Re: [ublas] indirect_array<> and resize()
From: Ralf Denzer (Ralf.Denzer_at_[hidden])
Date: 2008-09-21 06:11:09

Gunter, thanks for your fast reply!

> resizing an indirect_array is not supported. You have to use the 
> compose() functions to create a new indirect_array with more indices. 
> The reason is that the behavior of an indirect_vector or 
> indirect_matrix would be undefined when the size of the underlying 
> indirect_array is changes because the indirect_array loses some indices 
> to gets new indices which might be out of range.

Ok, but I can do 

ublas::matrix<double> M(5,5);
ublas::indirect_array<> ia(3);
ia(0) = 1; ia[1] = 2; ia[2] = 4;
ublas::matrix_indirect<ublas::matrix<double> > Mi(M, ia,ia);
M.resize(4,4);
std::cout << "after M.resize() Mi=\n" << Mi << std::endl; // bad index "bombs"

In any case, I have to take care about ranges during runtime.

The real cause for my question is that I like to initialize something like 

ublas::vector<ublas::indirect_array<> > vec_ia;
vec_ia.resize(10);
for (unsigned i = 0; i < 10; ++i)
{
    vec_ia(i).resize( n(i) );  // n(i) are user given unsigned int's (size_t's)
    vec_ia(i)(1:n(i)) = ...     // ok, matlab notation
}

With this data structure I'd like to do the usual gathering/scattering 
process in a finite element code. For element "i" the corresponding indices are 
stored in vec_ia(i) . Sorry, that I didn't mentioned this earlier.

Bye 

Ralf 

____________________________________________________________________
Psssst! Schon vom neuen WEB.DE MultiMessenger gehört? 
Der kann`s mit allen: http://www.produkte.web.de/messenger/?did=3123