HDTenEightyP/NetNews huggingface dataset:
From: agnew@gems.vcu.edu
Newsgroups: comp.lang.fortran
Subject: Re: Public domain F77 compiler for DOS?
Date: 4 Jan 93 08:35:31 -0400
Organization: Medical College of Virginia

In article <bibhas.725658992@femto.engr.mun.ca>, bibhas@pico.engr.mun.ca (Bibhas Bhattacharya) writes:
> Hi!
> Is there any free FORTRAN compiler for DOS?
> Thanks.
> Bibhas.
> --
> 	+==========================+=========================+
> 	| login name: DR. Jackyl   |  In real life: Mr. Hyde.|
> 	+==========================+=========================+
> 		          bibhas@pico.engr.mun.ca  


Try looking up bcf77 on Simtel20  dunno if it's free, but it's shareware,
so if you are hardhearted, it's the same thing...

         /^^^\   \ /   Jim Agnew         | AGNEW@RUBY.VCU.EDU  (Internet)
        /      >  ||   Neurosurgery,     | AGNEW@VCUVAX        (Bitnet)
   /\_/     '   \  /   MCV-VCU           | This disc will self destruct in
 /________________>    Richmond, VA, USA | five seconds.  Good luck, Jim..."

Newsgroups: comp.lang.fortran
From: jlg@cochiti.lanl.gov (J. Giles)
Subject: Re: Is FORTRAN a viable language?
Organization: Los Alamos National Laboratory
Date: Mon, 4 Jan 1993 20:12:29 GMT

In article <1hqnn7INNfac@shelley.u.washington.edu>, chew@accel.aa.washington.edu (Gilbert Chew) writes:
|> [...]
|> However, while perusing the newsgroups for help
|> and information in buying said compiler, I noticed
|> that the C/C++ newsgroup seems to have about
|> 100 times the traffic of the FORTRAN newsgroup!

If you actually look at those newsgroups, you will find that most of
the traffic consists of trivial questions and long discussions about
issues which *should be* well defined and easy to learn/use/describe
in any really well designed language.  In other words, it's not clear
whether the volume of traffic in those groups is because of greater
popularity of the language or because of poorer design.

-- 
J. Giles
Newsgroups: comp.lang.fortran
From: Svante Signell
Subject: GNU Fortran compiler ? 
Reply-To: Svante.Signell@era_t.ericsson.se
Organization: Ericsson
Date: Tue, 12 Jan 1993 09:31:46 GMT

 What is the current status of the GNU Fortran (F77) compiler?
 Does anyone know if it is still under development?
Newsgroups: comp.lang.fortran
From: weiping@eng2.eng.monash.edu.au (Weiping Hu Materials)
Subject: Re: Block Data
Organization: Faculty of Engineering, Monash University 
Date: Tue, 12 Jan 1993 00:17:48 GMT

In article <C0pDpw.AxK@fs7.ece.cmu.edu> winstead@faraday.ece.cmu.edu (Charles Holden Winstead) writes:
>From: winstead@faraday.ece.cmu.edu (Charles Holden Winstead)
>Subject: Block Data
>Date: Mon, 11 Jan 1993 18:51:30 GMT
>This is probably an easy question, but I didn't find it in either of my
>FORTRAN 77 texts.  I have a program which has a section labelled BLOCK DATA
>where it defines a bunch of variables as common and fills them with data.
>My question is when is this done, before execution?  Is there any official 
>description of this?  Also this program contains a bunch of subroutines

I *think* the BLOCK DATA are assigned there values at compilation, i.e., 
before execution.

>that I'd like to place into separate files and link their object files at
>compilation, and I was wondering how to handle common blocks and this 
>block data section.
>
>Thanks.
>
>-Chuck
>

Just treat the block data section as a subroutine, or keep it with the main 
program.  The variable reference will be resolved at link time.

Hope this will be of some help.

>
>P.S.  I don't see how multiplying an inverse matrix by a column vector is 
>faster than back substitution. (Previous Thread).  You don't gain the 
>advantage of solving 1 at a time as with back substitution.

Weiping Hu
Newsgroups: comp.lang.fortran
From: winstead@faraday.ece.cmu.edu (Charles Holden Winstead)
Subject: Block Data
Organization: Electrical and Computer Engineering, Carnegie Mellon
Date: Mon, 11 Jan 1993 18:51:30 GMT

This is probably an easy question, but I didn't find it in either of my
FORTRAN 77 texts.  I have a program which has a section labelled BLOCK DATA
where it defines a bunch of variables as common and fills them with data.
My question is when is this done, before execution?  Is there any official 
description of this?  Also this program contains a bunch of subroutines
that I'd like to place into separate files and link their object files at
compilation, and I was wondering how to handle common blocks and this 
block data section.

Thanks.

-Chuck


P.S.  I don't see how multiplying an inverse matrix by a column vector is 
faster than back substitution. (Previous Thread).  You don't gain the 
advantage of solving 1 at a time as with back substitution.
From: cabhhmm@cli74aa.NoSubdomain.NoDomain (Michele Morin 5115 A408)
Newsgroups: comp.lang.fortran
Subject: Re: Name of Fortran Standard again (was Dynamic Memory Allocation)
Date: 5 Jan 93 14:12:26 GMT
Organization: EDF Direction des Etudes et Recherches

In the ISO/IEC 1539 : 1991 (E) document, I find, at the beginning, 
in the Foreword :
FORTRAN 77 is an archival standard in the United States.
Fortran 90 was developed jointly by the X3 tecnical subcommittee X3J3 
and IEC/ISO JTC1/SC22/WG5. It is identical to ISO/IEC 1539 : 1991, the
International Fortran Standard.

Michele Morin.
From: corbett@lupa.Eng.Sun.COM (Robert Corbett)
Newsgroups: comp.lang.fortran
Subject: Re: Name of Fortran Standard again (was Dynamic Memory Allocation)
Date: 9 Jan 1993 01:22:16 GMT
Organization: Sun
NNTP-Posting-Host: lupa

In article <1993Jan5.141226.25607@edf.fr> cabhhmm@cli74aa.NoSubdomain.NoDomain (Michele Morin 5115 A408) writes:
>In the ISO/IEC 1539 : 1991 (E) document, I find, at the beginning, 
>in the Foreword :
>FORTRAN 77 is an archival standard in the United States.
>Fortran 90 was developed jointly by the X3 tecnical subcommittee X3J3 
>and IEC/ISO JTC1/SC22/WG5. It is identical to ISO/IEC 1539 : 1991, the
>International Fortran Standard.
>
>Michele Morin.

It seemed odd to me that ISO/IEC 1539 : 1991 (E) would claim to be identical
to ISO/IEC 1539 : 1991, so I pulled out my copy of ISO/IEC 1539 : 1991 (E)
and checked.  I found no such claim in the foreword.  There also is no reference
to FORTRAN 77, X3J3, or IEC/ISO JTC1/SC22/WG5.  It does say

     International Standard ISO/IEC 1539 was prepared by
     Joint Technical Committee ISO/IEC 1, Information technology.

My guess is that the quotes in question came from a draft of the standard,
not the standard.

					Yours truly,
					Robert Corbett
From: cse0603@desire.wright.edu
Newsgroups: comp.lang.fortran
Subject: Compilier for SUN
Date: 11 Jan 93 12:42:19 EST
Organization:  Wright State University 

Hello

Does anyone know of a Fortran compilier for the SUN Sparc?  I need one that
will allow me to download the executable to a 68030 Single Board Computer.
Any pointers would be appreciated.

Thanks.

-Steve
cse0603@desire.wright.edu
Newsgroups: comp.lang.fortran
Organization: Penn State University
Date: Fri, 8 Jan 1993 12:54:09 EST
From: H. D. Knoble <HDK@psuvm.psu.edu>
Newsgroups: comp.lang.fortran
Subject: Re: inverse matrix

Regardless of the kind of matrix that is being dealt with, the code
(e.g., IMSL, NAG) should check for computational singularity (or in
the case of generalized inverse check that the generalized inverse is
also a least squares inverse and set an appropriate flag). The code
should also effect as accurate an inverse as feasible on each platform
it is run on. Else one can get an "inverse" that is just plain wrong
too easily.

If this matrix inversion is for an important project, it might be wise
to stay away from Numerical Recipes algorithms (and textbook code
in general with a few exceptions) as they are not usually robust in
the above sense. What this amounts to is that good algorithms are most often
more than taking the algebra and applying FORmula TRANslation to it.
So know WHO you're getting code from, and what algorithm(s) are being
used.
Newsgroups: comp.lang.fortran
From: jayas@ee.umanitoba.ca (Rohitha P. Jayasinghe)
Subject: Re: HELP: calling c from fortran 77
Organization: Dept. of Elect. & Comp. Eng., U of Manitoba
Date: Fri, 8 Jan 1993 18:07:31 GMT

In article <C0IBxD.JL3@ccu.umanitoba.ca>, jayas@ee.umanitoba.ca (Rohitha P. Jayasinghe) writes:
> Is it possible to call a C subroutine from fortran 77?
> If so, how can that be done?
> 
> thanks.
> 
> Rohitha.

Many thanks to those who replied.
I should have included more information in my original posting.

 - I am working on a Sun workstaion.
 - operating system Unix / Sun 4.1.2 
 - Fortran compiler Sun FORTRAN 1.4
 - C compiler Sun C 1.1

I have to call several C subroutines.
Argument types vary from integer/real
 to arrays/single variables.

Thanks again.

Rohitha.
Newsgroups: comp.lang.fortran
From: michael@chpc.utexas.edu (Michael Lemke)
Subject: Re: A VMS Fortran I/O question
Organization: The University of Texas System - CHPC
Date: Wed, 6 Jan 93 17:45:59 GMT

In article <MAINE.93Jan5095028@altair.dfrf.nasa.gov> maine@altair.dfrf.nasa.gov (Richard Maine) writes:
>On Tue, 5 Jan 93 06:19:34 GMT, michael@chpc.utexas.edu (Michael Lemke) said:
>
>         [... about carriage control interpretation...]
>
[... and the standard's definition of it ...]

Thank you for your answer, Richard, that was something I've long been
wondering about.  (Thanks also to Keith Bierman <Keith.Bierman@Eng.Sun.com>
who sent me private e-mail to the same question.)

Michael
-- 
Michael Lemke
Astronomy, UT Austin, Texas
(michael@io.as.utexas.edu or UTSPAN::UTADNX::IO::MICHAEL [SPAN])
Newsgroups: comp.lang.fortran
From: Seokmin Ahn <smahn@vx.cis.umn.edu>
Subject: Help: Fortran programming guideline anywhere?
X-Xxmessage-Id: <A773472CD301C378@me51-3-120.me.umn.edu>
X-Xxdate: Fri, 8 Jan 93 14:58:52 GMT
Organization: University of Minnesota
X-Useragent: Nuntius v1.1.1d17
Date: Fri, 8 Jan 1993 20:55:20 GMT

Hi,  netters,

I am looking for a programming guideline for Fortran.
Would you let me know where I can download it?
Or would you send me one if you have one?

Thank you.
Newsgroups: comp.lang.fortran
From: walt@netcom.com (Walt Brainerd)
Subject: Re: Name of Fortran Standard again (was Dynamic Memory Allocation)
Organization: Netcom - Online Communication Services  (408 241-9760 guest) 
Date: Mon, 11 Jan 1993 18:31:03 GMT

In article <lksa68INNkpt@exodus.Eng.Sun.COM>, corbett@lupa.Eng.Sun.COM (Robert Corbett) writes:
> In article <1993Jan5.141226.25607@edf.fr> cabhhmm@cli74aa.NoSubdomain.NoDomain (Michele Morin 5115 A408) writes:
> >In the ISO/IEC 1539 : 1991 (E) document, I find, at the beginning, 
> >in the Foreword :
> >FORTRAN 77 is an archival standard in the United States.
> >Fortran 90 was developed jointly by the X3 tecnical subcommittee X3J3 
> >and IEC/ISO JTC1/SC22/WG5. It is identical to ISO/IEC 1539 : 1991, the
> >International Fortran Standard.
> >
> >Michele Morin.
> 
> It seemed odd to me that ISO/IEC 1539 : 1991 (E) would claim to be identical
> to ISO/IEC 1539 : 1991, so I pulled out my copy of ISO/IEC 1539 : 1991 (E)
> and checked.  I found no such claim in the foreword.  There also is no reference
> to FORTRAN 77, X3J3, or IEC/ISO JTC1/SC22/WG5.  It does say
> 
>      International Standard ISO/IEC 1539 was prepared by
>      Joint Technical Committee ISO/IEC 1, Information technology.
> 

The foreword pages of drafts, the ISO standard, and possibly the
ANSI version of the standard (if there does turn out to be one)
are DIFFERENT.  The rest of the document is the same.  So I think
everybody is right (or wrong) on this one.
-- 
Walt Brainerd                walt@netcom.com
Unicomp, Inc.                +1-415-949-4052
235 Mt. Hamilton Ave.        +1-415-949-4058 (fax)
Los Altos, CA 94022 USA      "F90--Fortran for the '90s"
From: agnew@gems.vcu.edu
Newsgroups: comp.lang.fortran
Subject: Re: Disk-based Arrays
Date: 6 Jan 93 08:34:37 -0400
Organization: Medical College of Virginia

In article <1993Jan4.115514.1@camins.camosun.bc.ca>, 
alumni107@camins.camosun.bc.ca writes:
> I would like to find source for an implementation of 
> disk-based giant arrays-for those of us running 
> on systems without virtual memory :-( 
> 
> Can anyone suggest an ftp-able  set of routines that do this ?
> 
> Jeff Johnson; .sig agnostic
> sjohnson@post.RoyalRoads.ca

Jeff, could you try simple direct reads/writes?  I know PDP-11s have 
virtual arrays, is there something similar for you?

Jim 
         /^^^\   \ /   Jim Agnew         | AGNEW@RUBY.VCU.EDU  (Internet)
        /      >  ||   Neurosurgery,     | AGNEW@VCUVAX        (Bitnet)
   /\_/     '   \  /   MCV-VCU           | This disc will self destruct in
 /________________>    Richmond, VA, USA | five seconds.  Good luck, Jim..."

Newsgroups: comp.lang.fortran
From: mccalpin@perelandra.cms.udel.edu (John D. McCalpin)
Subject: Re: inverse matrix
Organization: College of Marine Studies, U. Del.
Date: Mon, 11 Jan 1993 21:25:37 GMT

In article <EIJKHOUT.93Jan11153758@cupid.cs.utk.edu> eijkhout@cupid.cs.utk.edu (Victor Eijkhout) writes:
>In article <C0p5JD.M0L@news.udel.edu> mccalpin@perelandra.cms.udel.edu (John D. McCalpin) writes:
>
>  In article <1993Jan8.201645.14915@news.eng.convex.com> dodson@convex.COM (Dave Dodson) writes:
>  >I'd like to point out that it is almost never required or desirable to
>  >compute the inverse of a matrix.
>
>  The direct use of the inverse matrix is generally the fastest way to
>  solve a dense system of equations with multiple, consecutive right-hand-sides
>  (as in a time-dependent fluid dynamics problem).
>
>Where do dense systems come from in fluid dynamics? Usually
>differential equations (partial or otherwise) give sparse
>matrices, and then calculating the inverse is at least a major
>waste of space.

The matrices are dense for spectral integration or differentiation
using any basis functions except trig functions.   The last time I
checked, this was still the fastest way to solve separable Poisson-like
equations using Chebyshev discretization, for example.

Dense matrices alse arise in the capacitance matrix method for solving
the elliptic PDE's that often arise in fluid dynamics problems.
-- 
--
John D. McCalpin                        mccalpin@perelandra.cms.udel.edu
Assistant Professor                     mccalpin@brahms.udel.edu
College of Marine Studies, U. Del.      John.McCalpin@mvs.udel.edu
Newsgroups: comp.lang.fortran
From: michael@chpc.utexas.edu (Michael Lemke)
Subject: Re: pretty printers for fortran
Organization: The University of Texas System - CHPC
Date: Sat, 9 Jan 93 21:50:23 GMT

In article <1993Jan9.202750.21806@sol.ctr.columbia.edu> anselmo@thermsa.eng.sunysb.edu (Andrew P. Anselmo) writes:
>Anybody know of pretty printers for fortran source code (for a Sun SS2 box?)
>I'm looking for stuff I can get off the net. 
>

Check out toolpack.  It will do much more like restructuring Spaghetti 
Fortran (tm).  Archie says toolpack is available from

Host cs.dal.ca

    Location: /pub/comp.archives
      DIRECTORY drwxrwxrwx        512  Mar  8 1991  toolpack

Host evolution.bchs.uh.edu

    Location: /pub
      DIRECTORY drwxr-xr-x       1024  Oct 19 1991  toolpack

Host ee.utah.edu

    Location: /
      DIRECTORY drwxr-xr-x       1024  Apr 13 1991  toolpack

Host walhalla.germany.eu.net

    Location: /pub/utils
      DIRECTORY drwxrwxr-x        512  May  7 1992  toolpack

Host unix.hensa.ac.uk

    Location: /pub/misc/unix
      DIRECTORY drwxr-xr-x        512  Nov  9 14:34  toolpack


Michael
-- 
Michael Lemke
Astronomy, UT Austin, Texas
(michael@io.as.utexas.edu or UTSPAN::UTADNX::IO::MICHAEL [SPAN])
Newsgroups: comp.lang.fortran
From: ejh@khonshu.colorado.edu (Edward J. Hartnett)
Subject: Re: Is FORTRAN a viable language?
Organization: CIRES, University of Colorado
	<1993Jan4.201229.12073@newshost.lanl.gov> <29250@oasys.dt.navy.mil>
Date: 5 Jan 93 09:20:24

In article <29250@oasys.dt.navy.mil> roth@oasys.dt.navy.mil (Pete Roth) writes:

   In comp.lang.fortran, jlg@cochiti.lanl.gov (J. Giles) writes:
   >u (Gilbert Chew) writes:
   >|> that the C/C++ newsgroup seems to have about
   >|> 100 times the traffic of the FORTRAN newsgroup!
   >
   >If you actually look at those newsgroups, you will find that most of
   >the traffic consists of trivial questions and long discussions about
   >issues which *should be* well defined and easy to learn/use/describe
   >in any really well designed language.  In other words, it's not clear
   >whether the volume of traffic in those groups is because of greater
   >popularity of the language or because of poorer design.

   A third possibility: the concepts of the language are new, and so the
   terms "easy to learn/use/describe" are not yet applicable. When we all
   understand what virtual methods, templates, overloading, private|protected|
   public, friend, class, etc. mean, *then* the issues will be well defined.

   Me, I'm just lurking on the C++ fora, tryina figure it out...

Actually, the real reason is that FORTRAN programmers are a lot
smarter than C programmers, who often don't even know how to use
modern plumbing or pick out their clothes in the morning. Studies show
that 72% of all C programmers have "severe social, if not sexual,
problems, and few problem solving skills, but an abundance of problem
causing skills." (from "Studies of Programmers: Laguage and Body
Odour, is there a connection?" British Journal of the Social
Sciences).

				Ed


--
	Don't blame me, I voted against Amendment 2!

Edward Hartnett				ejh@khonshu.colorado.edu
From: kvetzal@watarts.uwaterloo.ca (Ken Vetzal)
Subject: Reaching NAG via e-mail?
Originator: kvetzal@watarts.uwaterloo.ca
Organization: University of Waterloo
Date: Fri, 8 Jan 1993 21:39:34 GMT

Does anyone know if NAG has an e-mail address which accepts reports of 
suspected bugs?  I find that calling one of their random number 
generator routines works fine if I request an even number of random
numbers but blows up (produces Infinity - NaN) if an odd number of
random numbers is requested.  A sample program which does this is
included below:


	implicit double precision(a-h,o-z) 
	parameter(maxobs=365)
	double precision e(maxobs)
	external g05cbf, g05fdf
	call g05cbf(170565)
	do 100 n=1,2
	call g05fdf(0.d0,1.d0,maxobs,e)
	do 101 i=1,maxobs
	write(8,*)i,e(i)
101	continue
100	continue
	stop
	end

This will blow up with maxobs set to 365 or 367 or any other odd
number I've tried > 100, but seems to work fine at 364 or 366 or
any other even number > 100.  Thanks in advance for any advice.

Ken Vetzal (kvetzal@watarts.uwaterloo.ca)
From: eijkhout@cupid.cs.utk.edu (Victor Eijkhout)
Newsgroups: comp.lang.fortran
Subject: Re: inverse matrix
Date: 11 Jan 93 15:37:58
Organization: /pearl/homes/eijkhout/.organization
	<1993Jan8.201645.14915@news.eng.convex.com> <C0p5JD.M0L@news.udel.edu>
NNTP-Posting-Host: cupid.cs.utk.edu
In-reply-to: mccalpin@perelandra.cms.udel.edu's message of Mon, 11 Jan 1993 15:54:48 GMT

In article <C0p5JD.M0L@news.udel.edu> mccalpin@perelandra.cms.udel.edu (John D. McCalpin) writes:

   In article <1993Jan8.201645.14915@news.eng.convex.com> dodson@convex.COM (Dave Dodson) writes:
   >I'd like to point out that it is almost never required or desirable to
   >compute the inverse of a matrix.  Almost without exception, you can do
   >any computation in which you would use the inverse in a better way that
   >does not use the inverse.  By 'better' I mean faster, uses less memory,
   >more accurate, etc.

   The direct use of the inverse matrix is generally the fastest way to
   solve a dense system of equations with multiple, consecutive right-hand-sides
   (as in a time-dependent fluid dynamics problem).

Where do dense systems come from in fluid dynamics? Usually
differential equations (partial or otherwise) give sparse
matrices, and then calculating the inverse is at least a major
waste of space.

--
Victor Eijkhout ................................ `There are also a few bugs,
Department of Computer Science .......... though not as many as I've come to
University of Tennessee ............... expect in new [MS]Windows products.'
Knoxville TN 37919 ................................ (from a software review)
Newsgroups: comp.lang.fortran
From: mccalpin@perelandra.cms.udel.edu (John D. McCalpin)
Subject: Re: inverse matrix
Organization: College of Marine Studies, U. Del.
Date: Mon, 11 Jan 1993 15:54:48 GMT

In article <1993Jan8.201645.14915@news.eng.convex.com> dodson@convex.COM (Dave Dodson) writes:
>I'd like to point out that it is almost never required or desirable to
>compute the inverse of a matrix.  Almost without exception, you can do
>any computation in which you would use the inverse in a better way that
>does not use the inverse.  By 'better' I mean faster, uses less memory,
>more accurate, etc.

The direct use of the inverse matrix is generally the fastest way to
solve a dense system of equations with multiple, consecutive right-hand-sides
(as in a time-dependent fluid dynamics problem).

On the Cray Y series, for example, multiplying by the inverse matrix is
always (?) faster than performing the forward-backward substitution
step for each new RHS.  This is true even though both operations have
the same computational complexity in terms of operation counts (=2 N^2).
The difference is that matrix multiplication is a significantly
simpler algorithm with less overhead than forward-backward substitution.   

I suspect that the same sorts of speedups would be observed on Dodson's
Convex machines....
-- 
--
John D. McCalpin                        mccalpin@perelandra.cms.udel.edu
Assistant Professor                     mccalpin@brahms.udel.edu
College of Marine Studies, U. Del.      John.McCalpin@mvs.udel.edu
From: ross@turock.psych.upenn.edu (Ross Porter)
Newsgroups: comp.lang.fortran
Subject: Fortran for RS/6000
Date: 12 Jan 93 22:28:55 GMT
Organization: University of Pennsylvania

Anyone know of RS/6000 Fortran compilers other than IBM's?
Impressions appreciated.

Ross Porter
ross@psych.upenn.edu
From: nagrwb@vax.oxford.ac.uk
Newsgroups: comp.lang.fortran
Subject: Re: Reaching NAG via e-mail?
Date: 11 Jan 93 09:08:17 GMT
Organization: Oxford University VAX 6620

In article <C0K1Hy.BMJ@watserv2.uwaterloo.ca>, kvetzal@watarts.uwaterloo.ca (Ken Vetzal) writes:
> Does anyone know if NAG has an e-mail address which accepts reports of 
> suspected bugs?  
Certainly:    infodesk@nag.co.uk
              infodesk@nag.com   for North American customers
> I find that calling one of their random number 
> generator routines works fine if I request an even number of random
> numbers but blows up (produces Infinity - NaN) if an odd number of
> random numbers is requested.  A sample program which does this is
> included below:
> 
> 
> 	implicit double precision(a-h,o-z) 
> 	parameter(maxobs=365)
> 	double precision e(maxobs)
> 	external g05cbf, g05fdf
> 	call g05cbf(170565)
> 	do 100 n=1,2
> 	call g05fdf(0.d0,1.d0,maxobs,e)
> 	do 101 i=1,maxobs
> 	write(8,*)i,e(i)
> 101	continue
> 100	continue
> 	stop
> 	end
> 
> This will blow up with maxobs set to 365 or 367 or any other odd
> number I've tried > 100, but seems to work fine at 364 or 366 or
> any other even number > 100.  Thanks in advance for any advice.
> 
> Ken Vetzal (kvetzal@watarts.uwaterloo.ca)
This error in G05FDF is known to us. More precisely, the error only occurs
if two or more calls are made to G05FDF without an intervening call of G05CBF
or G05CCF, and if at least one of the values of N in those calls is odd.
The corrected version of the routine will be in the next release of
the NAG Fortran Library and in a future error bulletin.

Our apologies for the inconvenience that this has caused.

Richard Brankin
-- 
** R.W. Brankin -- Numerical Algorithms Group Ltd. -- nagrwb@vax.ox.ac.uk **
**   suggestions for sig. welcome   --   will need   include<std.discl>   **
Newsgroups: comp.lang.fortran