Commit 6a979e89 authored by Robert Schweppe's avatar Robert Schweppe
Browse files

- cleaned up code and fixed bug

parent 70eaf174
Pipeline #95863 failed with stages
in 20 minutes
......@@ -1557,16 +1557,15 @@ contains
point(1) = mod_shift(point(1))
end if
end if
lastiPolygon = currentPolygon
call inpoly(point, polygonNodes(1:nNodes,:), result)
if (result > 0_i4) then
log_trace(*) 'compute_weights_poly: subcell ', iSubcell, ' in target polygon ', currentPolygon
indices(iSubcell) = currentPolygon
! advance counter
self%subcells(currentPolygon) = self%subcells(currentPolygon) + 1_i8
lastiPolygon = currentPolygon
exit
end if
lastiPolygon = currentPolygon
end do
end do
maxSubcells = maxval(self%subcells)
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment