s);
for( unsigned int i=0; iset_properties(quantities);
}
}
On Tuesday, June 16, 2020 at 10:25:01 PM UTC-4, Wolfgang Bangerth wrote:
>
> On 6/16/20 1:46 PM, Andrew Davis wrote:
> >
> > and I have gotten what I expect. I have also tried attaching the
> particle
>
since it indicates an error in deal.II, not in your own program.
Try to come up with the smallest possible program that still demonstrates
the error and contact the deal.II mailing lists with it to obtain help.
On Tuesday, June 16, 2020 at 3:46:36 PM UTC-4, Andrew Davis wrote:
>
> I
I have a question about how to set the quantities for dealii::Particles
that are stored in a dealii::ParticleHandler. I have successfully created a
dealii::ParticleHandler and randomly placed particles. I have also
successfully solved a PDE, whose solution is stored in the vector solution.
Then
Hello,
I have not managed to solve this issue. I have gotten as far as diagnosing
that the error occurs when the refinement needs to happen "across" the
period boundary. For some reason, this does not happen and then when you go
to solve the problem you end up with an inconsistent mesh. For the
obal(3);
Attached is the full code.
On Thursday, February 27, 2020 at 5:45:30 PM UTC-5, Andrew Davis wrote:
>
> HI all,
>
> I'm trying to implement a time-dependent convection equation using DG
> elements with periodic boundary conditions on and adaptive mesh.
>
> However, wh
stem? Are you using MeshWorker,
> MatrixFree or do you assemble matrices on your own using FEFaceValues?
> Is your code working in serial and for global refinement?
>
> Best,
> Daniel
>
>
> Am Fr., 28. Feb. 2020 um 08:29 Uhr schrieb Andrew Davis >:
>
>> I
I just say a slight problem and changed this line the declaration of
periodicityVector to:
std::vector::cell_iterator>
> periodicityVector;
But I'm getting the same behavior.
On Friday, February 28, 2020 at 8:08:19 AM UTC-5, Andrew Davis wrote:
>
> Yes, here is the code I
/doxygen/deal.II/step_45.html.
>
> Best,
> Daniel
>
> Am Do., 27. Feb. 2020 um 17:45 Uhr schrieb Andrew Davis >:
>
>> HI all,
>>
>> I'm trying to implement a time-dependent convection equation using DG
>> elements with periodic boundary conditio
HI all,
I'm trying to implement a time-dependent convection equation using DG
elements with periodic boundary conditions on and adaptive mesh.
However, when I try to adapt the mesh using periodic using this code:
// estimate the error in each cell
Vector estimatedErrorPerCell(triangulation
n Friday, January 31, 2020 at 4:27:26 AM UTC-5, Wolfgang Bangerth wrote:
>
> On 1/30/20 2:48 PM, Andrew Davis wrote:
> > I found a fix for this by editing the deal.ii source code. I am very new
> to
> > deal.ii so it is extremely likely that there is a way to accomplish
dy
On Thursday, January 30, 2020 at 3:39:28 PM UTC-5, Andrew Davis wrote:
>
> Fair point. I reinstalled in Debug mode and now I'm failing at that
> Assert. Here is the error message (using the same code I attached
> previously):
>
> An error occurred in line <2891> of
20 at 2:28:55 PM UTC-5, Andrew Davis wrote:
>>
>>
>> from the file fe_values.impl.2.inst.in. The Assert(false) seems strange
>> to me---why doesn't it crash? (I'm in Release mode so that could answer
>> that). However, should it be calling a different functio
20 at 11:31:38 AM UTC-5, Andrew Davis wrote:
>
> As another diagnostic---I added cout statements in the function `
> get_function_values` in fe_values.cc. It looks like the values of the
> input variable `fe_function` are correct but the `dof_values` are always
> set to zero. Spec
lways prints zero.
>
> On Thursday, January 30, 2020 at 10:49:18 AM UTC-5, Wolfgang Bangerth
> wrote:
>>
>> On 1/30/20 8:45 AM, Andrew Davis wrote:
>> > I thought the same thing---that 'old_solution' would be a zero
>> vector---but
>> > w
Sure, attached is a *.cpp file that should run. The cout statement on line
87 prints non-zero values for `old_solution` but the cout statement online
92 always prints zero.
On Thursday, January 30, 2020 at 10:49:18 AM UTC-5, Wolfgang Bangerth wrote:
>
> On 1/30/20 8:45 AM, Andrew Davis
vector
that I expect.
On Thursday, January 30, 2020 at 10:39:27 AM UTC-5, Wolfgang Bangerth wrote:
>
> On 1/29/20 1:14 PM, Andrew Davis wrote:
> > *
> > *
> > *For some reason the feValues.get_function_values(old_solution,
> > old_solution_values
I'm not sure if this is helpful but as a diagnostic, replacing
MeshWorker::loop with this loop that I wrote myself:
const QGauss quadrature_formula(fe.degree + 1);
const unsigned int n_q_points = quadrature_formula.size();
FEValues fe_values(fe, quadrature_formula, update_values |
update_Jx
I'm trying to implement a time-dependent solver that assembles the system
using the MeshWorker::loop tool. However, I cannot figure out how get the
values of the solution at the previous timestep at each quadrature point.
Currently, I store the old solution in a class called "AdvectionProblem"
18 matches
Mail list logo
