Skip to content

Commit

Permalink
"BenchEllispe": fix diff clang/gcc
Browse files Browse the repository at this point in the history
  • Loading branch information
jmmuller committed Jan 23, 2024
1 parent f782470 commit a8ec3c5
Showing 1 changed file with 15 additions and 13 deletions.
28 changes: 15 additions & 13 deletions MMVII/src/Geom2D/cEllipse2D.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -380,19 +380,21 @@ void cEllipse::BenchEllispe()
}

for (int aKp=0 ; aKp<10 ; aKp++)
{
cPt2dr aP1 = anEl.PtOfTeta(RandInInterval(0.1,10),RandInInterval(-M_PI,M_PI));
cPt2dr aPP = anEl.ProjOnEllipse(aP1);
if (Bug)
{
StdOut() << anEl.ToRhoTeta(aP1) << " " << Scal(anEl.Tgt(aPP),aP1-aPP) << aP1-aPP << std::endl;
StdOut() << "Ratiooo=" << aLSa/aLGa << std::endl;
}

MMVII_INTERNAL_ASSERT_bench(std::abs(anEl.ToRhoTeta(aPP).x()-1.0)<1e-5,"BenchEllispe");
MMVII_INTERNAL_ASSERT_bench(std::abs(Scal(anEl.Tgt(aPP),aP1-aPP))<1e-3,"BenchEllispe");
// getchar();
}
{
auto v1 = RandInInterval(0.1,10);
auto v2 = RandInInterval(-M_PI,M_PI) ;
cPt2dr aP1 = anEl.PtOfTeta(v1,v2);
cPt2dr aPP = anEl.ProjOnEllipse(aP1);
if (Bug)
{
StdOut() << anEl.ToRhoTeta(aP1) << " " << Scal(anEl.Tgt(aPP),aP1-aPP) << aP1-aPP << std::endl;
StdOut() << "Ratiooo=" << aLSa/aLGa << std::endl;
}

MMVII_INTERNAL_ASSERT_bench(std::abs(anEl.ToRhoTeta(aPP).x()-1.0)<1e-5,"BenchEllispe");
MMVII_INTERNAL_ASSERT_bench(std::abs(Scal(anEl.Tgt(aPP),aP1-aPP))<1e-3,"BenchEllispe");
// getchar();
}

/* Check experimentaly that ApproxDist ~ EuclidDist and that the approximation can be bounded as a fonction
* of ellipse excentricity*/
Expand Down

0 comments on commit a8ec3c5

Please sign in to comment.