Re: [glas] unallocated fill
From: Toon Knapen (toon.knapen_at_[hidden])
Date: 2005-10-03 15:56:35
Neal Becker wrote:
> It is important to separate allocation and initialization for HPC. Here is
> what I believe. The default construction should not perform initialization.
> It looks to me that glas already follows this.
However I'm not sure the current implementation is correct. Currently
the constructor of the dense_vector just calls the 'allocate(n)' member
of the allocator. What if an element in the vector will be assigned to
before fill is being called on the vector? For example following use-case:
glas::dense_vector<double> v(3) ;
v = 5 ;
So in this scenario an array that can contain 3 double's is allocated.
However the allocate member function only allocates but does not
construct any objects in the array (see table 32 in paragraph 20.1.5).
The second statement in the use-case above however relies on an
assignment-operator. But the assignment-operator is called on a
memory-address where no object is constructed yet.
For POD's I guess this will usually work fine but... IIUC uBLAS solves
this by calling an array-placement-new on the allocate'd array. Is this
indeed the best way to go ?