Hello Marc,
as far as I understood, the FilteredMatrix class applies the Dirichlet
boundary conditions on the fly when it performs the vmult and Tvmult
functions. However in the function factorize of SparseDirectUMFPACK you
copy the matrix into the memory format of UMFPACK first and then compute
the LU-decomposition of the matrix.
Therefore I believe one way to go is the following:
1. Write a new function factorize for the FilteredMatrix. Here you have
to apply the Dirichlet boundary values when copying over to UMFPACK
memory. This should work similar to the
MatrixTools::apply_boundary_values function, except that you store the
entries in the memory for UMFPACK instead of the matrix.
2. Write a new function solve, which takes into account the modification
of the right-hand side.
Sounds good?
Best Regards,
Markus
Am 24.11.2011 16:27, schrieb Marc Secanell Gallart:
Dear deal.ii users,
I am trying to solve a problem using the SparseDirectUMFPACK class.
Currently, I am applying the Dirichlet boundary conditions to my
matrix and right hand side using the class FilteredMatrix in order to
be able to re-use my system matrix. Unfortunately, I found out that
SparseDirectUMFPACK has not been instantiated for MATRIX templates of
type FilteredMatrix. I was wondering, could anyone please give me some
information on what I would need to do to extend the
SparseDirectUMFPACK interface to be able to use objects of class
FilteredMatrix with SparseDirectUMFPACK?
Thank you,
Marc
--
--------------------------------------------------------------------------------------------------------------------
Marc Secanell, Ph.D., P.Eng.
Assistant Professor
Energy Systems Design Laboratory
Department of Mechanical Engineering
University of Alberta
4-31F Mechanical Engineering Building
Edmonton AB T6G 2G8 Canada
Phone: (780) 492 6961
E-mail: [email protected] <mailto:[email protected]>
Confidentiality notice: Emails sent to [email protected]
<mailto:[email protected]> are forwarded to a Google-supported
email account. Emails, attachments, documents and calendaring data,
stored, received and sent through this University account may be
transmitted to and stored on servers outside of the University,
Alberta or Canada. The University of Alberta cannot and does not
guarantee protection against the possible disclosure of your data
including, without limitation, against possible secret disclosures of
data to a foreign authority as a consequence of foreign laws.
--------------------------------------------------------------------------------------------------------------------
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
