Skip to content

fix issue #58; set initial n to min. points#108

Open
manuelbb-upb wants to merge 2 commits intoJuliaMath:masterfrom
manuelbb-upb:perdioc_parametric
Open

fix issue #58; set initial n to min. points#108
manuelbb-upb wants to merge 2 commits intoJuliaMath:masterfrom
manuelbb-upb:perdioc_parametric

Conversation

@manuelbb-upb
Copy link

This should fix #58

There is a bug in https://www.netlib.org/dierckx/fpclos.f
Inspecting the goto-statements, we see that n is never properly initialized in case s>0.
In this case, it is assigned only after if(n.ge.nest) go to 620 (triggering an error code 1).

On our side, we can initialize the reference to n before calling the Fortran subroutine, this appears to fix the issue.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

Periodic ParametricSpline() does not accept an 's' argument

1 participant