Dear Jean-Paul,
I have checked Step-40. Lets say we have an* old_locally_relevant_solution
*vector
and a *new_locally_relevant_solution_update* vector (*extracted from the
completely distributed solution in the solve module*). How can we add these
two vectors in DEALII. Thanking you in anticipation.
On Thursday, February 21, 2019 at 3:32:08 PM UTC+5:30, mboro...@gmail.com
wrote:
>
> Dear all,
> I am Boron. I am a new to DEALII. I am currently trying to write a
> parallel code in DEALII for solving nonlinear Poisson's equation. The file
> is also attahed below. My doubt is "How do we pass history variable while
> constructing the cell_matrix?"
> A code snippet is (Line No 211-225) :
>
> for (; cell!=endc; ++cell)
> {
> if (cell->subdomain_id() == this_mpi_process)
> {
> cell_matrix = 0;
> cell_rhs = 0;
>
> fe_values.reinit (cell);
>
> fe_values.get_function_values(present_solution,
> old_solution);
> fe_values.get_function_gradients(present_solution,
> old_solution_gradients);
> for (unsigned int q_point=0; q_point<n_q_points; ++q_point)
> {
> // BUILD ELEMENTAL CELL MATRIX @ EACH GAUSS POINT
> }
>
> In the above code snippet, the line
> 'fe_values.get_function_values(present_solution, old_solution); ' throws
> an error. Is there a way to pass only the vectors relevant to the
> corresponding subdomains in DEALII?
>
>
--
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.
For more options, visit https://groups.google.com/d/optout.
