Skip to content

Wrong answer for derivatives #9

Closed
@hratcliffe

Description

@hratcliffe

Evaluating first derivatives gives erroneous results at the top (right hand end) of the domain.

Looks like in the conversion from f77 the ngo=2 or 8 case introduced a bug. Line 311 in splpak.F90 should read if (z>0.0_wp) then that is, greater-than to match the "IF z .LE. 0. GOTO end"

Really handy to have a single file spline package though, cheers for the conversion work.

Heather R

Metadata

Metadata

Assignees

No one assigned

    Labels

    bugSomething isn't working

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions