On 4/25/19 2:50 PM, Gary Uppal wrote:
>
> I am trying to solve the diffusion equation with Discontinuous Galerkin
> elements. The solution looks good with a regular structured mesh, but if I
> distort the mesh, the solution blows up and does not converge. Is there an
> obvious reason this
On 4/25/19 10:13 AM, 'Maxi Miller' via deal.II User Group wrote:
> After running some tests, I ended up reducing the problem to the transfer to
> and from the Epetra-vectors. I have to write an interface to the model
> (according to the instructions), and as shown in the code above. There I have
Dear Wolfgang, this was easier than I thought.
I followed your suggestion and used VectorTools::interpolate.
Everything works as expected.
Thank you very much
Final solution:
// initialConditionParameters.uvw is the parsed function
// newton_update is a local vector
const
Dear Wolfgang,
This was easier than I thought. I followed your suggestion and used
VectorTools::interpolate
Final solution to set the values of the DOF using a parsed function:
// initialConditionParameters.uvw is my parsed function
// newtonUpdate is my local vector that I use to alter the DOF
On 4/25/19 9:46 AM, Muhammad Mashhood wrote:
>
> As unlike the step-8 here in step-26 the rhs is being formed directly
> globally as "system_rhs" using "create_right_hand_side()" function so
> can I use the "create_boundary_right_hand_side()" function using
> following lines of code?
It is
After running some tests, I ended up reducing the problem to the transfer
to and from the Epetra-vectors. I have to write an interface to the model
(according to the instructions), and as shown in the code above. There I
have Epetra-Vectors created by my interface, with a size of
In the weak form after Neuman BC implementation for equation in step-26 i
got following additional terms contributing in rhs:
*for 1D domain x = [0, 1]*
*system_rhs + k_n * [ g1 * phi_(x1) - g0 * phi_(x0) ], where g0 and g1 are
the flux values at boundary nodes and phi_(x0), phi_(x1) are shape
On 4/25/19 7:56 AM, Gabriel Peters wrote:
> I wanted solution_o to be vector-valued.
> And I want to multiply it with a gradient function. So I want
>
> local_matrix(i,j) += solution_o * fe_values[velocities].shape_gradient(i,q)
> *
> fe_values[velocities].shape_gradient(j,q);
>
> Does this
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 assembling the system over all
On 4/25/19 6:58 AM, Bruno Blais wrote:
> I would like to be able to set-up an initial condition by fixing the value of
> the DOF themselves by using their x,y(,z) position with a ParsedFuction.
>
> Generally, I set the initial condition by using an L2 projection of the
> function (for complex
On 4/25/19 4:51 AM, gabriel.pet...@koeln.de wrote:
>
> std::vector > solution_old_values(n_q_points);
>
> While assembling the system over all cells I call
>
> fe_values.get_function_values (solution_old, solution_old_values);
>
> In the Step-21 tutorial, they call
>
> const double old_s
On 4/25/19 3:20 AM, MD wrote:
> I am trying to solve an nonlinear laplace problem with only dirichlet
> boundaries using dealii with distributed triangulation. I do have a solution
> vectorlocally_owned_sol and a newton_update. At the Dirichlet boundaries, I
> want torestrict the newton update
Hello everyone,
I believe this is a relatively dumb question, but I seem to struggle with
this basic concept.
I would like to be able to set-up an initial condition by fixing the value
of the DOF themselves by using their x,y(,z) position with a ParsedFuction.
Generally, I set the initial
It looks as if you have to recompile Trilinos using the -fPIC-switch to
create position independent code.
Am Donnerstag, 25. April 2019 12:10:52 UTC+2 schrieb gabrie...@koeln.de:
>
> Hey everyone,
> I have an error at compiling deal_ii_with_Trilinos=ON, which I dont
> understand.
>
> First I
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,
Hello all,
I am trying to solve an nonlinear laplace problem with only dirichlet
boundaries using dealii with distributed triangulation. I do have a
solution vector locally_owned_sol and a newton_update. At the Dirichlet
boundaries, I want to restrict the newton update to be zero and the
18 matches
Mail list logo