Hello all,
This code just solves the wave equation. Despite the fact that I have
provided an initial velocity in line 106, it does not appear to be reading
it in lines 1024-1062. My solution looks like the wave equation with 0
initial velocity, i.e. one hump turns into two moving away from each
Hi,
Thank you for taking a look. You are correct in that if I comment
out VectorTools::interpolate_boundary_values in line 320, the solution does
not change. If I understand you correctly, you are saying that the function
apply_boundary_values defined in line 553 and called in line 630 must tak
I've recently implemented a 1d standard wave equation and have been messing
around with different boundary conditions. For instance, I have run my code
with 1 dirichlet and 1 neumann, 2 neumann, and 1 absorbing condition and 1
neumann (as in step-24). When I run step-23 from the examples, which
I am using dealii version 9.1.1. The alternative you suggested has resolved
my issue.
Thank you,
Nikki
On Tuesday, March 9, 2021 at 4:16:36 AM UTC-5 Wolfgang Bangerth wrote:
> On 3/7/21 10:16 PM, Nikki Holtzer wrote:
> > Thank you for your input! That was my first attempt at tack
that
>> if statement above other than
>>
>> (face->boundary_id() == 1)
>>
>> or is this the correct approach and it needs to be amended somehow?
>> On Sunday, March 7, 2021 at 1:49:12 PM UTC-5 Wolfgang Bangerth wrote:
>>
>>> On 3/7/21 5:30 AM, N
2021 at 1:49:12 PM UTC-5 Wolfgang Bangerth wrote:
> On 3/7/21 5:30 AM, Nikki Holtzer wrote:
> >
> > For a 1d problem, I would like to implement a "boundary matrix" for an
> > absorbing boundary condition such as the one presented in Step 24 for
> only the
>
Hello everyone!
For a 1d problem, I would like to implement a "boundary matrix" for an
absorbing boundary condition such as the one presented in Step 24 for only
the left boundary. Is it necessary to have the 'if' statement be a
conditional over faces like in Step 24 even though it is a 1d prob
), this is a
> quadratic operation, not a linear one. You cannot represent this
> operation with a linear operator. Did you intend to simply scale with
> that value?
>
> Best,
> Matthias
>
>
> On Mon, Oct 26, 2020, at 13:40 CDT, Nikki Holtzer <
> nhol...@math.ari
Hello all,
I am struggling to turn a matrix scalar product into a linear operator. I
currently have:
temp = M.matrix_scalar_product(u,u);
and I need to turn 'temp' into a linear operator.
I have tried turning 'M' into a linear operator first and then performing
the multiplication and have
Yes it appears to be the system_matrix_inverse because I can print out the
residual and it is indeed not empty.
On Saturday, October 17, 2020 at 2:32:02 PM UTC-4 Wolfgang Bangerth wrote:
> On 10/16/20 10:35 AM, Nikki Holtzer wrote:
> >
> > The error I'm receiving i
ers the error, and why. There are a
> number
> of video lectures about running a program in a debugger.
>
> Best
> WB
>
>
> On 10/14/20 6:14 PM, Nikki Holtzer wrote:
> > The first error seems to be occurring at:
> >
> > affine_constraints.distribute_local_
,
Nikki
On Wednesday, October 14, 2020 at 10:14:34 AM UTC-4 Nikki Holtzer wrote:
> Ok thank you that has resolved most of the adding issue. I am still
> receiving the following error. It reads to me that there's still an error
> on the 'FullMatrix::add()' function. I have
ssemble()
#4 ./gravwave: OfflineData<1>::prepare()
#5 ./gravwave: main
On Tuesday, October 13, 2020 at 3:37:03 PM UTC-4 Wolfgang Bangerth wrote:
> On 10/13/20 10:28 AM, Nikki Holtzer wrote:
> > The issue that I have when all matrices are FullMatrix matrices is I
> can&#
convert the vector
without error either.
Thank you,
Nikki
On Monday, October 12, 2020 at 7:27:35 PM UTC-4 Wolfgang Bangerth wrote:
> On 10/11/20 6:09 PM, Nikki Holtzer wrote:
> > I was hoping to take the outer product which is formed in a FullMatrix,
> called
> > Kroneck
e to show us the code you use to
> be
> sure.
>
> I don't understand the connection between the code you show and the outer
> product. The outer product is a dense matrix, but the sparsity pattern you
> create is sparse. Were you hoping to put the outer product of two vectors
/20 2:51 PM, Nikki Holtzer wrote:
> >
> > When doing so I receive the following error:
> >
> >
> > An error occurred in line <195> of file
> > in function
> >
> > dealii::SparseMatrix&
> > dealii::SparseMatrix::oper
Hello everyone,
I was just curious if anyone had any suggestions for the above 'const'
qualifier problem I have stated above?
Thanks!
Nikki
On Thursday, October 8, 2020 at 9:16:24 AM UTC-4 Nikki Holtzer wrote:
> When I do the above, with the modifications you have provided, I st
::SparseMatrix&)
const’
I have not added the 'const' qualifier to my definitions of either of my
sparse matrices.
On Wednesday, October 7, 2020 at 4:58:45 PM UTC-4 Wolfgang Bangerth wrote:
> On 10/7/20 2:51 PM, Nikki Holtzer wrote:
> >
> > When doing so
Hello everyone,
I am trying to recast a current FullMatrix as a
SparseMatrix in order to add this matrix to another
SparseMatrix.
I have tried something like this:
SparseMatrix Mat_Sparse;
Mat_Sparse.copy_from(Mat_Full);
Mat_Sparse_2.add(1.,Mat_Sparse);
When doing so I receive the follow
Thank you!
On Tuesday, October 6, 2020 at 8:31:47 PM UTC-4 Wolfgang Bangerth wrote:
> On 10/6/20 5:50 PM, Nikki Holtzer wrote:
> >
> > I am trying to form a cross product/ outer product of two vectors of
> type
> > deallii:Vector. I have attempted to use some of t
Hello everyone,
I am trying to form a cross product/ outer product of two vectors of type
deallii:Vector. I have attempted to use some of the built in
functions for the outer product from the Tensor Class but have had no luck.
I can't seem to get anything other than
error: no matching function
21 matches
Mail list logo
