OK everybody! You can relax. :-) I managed to spot the loony. After mucking around with valgrind, and before trying gdb, I had one more look at my code and *finally* saw the stupid thing that I had been doing.

In the call to .Fortran() I had a line

    nphi=as.integer(nphi),

but "nphi" was nowhere defined (!!!) in the R code. The name "nphi" appeared as an argument in the Fortran subroutine in question, but was nowhere actually *used*!!!

It seems that passing a non-existent value as an argument to a Fortran subroutine can *sometimes* confuse it. Understandably.

I think that this "nphi" was a left-over from an earlier version of the code. I must have changed the code so that nphi was no longer needed, but then forgot to remove it from some places. Psigh! I hate myself sometimes.

Anyhow, thanks to all those who took the time and made the effort to try to help me.

cheers,

Rolf

--
Technical Editor ANZJS
Department of Statistics
University of Auckland
Phone: +64-9-373-7599 ext. 88276

______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Reply via email to