Google Groups no longer supports new Usenet posts or subscriptions. Historical content remains viewable.

Dismiss

Problem with overloaded operations

16 views

Skip to first unread message

arjenmarkus

unread,

Nov 28, 2011, 9:10:58 AM11/28/11

Hello,

the program below produces an ICE with gfortran 4.6 and an error message
that I simply do not comprehend with Intel Fortran 12. The latter is:

pde_objects.f90(375): error #8314: If the rightmost part-name is of
abstract type, data-ref shall be polymorphic [NABLA2]
deriv = diff * solution%nabla2() !+ solution%process() +
solution%source( time )
-----------------------------------^
compilation aborted for pde_objects.f90 (code 1)

I have no idea what it means. If I comment out this line, then
the linker complains about some undefined references:
real_times_obj, obj_plus_obj and obj_assing_obj.

Any suggestions? (I know it is pretty lengthy, but I do not think
I can get it very much shorter, so I have left in its original
form)

Regards,

Arjen

-----
! pde_objects.f90 --
! Solve a diffusion problem using an object-oriented approach
!
! Note:
! This could be turned into a more sophisticated program
! using the techniques described in the chapter on
! mathematical abstractions.
! (That would allow the selection of the time integration
! method in a transparent way)
!

! base_pde_objects --
! Module to define the basic objects
!
module base_pde_objects
implicit none

type, abstract :: base_pde_object
! No data
procedure(compute_process), pointer :: process => null()
procedure(compute_source), pointer :: source => null()
contains
procedure :: initialise
procedure :: nabla2
procedure :: print
procedure(real_times_obj), pass(obj), deferred :: real_times_obj
procedure(obj_plus_obj), deferred :: obj_plus_obj
procedure(obj_assign_obj), deferred :: obj_assign_obj
generic :: operator(*) => real_times_obj
generic :: operator(+) => obj_plus_obj
generic :: assignment(=) => obj_assign_obj
end type

abstract interface
function compute_process( obj )
import base_pde_object
class(base_pde_object), intent(in) :: obj
class(base_pde_object), allocatable :: compute_process
end function compute_process
end interface

abstract interface
function compute_source( obj, time )
import base_pde_object
class(base_pde_object), intent(in) :: obj
real, intent(in) :: time
class(base_pde_object), allocatable :: compute_source
end function compute_source
end interface

abstract interface
function real_times_obj( factor, obj ) result(newobj)
import base_pde_object
real, intent(in) :: factor
class(base_pde_object), intent(in) :: obj
class(base_pde_object), allocatable :: newobj
end function real_times_obj
end interface

abstract interface
function obj_plus_obj( obj1, obj2 ) result(newobj)
import base_pde_object
class(base_pde_object), intent(in) :: obj1
class(base_pde_object), intent(in) :: obj2
class(base_pde_object), allocatable :: newobj
end function obj_plus_obj
end interface

abstract interface
subroutine obj_assign_obj( obj1, obj2 )
import base_pde_object
class(base_pde_object), intent(inout) :: obj1
class(base_pde_object), intent(in) :: obj2
end subroutine obj_assign_obj
end interface

contains

! print --
! Print the concentration field
subroutine print( obj )
class(base_pde_object) :: obj

! Dummy
end subroutine print

! initialise --
! Initialise the concentration field using a specific function
subroutine initialise( obj, funcxy )
class(base_pde_object) :: obj

interface
real function funcxy( coords )
real, dimension(:), intent(in) :: coords
end function funcxy
end interface

! Dummy
end subroutine initialise

! nabla2 --
! Determine the divergence
function nabla2( obj )
class(base_pde_object), intent(in) :: obj
class(base_pde_object), allocatable :: nabla2

! Dummy
end function nabla2

function process_base( obj )
class(base_pde_object), intent(in) :: obj
class(base_pde_object), allocatable :: process_base

! Dummy
end function process_base

function source_base( obj, time )
class(base_pde_object), intent(in) :: obj
real, intent(in) :: time
class(base_pde_object), allocatable :: source_base

! Dummy
end function source_base

end module base_pde_objects

! cartesian_2d_objects --
! PDE object on a 2D cartesian grid
!
module cartesian_2d_objects
use base_pde_objects

type, extends(base_pde_object) :: cartesian_2d_object
real, dimension(:,:), allocatable :: c
real :: dx
real :: dy
contains
procedure :: initialise => initialise_cart2d
procedure :: nabla2 => nabla2_cart2d
procedure :: print => print_cart2d
procedure, pass(obj) :: real_times_obj => real_times_cart2d
procedure :: obj_plus_obj => obj_plus_cart2d
procedure :: obj_assign_obj => obj_assign_cart2d
end type cartesian_2d_object

interface grid_definition
module procedure grid_definition_cart2d
end interface

contains

! grid_definition --
! Initialises the grid
!
subroutine grid_definition_cart2d( obj, sizes, dims )
class(base_pde_object), allocatable :: obj
real, dimension(:) :: sizes
integer, dimension(:) :: dims

allocate( cartesian_2d_object :: obj )
select type (obj)
type is (cartesian_2d_object)
allocate( obj%c(dims(1),dims(2)) )
obj%c = 0.0
obj%dx = sizes(1)/dims(1)
obj%dy = sizes(2)/dims(2)
end select

end subroutine grid_definition_cart2d

! print_cart2d --
! Print the concentration field to the screen
!
subroutine print_cart2d( obj )
class(cartesian_2d_object) :: obj

character(len=20) :: format

write( format, '(a,i0,a)' ) '(', size(obj%c,1), 'f6.3)'

write( *, format ) obj%c

end subroutine print_cart2d

! initialise_cart2d --
! Initialise the concentration field using a specific function
!
subroutine initialise_cart2d( obj, funcxy )
class(cartesian_2d_object) :: obj

interface
real function funcxy( coords )
real, dimension(:), intent(in) :: coords
end function funcxy
end interface

integer :: i, j
real, dimension(2) :: x

obj%c = 0.0
do j = 2,size(obj%c,2)-1
x(2) = obj%dy * (j-1)
do i = 2,size(obj%c,1)-1
x(1) = obj%dx * (i-1)
obj%c(i,j) = funcxy( x )
enddo
enddo

end subroutine initialise_cart2d

! nabla2_cart2d
! Determine the divergence
function nabla2_cart2d( obj )
class(cartesian_2d_object), intent(in) :: obj
class(base_pde_object), allocatable :: nabla2_cart2d

integer :: m, n
real :: dx, dy

m = size(obj%c,1)
n = size(obj%c,2)

dx = obj%dx
dy = obj%dy

allocate( cartesian_2d_object :: nabla2_cart2d )
select type (nabla2_cart2d)
type is (cartesian_2d_object)
allocate( nabla2_cart2d%c(m,n) )

nabla2_cart2d%c = 0.0
nabla2_cart2d%c(2:m-1,2:n-1) = &
-(2.0 * obj%c(2:m-1,2:n-1) - obj%c(1:m-2,2:n-1) -
obj%c(3:m,2:n-1)) / dx**2 &
-(2.0 * obj%c(2:m-1,2:n-1) - obj%c(2:m-1,1:n-2) -
obj%c(2:m-1,3:n)) / dy**2

end select
end function nabla2_cart2d

function real_times_cart2d( factor, obj ) result( newobj )
real, intent(in) :: factor
class(cartesian_2d_object), intent(in) :: obj
class(base_pde_object), allocatable :: newobj

allocate( cartesian_2d_object :: newobj )
select type (newobj)
type is (cartesian_2d_object)
allocate( newobj%c(m,n) )

newobj%c = factor * obj%c
end select
end function real_times_cart2d

function obj_plus_cart2d( obj1, obj2 ) result( newobj )
class(cartesian_2d_object), intent(in) :: obj1
class(base_pde_object), intent(in) :: obj2
class(base_pde_object), allocatable :: newobj

allocate( cartesian_2d_object :: newobj )
select type (newobj)
type is (cartesian_2d_object)
allocate( newobj%c(m,n) )

select type (obj2)
type is (cartesian_2d_object)
newobj%c = obj1%c + obj2%c
end select
end select
end function obj_plus_cart2d

subroutine obj_assign_cart2d( obj1, obj2 )
class(cartesian_2d_object), intent(inout) :: obj1
class(base_pde_object), intent(in) :: obj2

select type (obj2)
type is (cartesian_2d_object)
obj1%c = obj2%c
end select
end subroutine obj_assign_cart2d

end module cartesian_2d_objects

! define_pde_objects --
! Module to bring all the PDE object types together
!
module define_pde_objects
use base_pde_objects
use cartesian_2d_objects

implicit none

interface grid_definition
module procedure grid_definition_general
end interface
contains

subroutine grid_definition_general( obj, type, sizes, dims )
class(base_pde_object), allocatable :: obj
character(len=*) :: type
real, dimension(:) :: sizes
integer, dimension(:) :: dims

select case (type)
case ("cartesian 2d")
call grid_definition( obj, sizes, dims )

case default
write(*,*) 'Unknown grid type: ', trim(type)
stop
end select

end subroutine grid_definition_general

end module define_pde_objects

! pde_specific --
! Module holding the routines specific to the PDE that
! we are solving
!
module pde_specific
implicit none

contains
real function patch( coords )
real, dimension(:), intent(in) :: coords

if ( sum((coords-[50.0,50.0])**2) < 40.0 ) then
patch = 1.0
else
patch = 0.0
endif
end function patch
end module pde_specific

! test_pde_solver --
! Small test program to demonstrate the usage
!
program test_pde_solver

use define_pde_objects
use pde_specific

implicit none

class(base_pde_object), allocatable :: solution, deriv

integer :: i
real :: time, dtime, diff

!
! Create the grid
!
call grid_definition( solution, "cartesian 2d", [100.0, 100.0],
[20, 20] )
call grid_definition( deriv, "cartesian 2d", [100.0, 100.0],
[20, 20] )

!
! Initialise the concentration field
!
call solution%initialise( patch )

!
! Perform the integration - explicit method
!

time = 0.0
dtime = 0.1
diff = 0.1

write(*,*) 'Time: ', time
call solution%print

do i = 1,100

deriv = diff * solution%nabla2() !+ solution%process() +
solution%source( time )
solution = solution + dtime * deriv

if ( mod(i, 10) == 0 ) then
write(*,*) 'Time: ', time
call solution%print
endif
enddo

write(*,*) 'End result: '
call solution%print

end program test_pde_solver

Tobias Burnus

unread,

Nov 28, 2011, 11:16:54 AM11/28/11

On 11/28/2011 03:10 PM, arjenmarkus wrote:
> the program below produces an ICE with gfortran 4.6

It also does with the current 4.7 developer version; I will fill a
bugreport.

> an error message that I simply do not comprehend with Intel Fortran 12.

[...]

> I have no idea what it means. If I comment out this line, then
> the linker complains about some undefined references:
> real_times_obj, obj_plus_obj and obj_assing_obj.

With ifort 12.1.1, I do not get any error message, but I also get
undefined references.

I hope it helps a bit; I have not looked at the code itself. By the way:

Using crayftn 7.1.4.111 it fails with:
ftn-451 crayftn: ERROR NABLA2_CART2D, File = arjen.f90, Line = 227,
Column = 18
This reference to function NABLA2_CART2D without an actual argument
list is illegal.
This looks bogus (the line is a SELECT TYPE with the function-name =
result-name as selector).

And with PGI 11.5-0 it fails with
PGF90-S-0155-Generic set contains a non-specific type bound procedure
for = in type base_pde_object (arjen.f90: 125)

Tobias

Ian Chivers

unread,

Nov 28, 2011, 12:34:07 PM11/28/11

i use the nag compiler. here is the output.

i'll try the others i have access to when i get time.

linux-hn4d:~/document/fortran # nagfor test01.f90 -f2003 NAG Fortran
Compiler Release 5.3(842)
Warning: test01.f90, line 85: Unused dummy variable OBJ
Warning: test01.f90, line 99: Unused dummy procedure FUNCXY
Warning: test01.f90, line 99: Unused dummy variable OBJ
Warning: test01.f90, line 108: Allocatable function NABLA2 has not been
allocated or assigned a value
Warning: test01.f90, line 108: Unused dummy variable OBJ
Warning: test01.f90, line 115: Unused dummy variable OBJ
Warning: test01.f90, line 115: Allocatable function PROCESS_BASE has not
been allocated or assigned a value
Warning: test01.f90, line 123: Unused dummy variable OBJ
Warning: test01.f90, line 123: Allocatable function SOURCE_BASE has not
been allocated or assigned a value
Warning: test01.f90, line 123: Unused dummy variable TIME
Warning: test01.f90, line 251: Symbol M referenced but never set
Warning: test01.f90, line 251: Symbol N referenced but never set
Warning: test01.f90, line 268: Symbol M referenced but never set
Warning: test01.f90, line 268: Symbol N referenced but never set [NAG
Fortran Compiler normal termination, 14 warnings]
linux-hn4d:~/document/fortran # ./a.out
Time: 0.0000000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 1.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 1.000 1.000
1.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 1.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 Runtime Error:
Cannot allocate ALLOCATABLE variable - out of memory Program terminated
by fatal error Aborted

Tobias Burnus

unread,

Nov 28, 2011, 1:18:27 PM11/28/11

On 11/28/2011 05:16 PM, Tobias Burnus wrote:
> On 11/28/2011 03:10 PM, arjenmarkus wrote:
>> the program below produces an ICE with gfortran 4.6
> It also does with the current 4.7 developer version; I will fill a
> bugreport.

Actually, that's a known issue - gfortran doesn't handle nested
operators which return a polymorphic ("CLASS(..)") result. Cf. PR 46328
and PR 51334. I hope we can fix this in the next months, but the current
focus is on supporting polymorphic arrays; hopefully the latter patch
will be ready at some point in this week.

>> an error message that I simply do not comprehend with Intel Fortran 12.
> [...]
>> I have no idea what it means. If I comment out this line, then
>> the linker complains about some undefined references:
>> real_times_obj, obj_plus_obj and obj_assing_obj.
>
> With ifort 12.1.1, I do not get any error message, but I also get
> undefined references.

I have now looked at the source code: I think ifort has a similar issue
as gfortran: If one has
a = b*c()
where "a" is polymorphic, "b" is a polymophic function and "*" and "="
are operators, which return a polymorphic result, the compiler needs to
do something like:
tmp1 = c() ! (1)
tmp2 = tmp1%vptr%mult (tmp1, b) ! (2)
call tmp2%vptr%assign (c, tmp2) ! (3)

(I probably mixed up somewhere the order.) The "vptr" is a pointer to a
derived type which contains the procedure pointers to the function
belonging to the actual type (cf. C++'s virtual function table).

For "b*c()", one can try to write it in one line, which gives:
???%vptr%mult(c(), b)
or
c()%vptr%mult(???, b)
or
c()%vptr%mult(c(), b) ! c() is evaluated twice!
That is, one cannot really write it in a single line as one needs the
type of the "c()" expression twice: Once to decide which mult() function
one actually wants to call and once as argument to that mult() function.

I think ifort solves the problem by replacing the "???" of the first
version with the declared type, which happens to be an abstract
interface. Thus, it fails with a linking error.

gfortran does something similar, except it does not move to the basic
type but has something similar to "???", which causes an internal
compiler error. The solution is - as indicated above - to generate a
temporary variable.

Tobias

PS: At least the ifort part is pure speculation; there are also other
methods possibly, but the following should be a rather common way to
represent a polymorphic variable: Namely, by the actual data (or a
pointer to it) and a pointer to a "table". This "virtual" table
typically contains the default initialization, the size of the
derived-type, procedure pointers to copy or free the data, and procedure
pointers to all the type-bound procedures. Other implementations are
possible, but I think that's the most common one.

Steve Lionel

unread,

Nov 28, 2011, 3:55:41 PM11/28/11

On 11/28/2011 11:16 AM, Tobias Burnus wrote:
> With ifort 12.1.1, I do not get any error message, but I also get
> undefined references.

We know of a problem where ifort sometimes puts out external references
for abstract interfaces. It is being fixed.

--
Steve Lionel
Developer Products Division
Intel Corporation
Merrimack, NH

For email address, replace "invalid" with "com"

User communities for Intel Software Development Products
http://software.intel.com/en-us/forums/
Intel Software Development Products Support
http://software.intel.com/sites/support/
My Fortran blog
http://www.intel.com/software/drfortran

Refer to http://software.intel.com/en-us/articles/optimization-notice
for more information regarding performance and optimization choices in
Intel software products.

arjenmarkus

unread,

Nov 29, 2011, 3:11:47 AM11/29/11

On 2011-11-28 18:34, Ian Chivers wrote:
> i use the nag compiler. here is the output.

Thanks for this. I have solved the errors regarding M and N.
Could you try it with this improved version? (I do not expect
you to solve all the bugs that are lurking in there of course,
I just want to make sure this type of programming is legal)

Regards,

Arjen

implicit none

integer :: m, n

m = size(obj%c,1)
n = size(obj%c,2)

allocate( cartesian_2d_object :: newobj )
select type (newobj)
type is (cartesian_2d_object)
allocate( newobj%c(m,n) )

newobj%c = factor * obj%c
end select
end function real_times_cart2d

function obj_plus_cart2d( obj1, obj2 ) result( newobj )
class(cartesian_2d_object), intent(in) :: obj1
class(base_pde_object), intent(in) :: obj2
class(base_pde_object), allocatable :: newobj

integer :: m, n

m = size(obj1%c,1)
n = size(obj1%c,2)

Rafik Zurob

unread,

Nov 29, 2011, 3:36:50 AM11/29/11

Hi Arjen

I tried your program with the XL Fortran compiler. Here is what I got:

> xlf2003 x.f
"x.f", 1513-083 (E) Internal or module function source_base was not set
within the function.
"x.f", 1513-083 (E) Internal or module function process_base was not set
within the function.
"x.f", 1513-083 (E) Internal or module function nabla2 was not set within
the function.
** base_pde_objects === End of Compilation 1 ===
** cartesian_2d_objects === End of Compilation 2 ===
** define_pde_objects === End of Compilation 3 ===
** pde_specific === End of Compilation 4 ===
** test_pde_solver === End of Compilation 5 ===
1501-510 Compilation successful for file x.f.
> ./a.out
Time: 0.0000000000E+00

Time: 0.0000000000E+00

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.004 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.008 0.988 0.008

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.004 0.988 1.000 0.988
0.004 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.008 0.988 0.008

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.004 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

Time: 0.0000000000E+00

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.008 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.016 0.977 0.016

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.008 0.977 1.000 0.977
0.008 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.016 0.977 0.016

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.008 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

Time: 0.0000000000E+00

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.012 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.023 0.965 0.023

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.012 0.965 0.999 0.965
0.012 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.023 0.965 0.023

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.012 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

Time: 0.0000000000E+00

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.015 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.030 0.954 0.030

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.015 0.954 0.999 0.954
0.015 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.030 0.954 0.030

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.015 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

Time: 0.0000000000E+00

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.019 0.001

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.037 0.943 0.037
0.001 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.019 0.943 0.998 0.943
0.019 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.037 0.943 0.037
0.001 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.019 0.001

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

Time: 0.0000000000E+00

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.022 0.001

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.044 0.933 0.044
0.001 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.022 0.933 0.997 0.933
0.022 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.044 0.933 0.044
0.001 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.022 0.001

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

Time: 0.0000000000E+00

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.025 0.001

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.051 0.922 0.051
0.001 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.025 0.922 0.996 0.922
0.025 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.051 0.922 0.051
0.001 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.025 0.001

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

Time: 0.0000000000E+00

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.029 0.001

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.057 0.912 0.057
0.001 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.029 0.912 0.995 0.912
0.029 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.057 0.912 0.057
0.001 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.029 0.001

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

Time: 0.0000000000E+00

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.002 0.032 0.002

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.002 0.064 0.902 0.064
0.002 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.032 0.902 0.993 0.902
0.032 0.001 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.002 0.064 0.902 0.064
0.002 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.002 0.032 0.002

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

Time: 0.0000000000E+00

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.002 0.035 0.002

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.002 0.070 0.892 0.070
0.002 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.035 0.892 0.992 0.892
0.035 0.001 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.002 0.070 0.892 0.070
0.002 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.002 0.035 0.002

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

End result:

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.002 0.035 0.002

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

Regards

Rafik
Visit the Fortran Cafe at http://ibm.com/rational/cafe

"arjenmarkus" <arjen.m...@gmail.com> wrote in message
news:jb2443$f3$1...@dont-email.me...

arjenmarkus

unread,

Nov 29, 2011, 4:01:17 AM11/29/11

Hi Rafik,

that looks promising - it is the sort of numbers I would expect.
It also confirms that the program will work as envisaged.

Thanks,

Arjen

Ian Chivers

unread,

Nov 29, 2011, 4:49:09 PM11/29/11

Hi Arjen

here is the latest run.
i was working today and didn't see your reply
until tonight

cheers

Ian

linux-hn4d:/home/ian # nagfor test02.f90 -f2003 NAG Fortran Compiler
Release 5.3(842)
Warning: test02.f90, line 85: Unused dummy variable OBJ
Warning: test02.f90, line 99: Unused dummy procedure FUNCXY
Warning: test02.f90, line 99: Unused dummy variable OBJ
Warning: test02.f90, line 108: Allocatable function NABLA2 has not been

allocated or assigned a value

Warning: test02.f90, line 108: Unused dummy variable OBJ
Warning: test02.f90, line 115: Unused dummy variable OBJ
Warning: test02.f90, line 115: Allocatable function PROCESS_BASE has not

been allocated or assigned a value

Warning: test02.f90, line 123: Unused dummy variable OBJ
Warning: test02.f90, line 123: Allocatable function SOURCE_BASE has not

been allocated or assigned a value

Warning: test02.f90, line 123: Unused dummy variable TIME [NAG Fortran
Compiler normal termination, 10 warnings] linux-hn4d:/home/ian #
inux-hn4d:/home/ian # ./a.out

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.004

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.008 0.988
0.008 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.004 0.988 1.000

0.988 0.004 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.008 0.988
0.008 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.004

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.008

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.016 0.977
0.016 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.008 0.977 1.000

0.977 0.008 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.016 0.977
0.016 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.008

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.012

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.023 0.965
0.023 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.012 0.965 0.999

0.965 0.012 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.023 0.965
0.023 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.012

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.015

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.030 0.954
0.030 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.015 0.954 0.999

0.954 0.015 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.030 0.954
0.030 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.015

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.019
0.001 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.037 0.943
0.037 0.001 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.019 0.943 0.998

0.943 0.019 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.037 0.943

0.037 0.001 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.019
0.001 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.022
0.001 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.044 0.933
0.044 0.001 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.022 0.933 0.997

0.933 0.022 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.044 0.933

0.044 0.001 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.022
0.001 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.025
0.001 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.051 0.922
0.051 0.001 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.025 0.922 0.996

0.922 0.025 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.051 0.922

0.051 0.001 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.025
0.001 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.029
0.001 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.057 0.912
0.057 0.001 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.029 0.912 0.995

0.912 0.029 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.057 0.912

0.057 0.001 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.029
0.001 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.002 0.032
0.002 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.002 0.064 0.902
0.064 0.002 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.032 0.902 0.993

0.902 0.032 0.001 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.002 0.064 0.902

0.064 0.002 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.002 0.032
0.002 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.002 0.035
0.002 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.002 0.070 0.892
0.070 0.002 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.035 0.892 0.992

0.892 0.035 0.001 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.002 0.070 0.892

0.070 0.002 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.002 0.035
0.002 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

End result:

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

0.892 0.035 0.001 0.000 0.000 0.000 0.000 0.000 0.000

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.002 0.070 0.892

0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

Warning: Floating underflow occurred
linux-hn4d:/home/ian #

arjenmarkus

unread,

Nov 30, 2011, 2:45:49 AM11/30/11

On 2011-11-29 22:49, Ian Chivers wrote:
> Hi Arjen
>
> here is the latest run.
> i was working today and didn't see your reply
> until tonight
>

Hi Ian,

thank you very much. I won't speculate on the sudden
linebreaks in the middle of the matrices, but I do
see one flaw: the simulation time is not updated ;).

Like I replied before, this makes me think this
type of programming can be done - it makes
for a very concise style (once you get beyond the
ground layer of the derived types' implementation).

Regards,

Arjen

0 new messages