Re: [deal.II] Question about complex value finite element

2020-08-16 Thread yuesu jin
Sure, I will add my progress following this email. Best regards, On Sun, Aug 16, 2020 at 7:47 PM Wolfgang Bangerth wrote: > On 8/16/20 6:38 PM, yuesu jin wrote: > > Thank you very much! I found step-62 yesterday and I have > compared the > > difference between them. I want to write two

Re: [deal.II] PETSc iteration does not converge

2020-08-16 Thread yuesu jin
Dear Dr.Bangerth, I tried both, first I tried 1e-4*system_rhs.l2_norm(), it failed. Then I tried with exact precision limits, from 1e-3 to 1e-8, all of them cannot converge. If I set the precision level as 1e-4, with or without L2 norm it stopped at step-42 with the residue 0.0113395. Best

Re: [deal.II] Question about complex value finite element

2020-08-16 Thread Wolfgang Bangerth
On 8/16/20 6:38 PM, yuesu jin wrote:    Thank you very much! I found step-62 yesterday and I have compared the difference between them.  I want to write two versions and compare their CPU time to get an exact result. That would be interesting to hear about -- please let us know what you

Re: [deal.II] Question about complex value finite element

2020-08-16 Thread yuesu jin
Dear Dr. Bangerth, Thank you very much! I found step-62 yesterday and I have compared the difference between them. I want to write two versions and compare their CPU time to get an exact result. Best regards, Yuesu On Sun, Aug 16, 2020 at 7:19 PM Wolfgang Bangerth wrote: > Yuesu, > > >I

Re: [deal.II] PETSc iteration does not converge

2020-08-16 Thread Wolfgang Bangerth
On 8/14/20 9:20 PM, yuesu jin wrote: I tried a few levels of the tolerance, 1e-3, 1e-4 ,1e-6 and 1e-8. 1e-3 converged because that is in the same level with the solution. Error occurs below 1e-4.  I have tried the Jacobi preconditioner instead of SSOR in the single thread CG version, it

Re: [deal.II] Question about complex value finite element

2020-08-16 Thread Wolfgang Bangerth
Yuesu,   I am adding an absorbing boundary condition like step-24 did in the frequency domain, in which the time derivative gives a complex term.  I also found the complex acoustic wave problem step-29 which splits the complex wave function into two real parts.   What I want to know is what

Re: [deal.II] Re: Memory loss in system solver

2020-08-16 Thread Richard Schussnig
> > Hi again, it seems my last comment got lost, but I wanted to post it here for others researching that issue: I constructed a MWE and found out, that actually solving the system via GMRES, CG or BiCGstab using any of the implementations provided by PETSc, Trilinos or the dealii versions