Skip to content
Snippets Groups Projects
Commit de94c6e4 authored by Elias Pipping's avatar Elias Pipping Committed by Elias Pipping
Browse files

Use a time-dependent neumann condition

parent 2f3047e4
Branches
No related tags found
No related merge requests found
......@@ -101,8 +101,9 @@ void assemble_neumann(GridView const &gridView, FEBasis const &feBasis,
int run) { // constant sample function on neumann boundary
BoundaryPatch<GridView> neumannBoundary(gridView, neumannNodes);
LocalVectorType SampleVector(0);
// FIXME: random values (time-independent)
SampleVector[0] = 1;
// FIXME: random values (time-dependent)
SampleVector[0] =
run / 10.0; // If we leave this fixed, we get into trouble for some reason
SampleVector[1] = 0;
ConstantFunction<LocalVectorType, LocalVectorType> fNeumann(SampleVector);
NeumannBoundaryAssembler<GridType, LocalVectorType> neumannBoundaryAssembler(
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment