Skip to content

Resolve "Reformulate penalty factor" and "Fiber intersection algorithm fails if fiber ends at intersection inside grid"

Lukas Riedel requested to merge 15-reformulate-penalty-factor into master

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

Merge request reports