passing assumed-shape arrays in two levels of subroutines (Fortran 90) -

-

i have had problems calling successive subroutines assumed-shape arrays in fortran 90. more specifically, call 2 levels of subroutines, passing assumed-shape array parameter, in end array lost. demonstrate it, 1 can follow code below.

  program main    interface      subroutine sub1(x)        real, dimension(:):: x        real c      end subroutine sub1       subroutine sub2(x)        real, dimension(:):: x        real c      end subroutine sub2   end interface    real, dimension(:), allocatable:: x    allocate(x(1:10)) ! first executable command in main   x(1) = 5.   call sub1(x)   write(*,*) 'result = ',x(1)   deallocate(x)   end program main    subroutine sub1(x) ! first subroutine   real, dimension(:):: x   real c   call sub2(x)   end subroutine sub1    subroutine sub2(x) ! second subroutine   real, dimension(:):: x   real c   c2=x(1)   end subroutine sub2

very shortly, main allocates x call sub1(x). sub1 calls sub2(x). means allocated array passed subroutine passes subroutine. expect have in sub2 same array i've created in main, no. using gdb tool explore it, this:

1) in main, before calling sub1, array x defined:

(gdb) p x
$1 = (5, 0, 0, 0, 0, 0, 0, 0, 0, 0)

2) within sub1, before calling sub2, x defined:

(gdb) p x
$2 = (5, 0, 0, 0, 0, 0, 0, 0, 0, 0)

3) inside sub2, however, x have unexpected value , dimension absolutely wrong:

(gdb) p x
$3 = ()
(gdb) whatis x
type = real(4) (0:-1)

so, x has been passed main sub1, not sub1 sub2. i've been using intel fortran gfortran same results.

i've struggling long time. appreciated.
g.oliveira.

the use of assumed-shape dummy arguments requires explicit interface. in main program you've provided explicit interfaces 2 subroutines, these don't propagate subroutines themselves. subroutines compiled separate units, if you've put code 1 source file.
means sub1 doesn't have explicit interface available sub2, , uses implicit interface, argument x assumed real scalar.

all avoided putting 2 subroutines in module , use module in main program, automatically making explicit interfaces available. way don't have provide interfaces yourself, error prone.

as side note, advise use of implicit none in code.


Comments

Post a Comment

Popular posts from this blog

Javascript line number mapping -

-
magically formatted comments change reported line number of javascript errors in browsers; this: //@line n "f" n line number , f file name. unfortunately, //@line appears ungoogleable. know there documentation on feature, , browsers support it? (i found references here , here .) from can find available in mozilla based browsers. it using js_setoptions command can turn these options on or off depending on values passed. https://developer.mozilla.org/en/js_setoptions
Read more

c# - Is it possible to remove an existing registration from Autofac container builder? -

-
something along lines: builder.registertype<mytype>().as<itype>(); builder.registertype<mytype2>().as<itype>(); builder.deregistertype<mytype>().as<itype>() var container = builder.build(); var types = container.resolve<ienumerable<itype>>(); assert.istrue(types.count == 1); assert.istrue(types[0].gettype == typeof(mytype2)); scenario: go through bunch of assemblies , go register types want make sure have 1 implementation of given type. need before create container. track on own nice if autofac me bit. this cannot done directly using containerbuilder , unless start on new one. mind you, having first built container should able construct new container filtering away unwanted types , reusing registrations first container. this: ... var container = builder.build(); builder = new containerbuilder(); var components = container.componentregistry.registrations .where(cr => cr.activator.limittype != typ...
Read more

php - Mysql PK and FK char(36) vs int(10) -

-
i have read mysql can handle indexes , searches better if tables keys integers rather char(36) uuid's. is worth converting int(10)? note: use 1 server , have have plans use multiple databases concurrently therefore removing 1 of big reasons using uuid. our tables innodb if makes difference. thanks in advance. usually rdbms can handle integer keys pk more efficient, other datatypes. reason how build index column, yes: long dont require string (or other typed) keys, should use integers. however: char(36) , int(10) far away being equal, because int(10) much smaller, char(36). dont know, if require many different keys, should keep in mind. update, complete last paragraph: int(10) 32-bit, char(36) 36 - ^byte^ (= 288-bit). not mean, int consumes less space, means, char(36) has 4 times more different keys.
Read more