Hi Greg,
I guess what you could do is to treat the vector as vector associated with
a DoFHandler with FE_DGQ(0). During each refinement you "interpolate" the
vector to the new mesh and decrease the value by one. This approach
naturally extends to a parallel setting.
Alternately, you could crea
Dear all,
This might be a somewhat odd feature request but I was wondering whether
the functionality described below can be readily achieved in deal.II
already:
Say we are given a coarse mesh with every cell never refined, and a vector
of non-negative integers associated with each cell. The
Alright, I've had some decent success, though I've hit another snag in the
solver section. Your suggestions have opened the doors and actually
allowed me to shift gears to my preferred approach.
I've changed the end of Setup_System() to be:
dsp.collect_sizes();
*DoFTools::make_spars
On 11/8/23 02:44, Vinayak Vijay wrote:
Can someone help me with this? The main objective is to obtain the global dof
indices using the -dealii::TriaActiveIterator>
The error is generated in a place that says that you shouldn't be calling this
constructor. I don't see the backtrace, so can't
Dear Luca and Wolfgang,
Thanks for your pointers; I will work to implement them.
My github account is: https://github.com/AlexQuinlan
Thanks,
Alex
On Wednesday, November 8, 2023 at 3:36:43 AM UTC-5 luca@gmail.com wrote:
> Dear Alex,
>
> if you send me your github user, I’ll give you acce
Dear all,
attached is a minimal example in which I create two simple triangulations
(tria1, tria2),
merge them, and finally anisotropically refine all cells
(RefinementCase<3>::cut_xz).
I also write the merged triangulation to a file for visualizations sake.
My objective is to refine only in th
Hello,
I am trying to obtain the vertex global dofs using the DOFAccessor. I
employ the following code or the same:
const auto vertex_end = tria.end_vertex();
for (auto vertex = tria.begin_active_vertex(); vertex != vertex_end;
++vertex)
{
DoFAccessor<0, 1, 3, false> vertex_dofs(&tria,
vertex->
Dear Alex,
if you send me your github user, I’ll give you access to a code that does
exactly what you are trying to do, so that you can get some inspiration on how
to assemble the coupling matrices.
There is an open PR (the review has been stalling for some time on that)
https://github.com/de
