Hello Daniel,
thank you for the answer.
I have read in the documention only
This function is called by some GridIn::read_* functions.
>
fully forgetting the following
Only the vertices with indices in considered_vertices are tested for
> equality. This speeds up the algorithm, which is quadratic
Hello Daniel,
I have fixed the problem using Your proposed algorithm of mesh creation.
The strange behaviour of velocity-pressure field was really caused by the
bad quality of the mesh (possibly doubled vertices at the boundary between
merged entities).
I want to ask a question, which could be inte
016 02:18 PM, Vinetou Incucuna wrote:
>
>>
>> I have tried to output the solutions from the step one as
>> well as from the step three. There is some development, but
>> reasonable,i.e. solution from the first step and the third step look
>> similar (except the fir
Sorry, I wanted to say, that the imposition of boundary
conditions in the setup is apparently dominant.
Is this the general case in your applications?
Thank You
Marek
2016-07-30 13:09 GMT+02:00 Vinetou Incucuna :
> ah, sorry for interruption.
> From the log it is evident, than the
:00 Vinetou Incucuna :
> Hello,
> problem was,
>
>- that I have some meaningless calls for
>
> /*current_solution_phase.reinit(locally_owned_dofs_phase,
> locally_relevant_dofs_phase, mpi_communicator);
>
Hello,
problem was,
- that I have some meaningless calls for
/*current_solution_phase.reinit(locally_owned_dofs_phase,
locally_relevant_dofs_phase, mpi_communicator);
current_solution_phase = 0;
*/
/*solution_update.reinit(locally_owned_dofs_phase,
Hello,
I believe I solved partially the problem, do
not bother to answer the thread without my further
specification.
M
2016-07-29 10:17 GMT+02:00 Vinetou Incucuna :
> Hello,
> thank You for response and the advices.
>
>
> >> LA::Vector vec_old_solution(d
Hello,
thank You for response and the advices.
>> LA::Vector vec_old_solution(dof_handler_nse.n_dofs());
>>
>> VectorTools::interpolate(dof_handler_nse, ZeroFunction<3>(dim +
1),
>> vec_old_solution);
>>
>> old_solution_nse = vec_old_solution;
>
> > Are yo
Matthias, just to make things clears. What everything
should be build in release mode - deal.ii, trilinos (and
dependency libraries Blas, Lapack)
, SuperLU_Dist ?
I will give a try on Dofs renumbering.
Thank You
M
2016-07-22 13:43 GMT+02:00 Jean-Paul Pelteret :
> Also, have you tried optimisin
ctive design of the distributed
computation like the one in Step-40.
Any advice will be appreciated
Thank You
Marek
2016-07-22 8:42 GMT+02:00 Vinetou Incucuna :
> Hello,
> once more again
>
>>
>> 71874 dofs for Cahn-Hilliard part of system, cca 2000-3000 dofs per computer
>
Turcksin :
>
> Marek,
>
> 2016-07-21 15:59 GMT-04:00 Vinetou Incucuna :
> > I suppose, that the next
> > step will be an solver for general systems,
> > like
> > TrilinosWrappers::SolverBicgstab or TrilinosWrappers::SolverGMRES
> > ?
> Yes that'
Hello,
thank you for the answer.
I suppose, that the next
step will be an solver for general systems,
like
TrilinosWrappers::SolverBicgstab or TrilinosWrappers::SolverGMRES
?
M
2016-07-21 21:39 GMT+02:00 Bruno Turcksin :
> Marek,
>
> On Thursday, July 21, 2016 at 3:17:48 PM UTC-4, Marek Čap
>
> One more thing. Because the domain is curved, you will need to attach
> a manifold to the Triangulation
Thank You Bruno for the warning, I have never used it (I only assigned
boundary objects
to mesh). This feature is assumably more general
2016-07-15 17:45 GMT+02:00 Vinetou
are in dealii
some helpful features, which diminish this burden?
2016-07-15 17:37 GMT+02:00 Bruno Turcksin :
> Marek
>
> 2016-07-15 11:27 GMT-04:00 Vinetou Incucuna :
> > i have redesigned my *.geo files and got
> > much coarser meshes, which work in dealii GridIn
> Grea
14 matches
Mail list logo
