Re: [deal.II] modifying of only one component of solution

2018-02-07 Thread Marek Čapek
Hi, the source code of the method void Main::postprocess_solution(){ was correct. The function truncate_conc(double conc ) had a faulty implementation. It was do-nothing implementation on argument. Greets Marek Dne pondělí 5. února 2018 19:58:23 UTC+1 Marek Čapek napsal(a): > > Hi, > I

Re: [deal.II] modifying of only one component of solution

2018-02-05 Thread Marek Čapek
Hi, I have Cahn-Hilliard system - scalar equation for phase field and chemical potential. I tried to modify only the phase field part of solution in distributed code like step-40: template > void Main::postprocess_solution(){ > > nsSystem.computing_timer.enter_subsection("Solution postprocess

Re: [deal.II] modifying of only one component of solution

2017-08-28 Thread Wolfgang Bangerth
On 08/27/2017 11:53 AM, Marek Čapek wrote: I am getting the solution in PetscWrappers::MPI::Vector solution_phase_n; I am interested in modifying the vector of the solution - I want to cut the undershoots and overshoots in phase-field component of the solution. Namely I have values of phas

[deal.II] modifying of only one component of solution

2017-08-27 Thread Marek Čapek
Hello, I solve phase-field equation, i.e. Cahn-Hilliard equation. I have the following finite elements for the phase field and the chemical potential components of solution FESystem fe_phase; in the initialization list of constructor: fe_phase (FE_Q (degree), 2) (The first component - pha