dolfin team mailing list archive

[Johan Hake <johan.hake@xxxxxxxxx>]

Isn't there a way to initialize a paralell vector with arbitrary distribution 
pattern?

There is a way to general a local range from a given entity size using 
MPI.local_range. One should be able to generate a distributed Vector with this 
information.

Johan

On Monday January 24 2011 05:17:21 Joachim Haga wrote:
> Public bug reported:
> 
> The python code for __mul__() assumes an empty vector suitable for
> mult() can be created from just the dimensions of the matrix:
> 
>     if type(other) == Vector:
>         ret = Vector(self.size(0))
> 
> This is not true in parallel. And it won't help to copy the input
> vector, as that only works for square matrices.
> 
> How to create a vector which is compatible with a given matrix at this
> point (only the matrix available, not the sparsity pattern)? Would it
> perhaps be acceptable to add methods create_rhs_vector(),
> create_lhs_vector() to GenericMatrix?
> 
> ** Affects: dolfin
>      Importance: Undecided
>          Status: New