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

New Neumann boundary conditions

parent 64eb2b6d
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 1 addition and 1 deletion
src/one-body-sample.cc
+ 1
− 1
View file @ 50a5b95e
...@@ -108,7 +108,7 @@ void assemble_neumann(GridView const &gridView, FEBasis const &feBasis, ...@@ -108,7 +108,7 @@ void assemble_neumann(GridView const &gridView, FEBasis const &feBasis,
BoundaryPatch<GridView> neumannBoundary(gridView, neumannNodes); BoundaryPatch<GridView> neumannBoundary(gridView, neumannNodes);
LocalVectorType SampleVector(0); LocalVectorType SampleVector(0);
// FIXME: random values (time-dependent) // FIXME: random values (time-dependent)
SampleVector[0] = time; SampleVector[0] = 10 * sqrt(time);
SampleVector[1] = 0; SampleVector[1] = 0;
ConstantFunction<LocalVectorType, LocalVectorType> fNeumann(SampleVector); ConstantFunction<LocalVectorType, LocalVectorType> fNeumann(SampleVector);
NeumannBoundaryAssembler<GridType, LocalVectorType> neumannBoundaryAssembler( NeumannBoundaryAssembler<GridType, LocalVectorType> neumannBoundaryAssembler(
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment