dolfin team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

On Mon, Sep 13, 2010 at 06:11:13AM -0000, Evan Lezar wrote:
> Question #125019 on DOLFIN changed:
> https://answers.launchpad.net/dolfin/+question/125019
>
>     Status: Answered => Open
>
> Evan Lezar is still having a problem:
> Hi
>
> I think Anders is correct. I have just had a look at the SystemAssembler
> code, and it does seem to also place ones on the diagonals of the matrix
> to which the boundary conditions are being applied. In addition, I am
> not solving a linear system, but eigenvalue problems and as such b is a
> Matrix not a Vector.
>
> At present I am manually removing the rows and columns for the matrices
> that correspond to (zero) Dirichlet conditions and then reconstructing
> the correct solution vector from the eigenvector that I obtain. This is
> fine for testing and getting the rest of the system running.
>
> Thanks for the help so far.

How do you remove the rows and columns?

If there is good support for doing this in the backends (PETSc,
Trilinos), we might consider adding a wrapper for it.

--
Anders