Re: [deal.II] RHS function represented by an interpolation

2022-07-15 Thread Jean-Paul Pelteret
Hi Alexander, The answer that Wolfgang gave is the one that you want to follow. But I can give some input into your one question. > The problem here is the step 3: somehow I need to specify the same geometric > cell for two FEValues with different finite element systems while preserving > the

Re: [deal.II] RHS function represented by an interpolation

2022-07-15 Thread Wolfgang Bangerth
On 7/15/22 08:36, Alexander Kiselyov wrote: Which tools could be used to overcome this problem? Or is my approach deficient in general? Alexander: You want to evaluate the solution obtained on one DoFHandler at quadrature points so that you can form the right hand side for a system that live

Re: [deal.II] Suggestions for finding local minima on distributed, adaptively-refined meshes

2022-07-15 Thread Wolfgang Bangerth
Lucas: I suspect this is a simple task, but I can't find anything which does it outright. I'd like to find the local minima associated with some quantity derived from a finite element field. The features that I'm looking at are reasonably well-behaved, in that they're smooth except exactly a

Re: [deal.II] How to use iterative solver to solve complex value problem?

2022-07-15 Thread Wolfgang Bangerth
On 7/15/22 07:35, 沈键 wrote: **   In step58, it suggested to use iterative solver to improve performance. And I tried gmres, but program threw error: "solver_gmres.h:676:13: error: cannot convert ‘std::complex’ to ‘double’ in assignment". My code is as follows: void TEST_Solver_Complex::solv

Re: [deal.II] Trilinos::Amesos_Superludist not scaling across multiple nodes

2022-07-15 Thread Wolfgang Bangerth
I'm not sure any of us have experience with Amesos:SuperLU, so I'm not sure anyone will know right away what the problem may be. I was wondering if, while writing the wrappers and testing them out, someone managed to figure out the requisite combination of installation time options.

[deal.II] Suggestions for finding local minima on distributed, adaptively-refined meshes

2022-07-15 Thread Lucas Myers
Hi folks, I suspect this is a simple task, but I can't find anything which does it outright. I'd like to find the local minima associated with some quantity derived from a finite element field. The features that I'm looking at are reasonably well-behaved, in that they're smooth except exactly a

[deal.II] RHS function represented by an interpolation

2022-07-15 Thread Alexander Kiselyov
Dear deal.II users, I'd like to ask a question concerning the use of RHS function represented by a DoF vector from another solver. I'm trying to build RHS of one solver (A) using solution obtained with another solver (B) , i.e. instead of Function there is a DoF vector. Both solvers have the same

[deal.II] How to use iterative solver to solve complex value problem?

2022-07-15 Thread 沈键
In step58, it suggested to use iterative solver to improve performance. And I tried gmres, but program threw error: "solver_gmres.h:676:13: error: cannot convert ‘std::complex’ to ‘double’ in assignment". My code is as follows: void TEST_Solver_Complex::solve() { LogStream::Prefix p("Solve"

Re: [deal.II] Trilinos::Amesos_Superludist not scaling across multiple nodes

2022-07-15 Thread Paras Kumar
Dear Wolfgang, Thank you for the response. > Paras: > I'm not sure any of us have experience with Amesos:SuperLU, so I'm not > sure > anyone will know right away what the problem may be. > > I was wondering if, while writing the wrappers and testing them out, someone managed to figure out the re