By the way, the solving procedure is done with the help of the direct
solver - SparseDirectUMFPACK
<https://www.dealii.org/current/doxygen/deal.II/classSparseDirectUMFPACK.html>
.
template <int dim>
void Poroelasticity<dim>::solve()
{
std::cout << "Solving linear system..";
SparseDirectUMFPACK A_direct;
A_direct.initialize(system_matrix);
A_direct.solve(system_matrix, system_rhs);
solution = system_rhs;
time.set_desired_next_step_size(1);
}
четверг, 14 октября 2021 г. в 11:15:44 UTC+2, Мария Бронзова:
> Dear Wolfgang,
>
> Thank you for the valuable hints! I did calculate a very simple system, I
> applied complex pressures (250,0) on certain boundaries and that is what I
> get after mapping and after using "apply_boundary_values":
>
> for (auto it = boundary_values.cbegin(); it!= boundary_values.cend(); ++it)
> {
> std::cout<<it->first<<" "<<it->second<<std::endl;
> }
>
> ------->
> 1 (0,0)
> 2 (0,0)
> 3 (250,0)
> 4 (0,0)
> 5 (0,0)
> 6 (0,0)
> 7 (0,0)
> 8 (250,0)
>
> Which makes perfect sense. But when I print out the solution values with
> the help of the DataOut for my pressures, I get the following values:
> (0,0)
> (0,0)
> (250,-3.92666e-06)
> (0,0)
> (0,0)
> (0,0)
> (0,0)
> (250,-5.66168)
>
> I am not quiet sure, what the problem might be. Would you possibly have an
> idea?
>
> Thank you a lot!
>
> Kind regards,
> Mariia
>
> среда, 13 октября 2021 г. в 17:35:12 UTC+2, Wolfgang Bangerth:
>
>> On 10/12/21 7:31 AM, Мария Бронзова wrote:
>> >
>> > The problem is that the real part gets perfectly assigned, whereas the
>> > imaginary value is ignored completely, even though I was trying to
>> > prescribe a zero value. In the results of calculation i get everything
>> > possible, rather than zero.
>>
>> Maria,
>> then it's time to debug :-) After the call to
>> interpolate_boundary_values(), are the values in the 'boundary_values'
>> map correct? After the call to apply_boundary_values(), are the values
>> in the 'solution' vector correct? (You can use DataOut to output the
>> solution vector at this point to look at it graphically.) What are the
>> values in 'solution' after solving?
>>
>> If you can't figure out what is wrong, can you shrink your program to
>> the minimal necessary to illustrate what you are seeing?
>>
>> Best
>> W.
>>
>>
>> --
>> ------------------------------------------------------------------------
>> Wolfgang Bangerth email: [email protected]
>> www: http://www.math.colostate.edu/~bangerth/
>>
>
--
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 [email protected].
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/4dc47872-b8fe-4e3f-ac34-0be5cc3b268an%40googlegroups.com.
