HDTenEightyP/NetNews huggingface dataset:
From: michael@chpc.utexas.edu (Michael Lemke)
Subject: Re: A VMS Fortran I/O question
Organization: The University of Texas System - CHPC
Date: Tue, 5 Jan 93 06:19:34 GMT

In article <C0CGGB.1tE@news.cso.uiuc.edu> ercolessi@uimrl3.mrl.uiuc.edu (furio ercolessi) writes:
>In article <1993Jan3.040453.20688@chpc.utexas.edu>, michael@chpc.utexas.edu (Michael Lemke) writes:
>|> I'd rather think the *Unix* stuff should honor
>|>FORTRAN carriage control.  Well, I know it doesn't but a PRINT
>|>statement after all creates a print file (or even sends something to a
>|>printer), so doesn't FORTRAN require print files to honor the carriage
>|>control character?  What is the excuse this time that the Unix joke is
>|>still standard conform?
>
>They have put there a filter called asa, and you just pipe the output
>through it.  But then you will have a bit of a problem for interactive
>programs (buffering will be your enemy!), and also for "batch"
>programs writing on a file if you like to peek at the output while the
>program is running.

Yes, I know that beggar and associated problems.  In my opinion this is
not a solution.  With the same argumentation they could output troff or
TeX or postscript.  Doesn't matter if the filter is called asa or
nroff.  In either case you get a file which does not print propoerly
without further action.  I am still waiting for an explanation if this
Unix behavior is standard conform or not (sorry, I don't have a copy of
the standard to check myself).

>
>|>  BTW, in Sun Fortran there is a special STATUS='PRINT' in
>|>the OPEN statement which causes the first character to be ignored if a
>|>blank or to be converted to a CTRL/L if a 1 etc.  `+' does not work
>|>though. This evne works for stdout but I don't know about PRINT.  Maybe 
>|>your Unix compiler has something similar.
>
>This sounds a neat idea and I wholeheartly suggest that all Unix compiler
>writers put this feature in (hopefully, with the same syntax ...).  
>I only have two questions: how can you OPEN the unit where the PRINTs
>go (usually mapped to the Unix standard output) ?  OPEN(*,) ???
>And, are you sure that the keyword is 'STATUS'? I think that STATUS
>has another meaning associated with it and cannot be overloaded.

Yes, they do overload STATUS and this is why I don't consider it a nice
idea.  Also, I think all files should have this attribute set by
default (see above) and second, I wonder if the standard says anything
about the external (file) representation of the carriage control
character -- does it have to be in the file literally?  If not,
something like the Sun solution would really be an acceptable way to
deal with it under Unix. asa et al. certainly is not.

I don't know right now how you can affect `channel *' or PRINTs.  PRINTs
might have a negative channel number or something else special.  * might
not work at all.

Michael
-- 
Michael Lemke
Astronomy, UT Austin, Texas
(michael@io.as.utexas.edu or UTSPAN::UTADNX::IO::MICHAEL [SPAN])
From: fred@poly2.nist.gov (Frederick R. Phelan Jr.)
Newsgroups: comp.lang.fortran
Subject: Re: Structure and Data Statement Usage
Date: 12 Jan 93 22:49:41 GMT

In article <1993Jan8.234201.13599@netcom.com>, walt@netcom.com (Walt Brainerd) writes:
|> In article <7876@dove.nist.gov>, fred@poly2.nist.gov (Frederick R. Phelan Jr.) writes:
|> > Heres a small code prototype of my problem:
|> > 
|> >       program test_35
|> > 
|> >       structure /geom/
|> >         integer ngprops
|> >       end structure
|> >       record /geom/ g1,g2
|> > 
|> > c     uncomment these for compiler havoc
|> > c      data g1.ngprops /3/
|> > c      data g2.ngprops /5/
|> > 

|> It is unfortunate that a person at  N  I  !!!S!!!  T  is not using
|> a STANDARD language 8^).

Ahem, just experimenting :-)
Fred
From: fred@poly2.nist.gov (Frederick R. Phelan Jr.)
Newsgroups: comp.lang.fortran
Subject: Structure and Data Statement Usage
Date: 8 Jan 93 18:04:15 GMT

Heres a small code prototype of my problem:

      program test_35

      structure /geom/
        integer ngprops
      end structure
      record /geom/ g1,g2

c     uncomment these for compiler havoc
c      data g1.ngprops /3/
c      data g2.ngprops /5/

c     these are ok when uncommented
c     g1.ngprops=3
c     g2.ngprops=5

      write(*,1000) g1.ngprops
 1000 format('g1=',i2)
      write(*,1001) g2.ngprops
 1001 format('g2=',i2)

      stop
      end

Ok, so I have a structure, and I want to initialize
g1.ngprops and g2.ngprops to differnent values.

Compiling with the data statements uncommented gives:

        f77 test_35.f -o test_35.exe
Error on line 6 of test_35.f: syntax error
Error on line 6 of test_35.f: syntax error
*** Error code 1

Stop.


Note: The record statement is line 6.  

The program works ok with the lines 

     g1.ngprops=3
     g2.ngprops=5

uncommented.

I know you can say:

      structure /geom/
        integer ngprops /5/
      end structure

but then g1.ngprops and g2.ngprops are initialized to the same 
value.

How can you use structure variables in data statements?

Fred
---
Fred Phelan
fred@poly2.nist.gov


PS: This is on a Silicon Graphics machine.

PPS: Anyone care to comment on how common f77 extensions
     like structures have become on other platforms? I'll
     summarize if there is sufficient replies and/or
     interest.


Newsgroups: comp.lang.fortran
From: ercolessi@uimrl3.mrl.uiuc.edu (furio ercolessi)
Subject: Re: Fortran on Windows NT Questionnaire
Reply-To: ercolessi@uimrl3.mrl.uiuc.edu (furio ercolessi)
Organization: MRL - UIUC
Date: Thu, 7 Jan 1993 01:45:50 GMT

In article <1993Jan6.215523.11278@nntpd.lkg.dec.com>, wild@tle.enet.dec.com (Joe Wild) writes:
|>   Subject: Questionnaire about Fortran
|>
|>	The Fortran team at Digital Equipment Corporation is 
|>   developing Fortran for the Windows NT operating system. You can help us 
|>   try to meet your Fortran needs by completing this 66-question survey.
|>
|>	Each person who responds electronically on or before January 12 will 
|>   receive an electronic copy of "Zen and the Art of the Internet" by Brendan 
|>   P. Kehoe as a token of our appreciation.

Come on! Think something better!!! We all have it already, and anyway, getting it
is much less work than filling the questionnaire!!!
Hints:
* 1% additional discount on all purchases for the next 3 years
* a PDP11/70 to put in the kitchen to fry eggs
* a VMS shell for IBM RS/6k 
* one night on the net alone with a 200Mhz alpha
* a paper copy of the "Zen" guide, but without the guide actually printed
  on it                        :-)
--
Furio Ercolessi
Materials Research Laboratory           |   Intl School for Advanced Studies
Univ. of Illinois at Urbana-Champaign   |   Trieste, Italy
furio@uiuc.edu                          |   furio@sissa.it
Newsgroups: comp.lang.fortran
From: art@cyclone.MITRE.org (Arthur T. McClinton)
Subject: Re: NT questionnaire, E-Mail Ad. ?????
Originator: art@cyclone
Organization: The MITRE Corp. McLean Va.
Date: Thu, 7 Jan 1993 20:34:01 GMT


In article <1993Jan7.183643.5224@news.clarkson.edu>, russell@sun.soe.clarkson.edu (David Russell) writes:
> From article by Dave DeWolfe:
> > Hi.  I filled out the questionnaire, but the address to which to mail
> > it doesn't exist as far as my mailer is concerned.
> > 
> > Pls respond to me via email. Thanks.
> 
>    Same problem here. What's the scoop.
Appears that one of the gateways into DEC.COM does not know the alias
tko.  I was able to send to the original poster who gave me the following
address.
raven@kobal.enet.dec.com
t worked when I sent to it.  Apparantly many people are getting through to
the original address.

Newsgroups: comp.lang.fortran
From: walt@netcom.com (Walt Brainerd)
Subject: Re: Structure and Data Statement Usage
Organization: Netcom - Online Communication Services  (408 241-9760 guest) 
Date: Fri, 8 Jan 1993 23:42:01 GMT

In article <7876@dove.nist.gov>, fred@poly2.nist.gov (Frederick R. Phelan Jr.) writes:
> Heres a small code prototype of my problem:
> 
>       program test_35
> 
>       structure /geom/
>         integer ngprops
>       end structure
>       record /geom/ g1,g2
> 
> c     uncomment these for compiler havoc
> c      data g1.ngprops /3/
> c      data g2.ngprops /5/
> 
It is unfortunate that a person at  N  I  !!!S!!!  T  is not using
a STANDARD language 8^).

In standard Fortran 90, it goes like this:

program test_35

type geom
  integer ngprops
end type geom

type (geom) g1, g2

data g1%ngprops /3/
data g2%ngprops /5/

print*, g1, g2
end

I believe this will run on SGI machines using the NAG compiler,
an execellent tool for code development and learning until SGI
gets their own compiler.
-- 
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"
Newsgroups: comp.lang.fortran
From: walt@netcom.com (Walt Brainerd)
Subject: Flowcharts
Keywords: Fortran flowchart
Organization: Netcom - Online Communication Services  (408 241-9760 guest) 
Date: Wed, 6 Jan 1993 23:58:03 GMT

There were some inquiries recently here about
Fortran flowcharting programs.  I just spent some
time talking with someone who claims to have
such a product.  It is written in Microsoft
Fortran and generates PostScript output.  I
don't know anything else about it yet, but if
you want to inquire:

   Sylvia Hart
   Scientific Programming Consultants
   5246 Valley Ridge Avenue
   Los Angeles, CA 90043 USA
   +1-213-293-7653
   [No e-mail and no fax!!!]

Let us know if you find out any more useful
info about this.
-- 
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: eijkhout@cupid.cs.utk.edu (Victor Eijkhout)
Newsgroups: comp.lang.fortran
Subject: Re: inverse matrix
Date: 12 Jan 93 16:16:20
Organization: /pearl/homes/eijkhout/.organization
	<EIJKHOUT.93Jan11153758@cupid.cs.utk.edu> <C0pKup.vJ@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 21:25:37 GMT

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

   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.

But only on domains that are a Cartesian product of intervals, right?
Or can you handle separable equations on arbitrarily shaped domains
(Chessapeake bay?) this way?

   Dense matrices alse arise in the capacitance matrix method for solving
   the elliptic PDE's that often arise in fluid dynamics problems.

Isn't that after eliminating most of a *sparse* matrix onto
some subset of the variables (corresponding to an interior line
in the domain) so that there is an outer problem that is sparse?

   John D. McCalpin
--
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)
From: cstewart+@pitt.edu (Carolyn A Stewart)
Newsgroups: comp.lang.fortran
Subject: Mac Fortran?
Date: 4 Jan 93 21:57:27 GMT

I am looking for a version of Fortran that will run on a Macintosh
LC.  Right now i have 2 megs of memory but am looking into upgrading that...
Does anyone know of a Program that will allow me to do Fortran on the mac?

					Carolyn

Please send mail to:
IN%"cstewart@vm2.cis.pitt.edu"

Subject: Re: Real Programmers
Summary: 
Expires: 
Followup-To: 
Distribution: 
Organization: University of Pittsburgh
Keywords:
From: burley@apple-gunkies.gnu.ai.mit.edu (Craig Burley)
Newsgroups: comp.lang.fortran
Subject: Re: GNU Fortran compiler ?
Date: 12 Jan 93 06:13:22
Organization: Free Software Foundation 545 Tech Square Cambridge, MA 02139
NNTP-Posting-Host: apple-gunkies.gnu.ai.mit.edu
In-reply-to: Svante Signell's message of 12 Jan 93 09:31:46 GMT

In article <1993Jan12.093146.14165@ericsson.se> Svante Signell writes:

    What is the current status of the GNU Fortran (F77) compiler?
    Does anyone know if it is still under development?

Okay, here is the definitive info (at least for now :-) on how to get
current status on info on g77 (nothing has changed, but I want to
put the info in one place):

If you can do finger -l, try:

    finger -l fortran@gate.gnu.ai.mit.edu

If that produces output of around 100 lines, then it probably works.

Otherwise, if you can do ftp, try:

    ftp gate.gnu.ai.mit.edu
    Username: ftp
    Password: (any)
    get g77.plan

The ftp copy of the file is updated every 3 hours from the live copy.

Otherwise, you could try sending email to user fortran@gnu.ai.mit.edu
asking to be sent the latest copy of the file, but _please_ make this
your last option, because it requires human intervention!
--

James Craig Burley, Software Craftsperson    burley@gnu.ai.mit.edu
Member of the League for Programming Freedom (LPF) lpf@uunet.uu.net
Newsgroups: comp.lang.fortran
From: alumni107@camins.camosun.bc.ca
Subject: Disk-based Arrays
Summary: where can I find source for this ?
Reply-To: sjohnson@post.RoyalRoads.ca
Organization: Camosun College
Date: Mon, 4 Jan 1993 19:55:14 GMT

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
Newsgroups: comp.lang.fortran
From: maine@altair.dfrf.nasa.gov (Richard Maine)
Subject: Re: inverse matrix
Organization: NASA Dryden, Edwards, Cal.
	<1993Jan8.134612.18468@news.unige.ch>
Date: Fri, 8 Jan 1993 17:01:10 GMT

On Fri, 8 Jan 1993 13:46:12 GMT, pfennige@scsun.unige.ch (PFENNIGER Daniel) said:


PFENNIGER> Is it true that a Fortran-90 programme inverting matrix A would be:
PFENNIGER>       B = 1./A

No.  That would be a legal statement, but it makes each element of B
the inverse of the corresponding element of A.  The f90 arrays are
better thought of as "arrays" than as "matrices".  In particular,
all the operations are element-by-element.  There is a standard
matmul subroutine for the usual matrix multiplication, but there is
no matrix inversion.

I wasn't on the committee at the time, but I believe the rationale was
that there was no single matrix inversion technique good for all
situations, and that numerically-oriented users were better off
choosing a particular, appropriate one than whatever the vendor
happened to implement (which might be different on different
platforms).

You (or a library writer) can make a "matrix" (or whatever you want to
call it) user-defined type and define division for it such that the
above sample line would do the matrix inversion, but such a type
isn't part of the base language.

--
--
Richard Maine
maine@altair.dfrf.nasa.gov

Newsgroups: comp.lang.fortran
From: kbanerje@grad12.math.ncsu.edu (Kaushik Banerjee)
Subject: Re: simultaneous equations...solutions
Organization: NCSU Mathematics Dept.
Date: Sun, 10 Jan 1993 21:17:32 GMT

Fact #1

An nxn linear system Ax=b has a unique solution if
and only if A is nonsingular.

Fact #2

An mxn (where m .ne. n) linear system Ax=b has either 
no solution (termed "inconsistent") or an infinite number 
of solutions.

Fact #3

The best you can do is a "least squares" solution.
Notation: A' = transpose of A
          A^ = inverse of A 
          ||x|| = Euclidean norm of the vector x
Say you want to solve Ax=b where A is nx3
==> A'Ax = A'b  (these are the "normal equations")
    If the columns of A are linearly independent, then the
    3x3 matrix A'A is nonsingular.
==> x = (A'A)^A'b, i.e., x is the unique solution to the 
    normal equations and the unique least squares solution
    to Ax=b.
    
Definition of a least squares solution:
 x is least squares solution if and only if
 ||Ax - b|| .le. ||z - b|| for all other vectors z.

You may want to refer to "Introduction to Matrix 
Computations" by Gilbert Stewart or some other text
on numerical linear algebra.

Hope this helps.

K. Banerjee
From: dtm@castle.ed.ac.uk (D Muxworthy)
Newsgroups: comp.lang.fortran
Subject: Name of Fortran Standard again (was Dynamic Memory Allocation)
Date: 5 Jan 93 12:57:56 GMT
Organization: Edinburgh University

In article <1993Jan4.195942.9993@newshost.lanl.gov> jlg@cochiti.lanl.gov 
(J. Giles) writes (amongst other things):
>FORTRAN 77 is the designation recommended by the standard document itself
>(ANSI X3.9-1978, Forward):

>Similarly, "Fortran 90" is actually the *title* of the ISO version
>of the new standard.  Whether these designations are *official*
>or merely *officially informal*, they are clearly the appropriate
>names to use.

This is nit-picking, but let's get things right.  The ISO document ISO/IEC
1539 : 1991 (E) has as its title "Information technology - Programming
languages - Fortran", i.e.  no "90" and capitalized as shown.  The
introduction starts 'This International Standard specifies the form and
establishes the interpretation of programs expressed in the Fortran language
(known informally as "Fortran 90").' Thereafter the document generally refers
to the new language simply as "Fortran" and to the previous standard as
"FORTRAN 77", actually using a font with "ORTRAN" in small capitals. 

There is only one international standard, the 1991 one which superseded
ISO/IEC 1539:1980 (that is FORTRAN 77).  The US on the other hand has both
X3.9-1978 FORTRAN and X3.198-1992 Fortran 90.  The "90" is part of the
official name, needed to distinguish the two standards. 
David Muxworthy
From: cabhhmm@cli74aa.NoSubdomain.NoDomain (Michele Morin 5115 A408)
Newsgroups: comp.lang.fortran
Subject: Re: fortran flowcharts?
Date: 4 Jan 93 13:32:34 GMT
Organization: EDF Direction des Etudes et Recherches

In EDF (Electricite de France), we use 2 tools :
1) the CRAY' tool "ftref"...  but it's work only with a CRAY computer !
2) QA-Fortran from Programming Research Limited : this tool work on many
workstations and computers.I give you the company address :
   Waynflete House
   74-76 High Street
   ESHER
   SURREY KT10 9QS
   England             Phone 0372 4621 30  Fax 0372 4683 53

Happy New Year

Michele Morin   
MorinM@cli74ak.der.edf.fr
Newsgroups: comp.lang.fortran
From: ercolessi@uimrl3.mrl.uiuc.edu (furio ercolessi)
Subject: Re: A VMS Fortran I/O question
Reply-To: ercolessi@uimrl3.mrl.uiuc.edu (furio ercolessi)
Organization: MRL - UIUC
Date: Mon, 4 Jan 1993 19:21:45 GMT

In article <1993Jan3.040453.20688@chpc.utexas.edu>, michael@chpc.utexas.edu (Michael Lemke) writes:
|> I'd rather think the *Unix* stuff should honor
|>FORTRAN carriage control.  Well, I know it doesn't but a PRINT
|>statement after all creates a print file (or even sends something to a
|>printer), so doesn't FORTRAN require print files to honor the carriage
|>control character?  What is the excuse this time that the Unix joke is
|>still standard conform?

They have put there a filter called asa, and you just pipe the output
through it.  But then you will have a bit of a problem for interactive
programs (buffering will be your enemy!), and also for "batch"
programs writing on a file if you like to peek at the output while the
program is running.

|>  BTW, in Sun Fortran there is a special STATUS='PRINT' in
|>the OPEN statement which causes the first character to be ignored if a
|>blank or to be converted to a CTRL/L if a 1 etc.  `+' does not work
|>though. This evne works for stdout but I don't know about PRINT.  Maybe 
|>your Unix compiler has something similar.

This sounds a neat idea and I wholeheartly suggest that all Unix compiler
writers put this feature in (hopefully, with the same syntax ...).  
I only have two questions: how can you OPEN the unit where the PRINTs
go (usually mapped to the Unix standard output) ?  OPEN(*,) ???
And, are you sure that the keyword is 'STATUS'? I think that STATUS
has another meaning associated with it and cannot be overloaded.

Thanks to all who replied to my original post.

furio
--
Furio Ercolessi
Materials Research Laboratory           |   Intl School for Advanced Studies
Univ. of Illinois at Urbana-Champaign   |   Trieste, Italy
furio@uiuc.edu                          |   furio@sissa.it
From: n8743196@henson.cc.wwu.edu (Jeff Wandling)
Newsgroups: comp.lang.fortran
Subject: couple of simple questions
Summary: left justification and case statements
Keywords: left justification case
Date: 11 Jan 93 00:24:03 GMT
Article-I.D.: henson.1993Jan11.002403.16544
Organization: Howling Dog

I come from a C programming background.. System is a SS2, Sun's f77.

In a format statement, how do I specify left justify when writing
floating point numbers?

Is there a f77 construct for "switch" statements?

Reply via email, please. Thanks.
-- 
perl -e '$a="Japhuneasorcttlk h eae r r ,";$]=0;$s=$];while($s<4){do{$_ .= pack
("c",vec($a,$],8));$]+=4;}while($]<length($a)&&vec($a,$],8)!=32);$s++;$]=$s;$_
.= ' ';}chop;print;' Jeff Wandling  <jdw@huxley.wwu.edu> "will hack 4 food"
Newsgroups: comp.lang.fortran
From: dnl@convex.is.macsch.com (David Lombard)
Subject: Re: inverse matrix
Organization: MacNeal-Schwendler Corp.
Date: Sat, 9 Jan 93 00:21:26 GMT

In article <1993Jan8.134612.18468@news.unige.ch> pfennige@scsun.unige.ch (PFENNIGER Daniel) writes:
>
>Is it true that a Fortran-90 programme inverting matrix A would be:
>      B = 1./A

I _think_ that will produce a term-by-term inversion of array A into B. 
i.e. B(i,j) = 1./A(i,j)


Regards,
DNL

                 MY_COMMENTS = MY_OPINIONS = NOBODY_ELSES;

David N. Lombard    The first thing we do,         The MacNeal-Schwendler Corp
dnl@macsch.com      Let's kill all the lawyers.    815 Colorado Blvd
(213) 259-4911      II Henry VI, IV.ii             Los Angeles, CA  90041
Newsgroups: comp.lang.fortran
From: cabhhmm@cli74aa.NoSubdomain.NoDomain (Michele Morin 5115 A408)
Subject: Re: inverse matrix
Organization: EDF Direction des Etudes et Recherches
Date: Fri, 8 Jan 1993 10:03:15 GMT

In article <C0I49C.Jrr@athena.cs.uga.edu>, DEANRE@bscf.uga.edu (ROBIN DEAN) writes:
|>       Does anyone know of a "good" fortran program for inverting matrices?
|> I was looking for a program to use on either a PC or mainframe.  Any help
|> would be very much appeciated.
|> 
|>                                         Thanks
|> 
|>                                         Rob
|> 
|> DEANRE@bscf.uga.edu

What type of matrix: full, sparse, symmetric, banded, 3-diagonal ???
what type of PC or mainframe ???
             Michele

michele.morin@der.edf.fr
Newsgroups: comp.lang.fortran