Hello, I have read several posts about this problem in this and other
newsgroups, but I can't quite wrap my head around a solution that fits
my problem.

I have a program with the basic skeleton:
Main Program
Call FileReader(one)
Call FileReader(two)
....
Call FileWriter
End Program

Main Program has 2 allocatable arrays (one and two), the sizes of
which is unknown until they are passed to FileReader. I need
FileReader to allocate the array "one" and fill it with data, and send
the array back to Main Program. FileReader will do the same with a
second file and array "two". From there, Main Program will call
FileWriter using the 2 allocatable arrays.

Does the BLAS library has the capicity to do two band matrixes
product?

And which function did it?

Hi, all
I have heard that the build-in function anint of fortran is kind
of slow, but I have to used it quite frequent in my code. Anyone know
a more efficient version of it?

Thanks a lot.

http://www.youtube.com/watch/v/as_b6...&feature=email

Lecturer: Bro Sheikh Feiz 4rm Austrailia

For how long, for how long are we going to be heedless of the signs of
Allah (SWT), for how long are we going to be busy with the Dunya
(World), busy making money, chasing our dreams, following our desires...

FREE SOFTWARE DOWNLOAD
http://freesoftware20.blogspot.com

Hi all,

Hi all,

Sorry for bother again, this time, I attach the fortran binary data
file I want to read use the fortran code.

I've a unformatted binary file named Si.cst_esp which include the
data as follows if these data have been read out into a regular format
file:

---------------
1 1 1 4.185302
2 1 1 2.377237
3 1 1 -0.053069
4 1 1 -1.013479
5 1 1 ...

On Jul 5, 3:18 pm, Damian rouson.net> wrote:

Dear all;

I'm going to use rlft3.f (3D Fast Fourier transform subroutine in
Numerical Recipes book )
I have constructed a 107*107*107 mesh for data.
I should call rlft3(data,speg,nn1,nn2,nn3,isign).
-----------------------------------------------------------------------------
SUBROUTINE rlft3(data,speq,nn1,nn2,nn3,isign)
INTEGER isign,nn1,nn2,nn3
COMPLEX data(nn1/2,nn2,nn3),speq(nn2,nn3)
-----------------------------------------------------------------------------
1- I don't understand why for data is written nn1/2 ? I have a
107*107*107 mesh!!!
2- Is it true to put nn1=nn2=nn3=107
3- I don't know what should I put instead of speq ?

I'm anxiously looking forward your reply and guidlines.
With my best regards,
Fatemeh

Dear all;

I'm going to use rlft3.f (3D Fast Fourier transform subroutine in
Numerical Recipes book )
I have constructed a 107*107*107 mesh for data.
I should call rlft3(data,speg,nn1,nn2,nn3,isign).
-----------------------------------------------------------------------------
SUBROUTINE rlft3(data,speq,nn1,nn2,nn3,isign)
INTEGER isign,nn1,nn2,nn3
COMPLEX data(nn1/2,nn2,nn3),speq(nn2,nn3)
-----------------------------------------------------------------------------
1- I don't understand why for data is written nn1/2 ? I have a
107*107*107 mesh!!!
2- Is it true to put nn1=nn2=nn3=107
3- I don't know what should I put instead of speq ?

I'm anxiously looking forward your reply and guidlines.
With my best regards,
Fatemeh

The regular fireworks are over. Let me see if I can light
up a few of my own. Following (in separate articles) are
the public review comments I'm sending in. I've tried
to keep the number down, and I've tried to keep the
comments short. I'll be sending them in sometime
Sunday.

--
J. Giles

"I conclude that there are two ways of constructing a software
design: One way is to make it so simple that there are obviously
no deficiencies and the other way is to make it so complicated
that there are no obvious deficiencies."