Re: [deal.II] Re: Read serial vector into a parallel vector

2016-06-20 Thread Daniel Arndt
Praveen, Am Montag, 20. Juni 2016 18:05:12 UTC+2 schrieb Praveen C: > > Hello Daniel > > I have a normal Triangulation and I solve on this using a Vector. > I could save this to file calling Vector::print. > > Is it now possible to read this into a TrilinosWrappers::MPI::Vector ? > No, there is n

[deal.II] Re: Reuse factorization in Trilinos direct solver

2016-06-20 Thread Denis Davydov
FYI, the doxygen field was wrong for the extra solve() function, thus it was not shown. Should be fixed in https://github.com/dealii/dealii/pull/2694 On Monday, June 20, 2016 at 3:48:56 PM UTC+2, Praveen C wrote: > > Hi > > I want to use the direct solver MUMPS via Trilinos. I need to solve the

[deal.II] Re: Reuse factorization in Trilinos direct solver

2016-06-20 Thread Denis Davydov
Hi Praven, I don't think you need to hack anything, this is AFAIK already possible, see https://github.com/dealii/dealii/pull/2493 You need to compile deal.II from Git to get this. Regards, Denis. On Monday, June 20, 2016 at 3:48:56 PM UTC+2, Praveen C wrote: > > Hi > > I want to use the direct

Re: [deal.II] WorkStream iterator

2016-06-20 Thread Daniel Shapero
So glad I asked this question instead of reinventing the wheel :) I was mostly working off the Boussinesq and elasticity examples, so I missed step-35. On Sunday, June 19, 2016 at 3:18:33 PM UTC-7, bangerth wrote: > > On 06/19/2016 04:12 PM, Daniel Shapero wrote: > > I'm working on a problem wh

Re: [deal.II] Re: Reuse factorization in Trilinos direct solver

2016-06-20 Thread 'Uwe Köcher' via deal . II User Group
Here you can find some class, which I'm using for that: https://github.com/koecher/meat/blob/master/source/DTM%2B%2B/core/lss/SolverAmesos.cc Best Uwe On Monday, June 20, 2016 at 6:22:40 PM UTC+2, Praveen C wrote: > > Thanks Uwe. I will try that. > > Best > praveen > > On Mon, Jun 20, 2016 at

Re: [deal.II] Re: Reuse factorization in Trilinos direct solver

2016-06-20 Thread Praveen C
Thanks Uwe. I will try that. Best praveen On Mon, Jun 20, 2016 at 7:54 PM, 'Uwe Köcher' via deal.II User Group < dealii@googlegroups.com> wrote: > No, this interface uses a solve function, in which the algebraic and > numeric factorisation is done each time the solve() function > is called. This

Re: [deal.II] Re: Read serial vector into a parallel vector

2016-06-20 Thread Praveen C
Hello Daniel I have a normal Triangulation and I solve on this using a Vector. I could save this to file calling Vector::print. Is it now possible to read this into a TrilinosWrappers::MPI::Vector ? If not, what other option do I have to achieve this result. Thanks praveen On Mon, Jun 20, 2016

[deal.II] Re: Reuse factorization in Trilinos direct solver

2016-06-20 Thread 'Uwe Köcher' via deal . II User Group
No, this interface uses a solve function, in which the algebraic and numeric factorisation is done each time the solve() function is called. This has something to do with the common interface for the direct and iterative solvers. But it is easy to implement yourself, take a look into lac/trilino

[deal.II] Re: Read serial vector into a parallel vector

2016-06-20 Thread Daniel Arndt
Praveen, I solve a PDE in serial and save the solution to file. > Do you mean that your triangulation is not a parallel::distributed::Triangulation or that you are simply running with one process? > Is there a way to now read this solution into a > TrilinosWrappers::MPI::Vector object ? > If

[deal.II] Read serial vector into a parallel vector

2016-06-20 Thread Praveen C
Hello I solve a PDE in serial and save the solution to file. Is there a way to now read this solution into a TrilinosWrappers::MPI::Vector object ? Thanks praveen -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/

[deal.II] Reuse factorization in Trilinos direct solver

2016-06-20 Thread Praveen C
Hi I want to use the direct solver MUMPS via Trilinos. I need to solve the same matrix problem repeatedly. How can I reuse the LU factorization ? In SparseDirectUMFPACK, calling vmult allows me to solve many problems with same factorization. Is such a thing possible with TrilinosWrappers::Solver