Hi everyone,
I am trying to implement an adaptive mesh refinement scheme for the Step-35
examle program.
I implemented the mesh-refinement acording to the Step-26 tutorial.
So I proceeded in the refine_mesh function as follows:
1. Create a Solution_Transfer object for each Vector;
2. prepare
Am Dienstag, 2. April 2019 01:48:41 UTC+2 schrieb Daniel Arndt:
>
> Alex,
>
> Hey I just added my changes to a github branch (
>> https://github.com/gabpeters/dealii/tree/sparsity )
>> I changed the make_sparsity_pattern(dof_row,dor_col,sparsity) function
>> and also changed the header in the
Hey, I just worked through the make_sparsity_pattern code.
I think I know how to fix the mentioned problem. But I dont understand the
how to use the "bool_dof_mask" in the correct way. It is used in he constraints.add_entries_local_to_global(...)
function.
In the make_sparsity_pattern
Hey everyone,
I am trying to apply boundary conditions on a Sparsematrix with two
Dof_handlers in a parallel::distributed::trinagulation setting for the
Step-35 tutorial.
( Same context as in this thread
https://groups.google.com/forum/#!topic/dealii/9ibgrQ0mFBs).
At initializing the
Hi. Thanks a lot for the advice.
I will Try it tomorrow and send a response whether it works.
Regards
Gabriel
Gabriel Peters
Endenicher Str. 310
53121 Bonn
00491525/5478185
gabriel.pet...@koeln.de
Am 26.03.19 um 13:32 schrieb Daniel Arndt
> Gabriel,
>
>
> > I am trying to
Hey Daniel,
thanks for the pull, with this instantiation everything compiles fine and I can use it.
But I think I did a mistake somewhere in the function. Using the definiion of the function from your pull request
I call
DoFTools::make_sparsity_pattern
Gabriel Peters
Endenicher Str. 310
53121 Bonn
00491525/5478185
gabriel.pet...@koeln.de
Am 05.04.19 um 15:55 schrieb Daniel Arndt
> Gabriel,
>
> to extend on my last post. Finally, I think that your approach with the two
> DoFHandler objects will give you some more prob
Am Mittwoch, 27. März 2019 13:58:47 UTC+1 schrieb Daniel Arndt:
>
> Gabriel,
>
> Am Mittwoch, 27. März 2019 06:16:18 UTC-4 schrieb Gabriel Peters:
>>
>> Hey, I just worked through the make_sparsity_pattern code.
>> I think I know how to fix the mentioned
Okay
Thanks a lot.
If it works now, I will contribute it.
Best
Gabriel
Gabriel Peters
Endenicher Str. 310
53121 Bonn
00491525/5478185
gabriel.pet...@koeln.de
Am 27.03.19 um 13:58 schrieb Daniel Arndt
> Gabriel,
>
> Am Mittwoch, 27. März 2019 06:16:18 UTC-4 schrieb Gabri
Am Donnerstag, 28. März 2019 21:37:19 UTC+1 schrieb Daniel Arndt:
>
> Gabriel,
>
> At first I extended the make_sparsity_pattern function to
>>
>> void DoFTools::make_sparsity_pattern ( const DoFHandlerType & dof_row,
>>const
>>
Am Montag, 18. Februar 2019 22:21:06 UTC+1 schrieb David Wells:
>
> Hi Gabriel,
>
> I believe that you are still using the standard deal.II SparseMatrix
> class; is this correct? If so, this class won't work with fully distributed
> calculations. You will need to use either the Trilinos or
Hey everybody,
I have a problem with initializing a Sparsematrix in a parallel program.
am using a parallel::distributed::triangulaton (as in the Step40 tutorial)
on
the Navier Stokes Projection solver from the step35 tutorial.
The problem arises in the function
Am Montag, 28. Januar 2019 15:19:36 UTC+1 schrieb Wolfgang Bangerth:
>
> On 1/28/19 6:38 AM, gabriel...@koeln.de wrote:
> >
> > But the problem is, that after one mesh-refinement step the time-step
> solving
> > of the Navierstokesprojection (in Step-35 it is in diffusion_step())
> >
I tried some things, but also if I insert a safety factor of 0.1 or less
the GMRES-procedure in the diffusion-step diverges.
I implemented the time step function such that
dt <= max_{T} h_T/(u_max(T))
where the max is taken over all cells and u_max is the maximal value of u
on the cell T.
Hey everyone,
I am wondering, whether there is a deal.ii function, which balances a
parallel::distributed::triangulation neighbouring elements by 2:1,
like the p8est_balance_ext function
(see
http://p4est.github.io/api/p8est__extended_8h.html#ad3a61dc2b4b5cc69b03c995f188bfd8d
)
Such a
: "Wolfgang Bangerth"
An: "dealii@googlegroups.com"
Betreff: Re: [deal.II] p4est balance_mesh
On 6/4/19 6:41 AM, Gabriel Peters wrote:
>
> I am wondering, whether there is a deal.ii function, which balances a
> parallel::distributed::triangulation neighbo
Ok thanks a lot.
Best
Gabriel
Gabriel Peters
Endenicher Str. 310
53121 Bonn
00491525/5478185
gabriel.pet...@koeln.de
Am 05.06.19 um 14:55 schrieb Wolfgang Bangerth
> On 6/5/19 3:36 AM, Gabriel Peters wrote:
> > I am setting refinement_flags on cells, based on the geometry I u
Hey everyone,
I am working with a parallel::distributed::triangulation.
In my program I load a coarse grid from a file.
Is there a function to manually initialize trees in my triangulation,
according to some regions of my domain (just as p4est does) or
is this automatically done by calling
Hey everyone,
I have a problem at solving a blocksystem.
I have a matrix called "system_matrix" in Blockform M = [[A,0],[B,C]]
and I want to solve classical Mx = F. (F = [f,0], x = [x1,x2])
I proceed according to the LinearSolvers class in the step-32 tutorial.
so by blockwise solving the
Hey Jean-Paul
Thanks a lot, Now it works finely.
Best regards
Gabriel
Gabriel Peters
Endenicher Str. 310
53121 Bonn
00491525/5478185
gabriel.pet...@koeln.de
Am 20.05.19 um 15:01 schrieb Jean-Paul Pelteret
> Dear Gabriel,
>
> It looks to me like you forgot to dereference the share
Hey everyone,
I have a problem with adding Vectors in a Multiparallel program.
Here I use TrilinosWrappers::MPI::BlockVector vector_1, vector_2;
which are initialized ion the same fashion, as in the step-32 tutorial
vector_1.reinit(stokes_relevant_partitioning, MPI_COMM_WORLD);
Hey everyone,
I got a problem regarding TrilnosWrappers::BlockSparseMatrix.
In my program I have a TrilnosWrappers::BlockSparseMatrix system_matrix
which consists of two
parts, one permanent part, which is computed at the beginning at the
programm, not every iteration, to save time (called
Hey Jean-Paul
Thanks a lot, that's what I was looking for, and now it works finely.
Best
Gabriel
Gesendet: Montag, 06. Mai 2019 um 14:27 Uhr
Von: "Jean-Paul Pelteret"
An: dealii@googlegroups.com
Betreff: Re: [deal.II] Adding TrilinosWrapper::BlockSparseMatrix
Dear Gabriel,
There
Hey everyone,
I am setting up a system matrix in Trilinos::Blocksparsematrix with a
DOfHandler, which has 3 velocity components and 1 pressure component
(analogously to step-32).
Somehow everything worked fine a long time, but now I get an error I dont
understand.
I call:
*std::vector
Hey everyone,
I have an error at compiling deal_ii_with_Trilinos=ON, which I dont
understand.
First I installed,then I called
cmake -DCMAKE_INSTALL_PREFIIX=/home/peters/Documents/installs
-DDEAL_II_WITH_TRILINOS=ON
-DTRILINOS_DIR=/home/peters/Documents/build_trilinos
This all worked finely,
Hey everyone,
I have an error at compiling deal_ii_with_Trilinos=ON, which I dont
understand.
First I installed Trilinos,then I called
cmake -DCMAKE_INSTALL_PREFIIX=/home/peters/Documents/installs
-DDEAL_II_WITH_TRILINOS=ON
-DTRILINOS_DIR=/home/peters/Documents/build_trilinos
This all
Hey everyone,
I want to set up a BlockMatrix, with entries, which are enforced by a
velocity components of aBlockVector "solution_old".
The Vector contains velocity entries in the first dim terms and pressure
entries in the last component.
Similar to the 'assemble_system' function of Step-21
Gabriel Peters
Endenicher Str. 310
53121 Bonn
00491525/5478185
gabriel.pet...@koeln.de
Am 25.04.19 um 15:28 schrieb Wolfgang Bangerth
> On 4/25/19 4:51 AM, gabriel.pet...@koeln.de wrote:
> >
> > std::vector > solution_old_values(n_q_points);
> >
> > While
Hey everyone,
I have the following problem
I have a classical velocity - pressure Blockvector "soluion" and a FEsystem
(initialized via
* fe(FE_Q(degree + 1), dim, FE_Q(degree), 1)*
For some testing I want to get the velocity compontents of solution in
seperate vectors,
i.e. a
Am Donnerstag, 18. Juli 2019 10:52:46 UTC+2 schrieb Wolfgang Bangerth:
>
> On 7/9/19 9:01 AM, Gabriel Peters wrote:
> >
> > After some debugging I found out that the term /(-1)*
> sol_p[i]*div_sol_u[j]/
> > makes this trouble. Can somebody tell me why?
>
&g
30 matches
Mail list logo