Add (ccr, lf) x (p3, p4) computation

diff --git a/src/psp.cc b/src/psp.cc
index 9b18fae9bb9b53a7b4a66fdaab89d2c820aabc05..99f5d3d71d1c1f0e460ce6fd287e291891cfee41 100644 (file)
--- a/src/psp.cc
+++ b/src/psp.cc
@@ -545,7 +545,17 @@ void PSPlanner::fe_parallel()
                                 std::get<1>(cli2), std::get<2>(cli2),
                                 std::get<3>(cli2), std::get<4>(cli2)
                         );
-                        // TODO (ccr, lf) x (p3, p4)
+                        auto cli3 = ::intersect(
+                                ccrx, ccry, ccr_lf,
+                                this->ps().x3(), this->ps().y3(),
+                                this->ps().x4(), this->ps().y4()
+                        );
+                        double a3 = smaller_angle_of_two(
+                                this->cc().lfx(), this->cc().lfy(),
+                                ccrx, ccry,
+                                std::get<1>(cli3), std::get<2>(cli3),
+                                std::get<3>(cli3), std::get<4>(cli3)
+                        );
 {
         double rf = sqrt(
                 pow(this->cc().lfy() - ccry, 2)