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

linux - Mailx and Gmail nss config dir -

-
i trying send mail linux command line using mailx command. can send local domain no problem want set mail send gmail account receive mail sent gmail account. after configuring mail.rc so: account gmail { set smtp=smtps://smtp.gmail.com:587 set smtp-auth=login set smtp-auth-user=username@gmail.com set smtp-auth-password=password set ssl-verify=ignore } i error: resolving host smtp.gmail.com . . . done. connecting 74.125.25.109 . . . connected. missing "nss-config-dir" variable. "/home/username/dead.letter" 11/354 . . . message not sent. after looking "nss-config-dir" here , located certn.db , keyn.db files , added mail.rc so: account gmail { set smtp=smtps://smtp.gmail.com:587 set smtp-auth=login set smtp-auth-user=username@gmail.com set smtp-auth-password=password set ssl-verify=ignore set nss-config-dir=/home/user/.mozilla/firefox/location.default } now when try send mail using command: echo "sent gmail account" | mailx -v -a
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