Handling vector valued
problems module? I was wondering when we might need to follow step-44 and
for example evaluate i_group and j_group s within the assembly function?
Thanks again,
Ali
On Wed, May 25, 2022 at 2:46 AM Wolfgang Bangerth
wrote:
> On 5/23/22 17:52, Ali Seddiq wrote:
> >
>
Dear all,
I am trying to use step 20 as a sample for assembling my linear system
which is composed of a 2*2 system matrix , a solution vector of two
scalars, and a right hand side vector of F and G. However, F in my case is
composed of two terms. Additionally in the system matrix which is define
Dear Wolfgang,
Thank you for your reply. In order to reconstruct the following final
discretized form as
(M+A)u + A'p=F+G,
a 2*2 matrix is subdivided into blocks in a way that (block)components
(0,0) and (1,1) are M+A and A' respectively and components (0,1) and (1,0)
set to zero.
Solution ve
Dear all,
I have a problem with imposing homogeneous Dirichlet boundary conditions
for pressure in a vector-valued problem.
I have defined interpolate_ and apply_boundary_values as following in the
assembly routine:
std::map boundary_values;
VectorTools::interpolate_boundary_values(dof_handler,
/11/22 09:54, Ali Seddiq wrote:
> >
> > Thanks very much for your advice. I definitely should and will follow
> that.
> > But may I still narrow my question for a clarification (with upcoming
> > complexity), and as a quick question ask if adding the boundary term
> thr
." to the cell_rhs for this purpose?
Thank you,
Ali
On Tue, Feb 8, 2022 at 10:19 PM Wolfgang Bangerth
wrote:
> On 2/8/22 09:12, Ali Seddiq wrote:
> > This is while if I change the pressure, naming it as U2 to another
> > scalar field variable , naming it as U1, it
dition of boundary term to cell_rhs leads to
the correct recognition as system_rhs.block(1) in the solve?
Would you or anyone else please help me with this problem?
Thanks very much,
Ali
On Monday, January 31, 2022 at 8:08:18 PM UTC+1 Wolfgang Bangerth wrote:
> On 1/31/22 11:41,
I guess. So not sure which
approach should be taken.
Thank you,
Ali
On Mon, Jan 31, 2022 at 6:29 PM Wolfgang Bangerth
wrote:
> On 1/31/22 10:22, Ali Seddiq wrote:
> >
> > I am trying to mimic step-7 implementation of Neumann boundary condition
> in
> > the case of
Hi everyone,
I am trying to mimic step-7 implementation of Neumann boundary condition in
the case of a vector-valued problem.
I did some modifications as shown in the below, but it doesn't work as
expected;
const FEValuesExtractors::Scalar pressure(0);
.
cell_rhs(i) +=
(
Dear Prof. Bangerth,
Thank you very much for your clear explanations.
So according to that, is it still correct if I sum up the equations in
discretized form? (Considering the fact that they're still containing trial
and test functions, though in discretized form).
Also could you please point me ou
Hello,
As my first question (which is an elementary math level question) and
considering step-20 as an example, I was wondering how two equations
(bilinear forms) containing two unknowns can be summed up to a single
equation? While for example in the case of a system of two linear equations
it is
Hello,
I'm very new to Deal.ii and trying to understand it through modifying
step-15. However I get the following error. Although it is descriptive I'm
unable to figure out the problem:
"error: no match for ‘operator*’ (operand types are ‘const
dealii::Tensor<1, 2>’ and
‘__gnu_cxx::__alloc_
12 matches
Mail list logo
