[Problem] New loading

src/sand-wedge-data/boundaryconditions.py

+import math
+
 class neumannCondition:
     def __call__(self, relativeTime):
         return 0
 
 class velocityDirichletCondition:
     def __call__(self, relativeTime):
-        if relativeTime <= 0.25:
-            factor = 4*relativeTime;
-        else:
-            factor = 1
-        return -5e-5 * factor
+        finalVelocity = -5e-5
+        if relativeTime <= 0.1:
+            return finalVelocity * ( 1.0 - math.cos(relativeTime * math.pi / 0.1) ) / 2.0
+        return finalVelocity
 
 Functions = {
     'neumannCondition' : neumannCondition(),

src/sand-wedge.cc

@@ -258,7 +258,7 @@ int main(int argc, char *argv[]) {
     {
       double v_initial_const;
       velocityDirichletFunction.evaluate(0.0, v_initial_const);
-      assert(v_initial_const == 0.0);
+      assert(std::abs(v_initial_const) < 1e-14);
     }
     Vector vr_initial(fineVertexCount);
     vr_initial = 0.0;