Skip to content
Snippets Groups Projects
Commit 31147260 authored by Lionel GUEZ's avatar Lionel GUEZ Committed by GUEZ Lionel
Browse files

Do not orient nor compute the area

Do not orient nor compute the area of `cont_list(n_cont)`. When the
area is too small, then we have one more copy than before, of
`cont_list(n_cont)` to `out_cont`, but, on the other hand, we have one
less call to `ccw_orient` than before. So this change should not
impact performance. Note that the orientation and area of
`cont_list(n_cont)` were not used after `get_1_outerm`.

Preparing to keep `cont_list` in projection coordinates.
parent b7ef4289
No related branches found
No related tags found
No related merge requests found
......@@ -158,14 +158,14 @@ contains
! Assertion: n_cont <= n_max_cont - 1
end if
cont_list(n_cont)%area = spher_polyline_area(cont_list(n_cont)%polyline)
call ccw_orient(cont_list(n_cont))
out_cont = cont_list(n_cont)
out_cont%area = spher_polyline_area(out_cont%polyline)
if (cont_list(n_cont)%area < min_area) then
if (abs(out_cont%area) < min_area) then
out_cont = null_ssh_contour()
n_cont = 0
else
out_cont = cont_list(n_cont)
call ccw_orient(out_cont)
end if
end if test_n_points
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment