Add (ccr, rf) x (p2, p3) computation

diff --git a/src/psp.cc b/src/psp.cc
index b9b26a5802696662e56f0e3541a8ea17d3e57b2e..9b18fae9bb9b53a7b4a66fdaab89d2c820aabc05 100644 (file)
--- a/src/psp.cc
+++ b/src/psp.cc
@@ -534,7 +534,17 @@ void PSPlanner::fe_parallel()
                                 std::get<1>(cli1), std::get<2>(cli1),
                                 std::get<3>(cli1), std::get<4>(cli1)
                         );
-                        // TODO (ccr, rf) x (p2, p3)
+                        auto cli2 = ::intersect(
+                                ccrx, ccry, ccr_rf,
+                                this->ps().x2(), this->ps().y2(),
+                                this->ps().x3(), this->ps().y3()
+                        );
+                        double a2 = smaller_angle_of_two(
+                                this->cc().rfx(), this->cc().rfy(),
+                                ccrx, ccry,
+                                std::get<1>(cli2), std::get<2>(cli2),
+                                std::get<3>(cli2), std::get<4>(cli2)
+                        );
                         // TODO (ccr, lf) x (p3, p4)
 {
         double rf = sqrt(