Resolve "Reformulate penalty factor" and "Fiber intersection algorithm fails if fiber ends at intersection inside grid"
What does this MR do?
-
Reformulate the penalty factor to not include fiber segment information
-
Fix the fiber intersection algorithm in conjunction with updates to the local operator. Information from adjacent fiber segments with element intersections are not needed anymore. The algorithm now explicitly allows fibers ending exactly at a grid intersection.
-
Make intersection algorithm allow for fibers that begin at a grid intersection
-
Add test
test-fiber-intersection.cc
for verifying the fiber intersection algorithm
Related issues
Closes #10 (closed), #15 (closed)
Edited by Lukas Riedel