Re: [deal.II] Massively parallel & inverse matrix solver

2021-09-01 Thread Jean-Paul Pelteret 
Dear Hermes,

The equivalent classes to SparseDirectUMFPACK for parallel linear algebra would 
be:
- TrilinosWrappers::SolverDirect 

 : You can choose the implementation through the AdditionalData 

 struct that is passed into the class constructor.
- PETScWrappers::SparseDirectMUMPS 


Best,
Jean-Paul

> On 1. Sep 2021, at 20:10, Daniel Arndt  wrote:
> 
> Hermes,
> 
> The SparseDirectUMFPACK solver doesn't work with MPI parallel matrices. You 
> will need to find another solver based on the matrix class you are using. 
> Standard choices would be GMRES or CG but I am assuming that your linear 
> system might also be non-symmetric or indefinite.
> 
> Best,
> Daniel
> 
> 
> 
> Am Mi., 1. Sept. 2021 um 11:53 Uhr schrieb Hermes Sampedro 
> mailto:hermesampe...@gmail.com>>:
> Dear all,
> 
> I am trying to solve a similar problem as step-29 but including the massive 
> parallel solution as presented in step-40.
> I would like to ask how the solve() function would be in this case.
> 
> In step-29:
> //LU decomposition and inverse matrix
> SparseDirectUMFPACK A_direct;
> A_direct.initialize(system_matrix);
> 
> //Solve with inverse matrix
> A_direct.vmult(solution, system_rhs);
> When applying massive parallel computations, my initial guess was:
> 
> 
> LA::MPI::Vector
> completely_distributed_solution(locally_owned_dofs,mpi_communicator);
> 
> //LU decomposition and inverse matrix
> SparseDirectUMFPACK A_direct;
> A_direct.initialize(system_matrix);
> 
> //Solve with inverse matrix
> A_direct.vmult(completely_distributed_solution, system_rhs);
>  locally_relevant_solution = completely_distributed_solution;
> 
> 
>  However, I get a no matching member function for call to 'vmult' error. 
> 
> How can I use this LU decomposition and inverse matrix solver using MPI?
> 
> 
> 
> Thank you,
> 
> Regards,
> 
> H
> 
> 
> 
> -- 
> The deal.II project is located at http://www.dealii.org/ 
> 
> For mailing list/forum options, see 
> https://groups.google.com/d/forum/dealii?hl=en 
> 
> --- 
> You received this message because you are subscribed to the Google Groups 
> "deal.II User Group" group.
> To unsubscribe from this group and stop receiving emails from it, send an 
> email to dealii+unsubscr...@googlegroups.com 
> .
> To view this discussion on the web visit 
> https://groups.google.com/d/msgid/dealii/ded2e156-7bd7-4474-aa5f-fb57e8540926n%40googlegroups.com
>  
> .
> 
> -- 
> The deal.II project is located at http://www.dealii.org/ 
> 
> For mailing list/forum options, see 
> https://groups.google.com/d/forum/dealii?hl=en 
> 
> --- 
> You received this message because you are subscribed to the Google Groups 
> "deal.II User Group" group.
> To unsubscribe from this group and stop receiving emails from it, send an 
> email to dealii+unsubscr...@googlegroups.com 
> .
> To view this discussion on the web visit 
> https://groups.google.com/d/msgid/dealii/CAOYDWbKG87reQPBrhZu0RwoZHrMXMtVWS9Ui6EoFCWtmemG1vQ%40mail.gmail.com
>  
> .

-- 
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see 
https://groups.google.com/d/forum/dealii?hl=en
--- 
You received this message because you are subscribed to the Google Groups 
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/6B045CAD-E179-4E5F-96B6-2CC00891CA2F%40gmail.com.

Re: [deal.II] Massively parallel & inverse matrix solver

2021-09-01 Thread Hermes Sampedro 
Thank you for the answer. I am solving for real and imaginary parts.
What solver would you suggest using? Is it any example in any of the steps?

Thank you
Hermes.

El mié, 1 sept 2021 a las 20:11, Daniel Arndt ()
escribió:

> Hermes,
>
> The SparseDirectUMFPACK solver doesn't work with MPI parallel matrices.
> You will need to find another solver based on the matrix class you are
> using. Standard choices would be GMRES or CG but I am assuming that your
> linear system might also be non-symmetric or indefinite.
>
> Best,
> Daniel
>
>
>
> Am Mi., 1. Sept. 2021 um 11:53 Uhr schrieb Hermes Sampedro <
> hermesampe...@gmail.com>:
>
>> Dear all,
>>
>> I am trying to solve a similar problem as step-29 but including the
>> massive parallel solution as presented in step-40.
>> I would like to ask how the solve() function would be in this case.
>>
>> In step-29:
>> //LU decomposition and inverse matrix
>> SparseDirectUMFPACK A_direct;
>> A_direct.initialize(system_matrix);
>>
>> //Solve with inverse matrix
>> A_direct.vmult(solution, system_rhs);
>>
>> When applying massive parallel computations, my initial guess was:
>>
>> LA::MPI::Vector
>> completely_distributed_solution(locally_owned_dofs,mpi_communicator);
>> //LU decomposition and inverse matrix
>> SparseDirectUMFPACK A_direct;
>> A_direct.initialize(system_matrix);
>>
>> //Solve with inverse matrix
>> A_direct.vmult(completely_distributed_solution, system_rhs);
>>  locally_relevant_solution = completely_distributed_solution;
>>
>>
>>  However, I get a *no matching member function for call to 'vmult'*
>> error.
>>
>> How can I use this LU decomposition and inverse matrix solver using MPI?
>>
>>
>> Thank you,
>>
>> Regards,
>>
>> H
>>
>> --
>> The deal.II project is located at http://www.dealii.org/
>> For mailing list/forum options, see
>> https://groups.google.com/d/forum/dealii?hl=en
>> ---
>> You received this message because you are subscribed to the Google Groups
>> "deal.II User Group" group.
>> To unsubscribe from this group and stop receiving emails from it, send an
>> email to dealii+unsubscr...@googlegroups.com.
>> To view this discussion on the web visit
>> https://groups.google.com/d/msgid/dealii/ded2e156-7bd7-4474-aa5f-fb57e8540926n%40googlegroups.com
>> 
>> .
>>
> --
> The deal.II project is located at http://www.dealii.org/
> For mailing list/forum options, see
> https://groups.google.com/d/forum/dealii?hl=en
> ---
> You received this message because you are subscribed to a topic in the
> Google Groups "deal.II User Group" group.
> To unsubscribe from this topic, visit
> https://groups.google.com/d/topic/dealii/ul4Aa_5t3n0/unsubscribe.
> To unsubscribe from this group and all its topics, send an email to
> dealii+unsubscr...@googlegroups.com.
> To view this discussion on the web visit
> https://groups.google.com/d/msgid/dealii/CAOYDWbKG87reQPBrhZu0RwoZHrMXMtVWS9Ui6EoFCWtmemG1vQ%40mail.gmail.com
> 
> .
>

-- 
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see 
https://groups.google.com/d/forum/dealii?hl=en
--- 
You received this message because you are subscribed to the Google Groups 
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/CAB%3DnHhYgcGWSv7aNs-9XzUL1C%3DBVeSFSecCEK1ZXocG-3qP%2B2A%40mail.gmail.com.

Re: [deal.II] Massively parallel & inverse matrix solver

2021-09-01 Thread Daniel Arndt 
Hermes,

The SparseDirectUMFPACK solver doesn't work with MPI parallel matrices. You
will need to find another solver based on the matrix class you are using.
Standard choices would be GMRES or CG but I am assuming that your linear
system might also be non-symmetric or indefinite.

Best,
Daniel



Am Mi., 1. Sept. 2021 um 11:53 Uhr schrieb Hermes Sampedro <
hermesampe...@gmail.com>:

> Dear all,
>
> I am trying to solve a similar problem as step-29 but including the
> massive parallel solution as presented in step-40.
> I would like to ask how the solve() function would be in this case.
>
> In step-29:
> //LU decomposition and inverse matrix
> SparseDirectUMFPACK A_direct;
> A_direct.initialize(system_matrix);
>
> //Solve with inverse matrix
> A_direct.vmult(solution, system_rhs);
>
> When applying massive parallel computations, my initial guess was:
>
> LA::MPI::Vector
> completely_distributed_solution(locally_owned_dofs,mpi_communicator);
> //LU decomposition and inverse matrix
> SparseDirectUMFPACK A_direct;
> A_direct.initialize(system_matrix);
>
> //Solve with inverse matrix
> A_direct.vmult(completely_distributed_solution, system_rhs);
>  locally_relevant_solution = completely_distributed_solution;
>
>
>  However, I get a *no matching member function for call to 'vmult'*
> error.
>
> How can I use this LU decomposition and inverse matrix solver using MPI?
>
>
> Thank you,
>
> Regards,
>
> H
>
> --
> The deal.II project is located at http://www.dealii.org/
> For mailing list/forum options, see
> https://groups.google.com/d/forum/dealii?hl=en
> ---
> You received this message because you are subscribed to the Google Groups
> "deal.II User Group" group.
> To unsubscribe from this group and stop receiving emails from it, send an
> email to dealii+unsubscr...@googlegroups.com.
> To view this discussion on the web visit
> https://groups.google.com/d/msgid/dealii/ded2e156-7bd7-4474-aa5f-fb57e8540926n%40googlegroups.com
> 
> .
>

-- 
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see 
https://groups.google.com/d/forum/dealii?hl=en
--- 
You received this message because you are subscribed to the Google Groups 
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/CAOYDWbKG87reQPBrhZu0RwoZHrMXMtVWS9Ui6EoFCWtmemG1vQ%40mail.gmail.com.

[deal.II] Massively parallel & inverse matrix solver

2021-09-01 Thread Hermes Sampedro 
Dear all,

I am trying to solve a similar problem as step-29 but including the massive 
parallel solution as presented in step-40.
I would like to ask how the solve() function would be in this case.

In step-29:
//LU decomposition and inverse matrix
SparseDirectUMFPACK A_direct;
A_direct.initialize(system_matrix);

//Solve with inverse matrix
A_direct.vmult(solution, system_rhs);

When applying massive parallel computations, my initial guess was:

LA::MPI::Vector
completely_distributed_solution(locally_owned_dofs,mpi_communicator);
//LU decomposition and inverse matrix
SparseDirectUMFPACK A_direct;
A_direct.initialize(system_matrix);

//Solve with inverse matrix
A_direct.vmult(completely_distributed_solution, system_rhs);
 locally_relevant_solution = completely_distributed_solution;


 However, I get a *no matching member function for call to 'vmult'* error. 

How can I use this LU decomposition and inverse matrix solver using MPI?


Thank you,

Regards,

H

-- 
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see 
https://groups.google.com/d/forum/dealii?hl=en
--- 
You received this message because you are subscribed to the Google Groups 
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/ded2e156-7bd7-4474-aa5f-fb57e8540926n%40googlegroups.com.