On 10/26/2017 05:53 PM, Yiyang Zhang wrote:
|
IndexSetglobal_active_cells(this->tria.n_active_cells());
global_active_cells.add_range(0,this->tria.n_active_cells());
IndexSetlocal_active_cells(this->tria.n_active_cells());
local_active_cells.clear();
intcount =0;
for(autocell =this->t
Hello,
I am using parallel::shared::Triangulation and
TrilinosWrappers::MPI::Vector. It seems the refine_mesh procedure in Step
-18 cannot be directly transferred to Trilinos... Is there any tutorials on
this object? I appreciate if anyone can tell me some references...
Meanwhile, there is so
Hi,
I am a beginner user. This post is not a question; my intention is to
share possibly helpful information with other beginners like me.
I installed deal.ii via spack, because the README for deal.ii made the
manual installation look too hard and/or time-consuming, and I like conda
in Python, a
Hi,
Test 4 failed after make test. I make my research and find the following
answer:
https://groups.google.com/d/msg/dealii/mV-juFpnybA/DSBP4DArCQAJ
I guess it may be caused by different version of MPI. And I get the
following return of suggested command:
*
Hello Prof. Bangerth,
Thank you for the answer! Now things are clear to me.
Best,
Yiyang
On Thursday, October 26, 2017 at 11:00:39 AM UTC-5, Wolfgang Bangerth wrote:
>
> On 10/26/2017 09:50 AM, Yiyang Zhang wrote:
> >
> > It seems to me that both ways are fine. Is there any particular reason
Frederik,
without looking at your code (sorry!), I think you can use
Triangulation::register_data_attach and notify_ready_to_unpack to work
with any kind of data that you want to transfer around (this is the
functionality that is used by SolutionTransfer and the
ContinuousQuadratureDataTransfer.
Hello!
I've got a question on the usage of CellDataStorage for refined cells for a
parallelized simulation. Below you'll find a sample of the code I use for
refinement.
I'm transferring the solution as stated
in
https://www.dealii.org/8.5.0/doxygen/deal.II/classparallel_1_1distributed_1_1Sol
On 10/26/2017 09:50 AM, Yiyang Zhang wrote:
It seems to me that both ways are fine. Is there any particular reason
that we prefer one than the other?
Yes, both ways are correct.
But you will need a constraint object to eliminate hanging node
constraints and boundary values from the newton_
Hello Prof. Bangerth,
Thank you for your reply. I am sorry I didn't really follow...
Here are the two ways.
(i)
set the inhomogeneity in the initial solution vector.
in the constraints object we only need to store homogeneous constraints:
constraints.add_lines()
Then in the solve() func
On 10/26/2017 09:19 AM, 'Maxi Miller' via deal.II User Group wrote:
No, but including it removes the compilation errors, thanks!
How do I know for the future which includes I should use?
This one was difficult, but in general when you get a compiler error for
a function that the compiler can't
On 10/26/2017 09:16 AM, Yiyang Zhang wrote:
I mean in the constraints objects, I also put the inhomogeneous
constraints in, that is:
constraints.add_lines()
constraints.add_inhomogeneity()
Then, if I use version A, every time I solve the equation, I will add
the inhomogeneity to the
No, but including it removes the compilation errors, thanks!
How do I know for the future which includes I should use?
Am Donnerstag, 26. Oktober 2017 17:03:40 UTC+2 schrieb Wolfgang Bangerth:
>
> On 10/26/2017 07:27 AM, Wolfgang Bangerth wrote:
> >
> > This looks reasonable. The code does not
Hello Bruno,
Thank you for your quick reply!
ii) ConstraintMatrix::distribute() and Newton's method
> In a Newton's method solver, we have:
>
> version A:
> //...get the new newton_update vector from solving the equation
> constraints.distribute(newton_update);
> new_sln = old_sln + alpha* newto
On 10/26/2017 08:17 AM, 曾元圆 wrote:
Is there a record on how the template/class/function names are changed,
so I can find it quickly in case one day in the future I found the code
I wrote today can't run?
Yes, we list changes from each version to the next here:
https://www.dealii.org/develo
On 10/26/2017 07:27 AM, Wolfgang Bangerth wrote:
This looks reasonable. The code does not compile because of this error:
/opt/dealII/include/deal.II/base/template_constraints.h:388:40: error:
no match for 'operator*' (operand types are 'Sacado::Fad::DFad'
and 'dealii::Tensor<1, 2, double>')
Wow, it's really surprising and exciting to see the original developer of
the code I'm now runing answering this question! Thank you so much Prof.
Praveen.This is greetings from remote China! I'm now a master student and
am running and learning the DG code developed by you. So I really want to
Sorry, I should have describe my problem more specificlly. The problem is,
I'm running a code writen several years ago. In the code it used the class:
NamedData< Vector *>, but Deal.II has deleted this template class
so I have no idea how it work and the way to use it. I searched in the
mailin
On 10/26/2017 02:18 AM, 'Maxi Miller' via deal.II User Group wrote:
I added the errors I get while compiling in the log file. My code for
producing it is:
|
std::vector>>grad_N_AD(n_q_points);std::vector>local_dof_values_AD(cell->get_fe().dofs_per_cell);
for(unsignedinti=0;i
This looks reaso
Yiyang,
On Wednesday, October 25, 2017 at 8:07:46 PM UTC-4, Yiyang Zhang wrote:
> locally_relevant_sln_with_ghost_new = locally_relevant_sln_with_ghost_old;
> //This and the following one are rarely seen in the tutorials. Does this
> mean the vector is exactly copied? i.e. both locally_owned_do
step-39 should be of some help to you.
best
praveen
> On 26-Oct-2017, at 5:45 PM, 曾元圆 <2012zg...@gmail.com> wrote:
>
> I am trying to run a open source code, but an error occurred.
> In the code it uses NamedData like this:
> NamedData< Vector *> solution_data;
> Info_box.initialize(fe, mappin
On 10/26/2017 06:15 AM, 曾元圆 wrote:
I am trying to run a open source code, but an error occurred.
In the code it uses NamedData like this:
NamedData< Vector *> solution_data;
Info_box.initialize(fe, mapping, solution_data);
What is the error you get? The first step in fixing an error is to care
I am trying to run a open source code, but an error occurred.
In the code it uses NamedData like this:
NamedData< Vector *> solution_data;
Info_box.initialize(fe, mapping, solution_data);
How should I use AnyData to replace it?
Actually I am a new learner, I tried to read the documentation relat
22 matches
Mail list logo
