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
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
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
: "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
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
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 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 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 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 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 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
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 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
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 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,
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
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
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
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 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
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
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 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
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
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.
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())
> >
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
30 matches
Mail list logo