dolfin team mailing list archive

[Joachim Berdal Haga <jobh@xxxxxxxxx>]

>
> > It seems that this is only implemented for rows, not for columns.
>
> Do you mean GenericMatrix::local_range? It takes a dim argument,
>
>  std::pair<uint, uint> GenericMatrix::local_range(uint dim);


Yes, but:

   if (dim == 1)
    error("Cannot compute columns range for Epetra matrices.");

...

      if (dim == 1)
        error("Cannot compute columns range for PETSc matrices.");

-j.