[Rd] DierckxSpline segfault

[Spencer Graves]

      Thanks for this bug report.  I've replicated the problem using 
your excellent example.  Unfortunately, this may not be an easy bug to 
fix. 

      Might it be feasible for you to use either the 'periodicSpline' 
function in the 'splines' package or a Fourier basis in the 'fda' package? 

      Best Wishes,
      Spencer Graves

Андрей Парамонов wrote:
> I've just encountered a segfault when using DierckxSpline::percur
> function. Below is the minimal example which triggers the error:
>
> ---
>
> library(DierckxSpline)
> x <- 1:10
> y <- rep(0, 10)
> pspline <- percur(x, y)
>
> ---
>
>  *** caught segfault ***
> address (nil), cause 'memory not mapped'
>
> Traceback:
>  1: .Fortran("percur", iopt = as.integer(iopt), m = as.integer(m),
> x = as.single(x), y = as.single(y), w = as.single(w), k = k,     s =
> as.single(s), nest = as.integer(nest), n = as.integer(n),     knots =
> Knots, coef = coef, fp = single(1), wrk = wrk, lwrk = lwrk,     iwrk =
> iwrk, ier = integer(1))
>  2: curfitSS(xyw, s = s, knots = knots, n = n, from = From, to = To,
>   k = k, periodic = periodic, ...)
>  3: curfit.default(x, periodic = TRUE, ...)
>  4: curfit(x, periodic = TRUE, ...)
>  5: percur(x, y)
>
> Possible actions:
> 1: abort (with core dump, if enabled)
> 2: normal R exit
> 3: exit R without saving workspace
> 4: exit R saving workspace
> Selection:
>
> ---
>
> Package:       DierckxSpline
> Type:          Package
> Title:         R companion to "Curve and Surface Fitting with Splines"
> Version:       1.1-2
> Date:          2007-7-31
> Author:        Sundar Dorai-Raj
> Maintainer:    Sundar Dorai-Raj <sundar.dorai-raj at pdf.com>
> Description:   This package provides a wrapper to the FITPACK routines
>                written by Paul Dierckx. The original Fortran is
>                available from http://www.netlib.org/dierckx
> License:       GPL (>= 2)
> Depends:       R (>= 2.4.0), stats, lattice, PolynomF
> Suggests:      fda, splines
> Packaged:      Sun Jan 11 13:08:12 2009; spencerg
> Built:         R 2.8.1; i486-pc-linux-gnu; 2009-01-17 10:21:01; unix
>
> ---
>
> I'm ready to provide any additional info,
> Andrey
>
> ______________________________________________
> R-devel at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-devel
>
>