Commit de86ca7d authored by Peter Bastian's avatar Peter Bastian
Browse files

Testversion OneDGrid

[[Imported from SVN: r58]]
parent 36e3750b
......@@ -129,7 +129,7 @@ int main (int argc , char ** argv)
// uc2.grid().globalRefine(3);
// timeloop(uc2.grid(),0.5,3,7);
uc0.grid().globalRefine(4);
timeloop(uc0.grid(),0.5,4,9);
timeloop(uc0.grid(),0.5,4,8);
#if HAVE_MPI
MPI_Finalize();
......
......@@ -14,7 +14,7 @@
template<int dimworld, class ct>
double c0 (const Dune::FieldVector<ct,dimworld>& x)
{
if (x.two_norm()<0.125)
if (x.two_norm()>0.125 && x.two_norm()<0.5)
return 1.0;
else
return 0.0;
......@@ -24,6 +24,7 @@ double c0 (const Dune::FieldVector<ct,dimworld>& x)
template<int dimworld, class ct>
double b (const Dune::FieldVector<ct,dimworld>& x, double t)
{
return 0.0;
if (x.two_norm()<t+0.125)
return 1.0;
else
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment