operator overloading - Fortran: Accessing values from the return type of a function -


i have created derived type v, , included add function can use operator +.

however when 

z = u + v

the operation not performed. think because z%kn not being accessed.

however when 

call vsum(z, u, v)

everything works expected

below declaration of derived type , overloading function vadd.

module vtest  type :: v    character (len=8) :: kn   real,  allocatable :: vc(:)    contains      procedure :: vadd     generic :: operator (+) => vadd   end type vtest  contains  function vadd (b, c) result (a)   type (v) ::   class (v), intent (in) :: b, c    !!$ in vsum, use made of a% kn    call vsum (a, b, c)  end function vadd   subroutine vsum (ta, tb, tc)    type (v), intent (inout)  :: ta    type (v), intent (in) :: tb, tc     logical :: la, lb, lc     la = .false.;  lb = .false.;  lc = .false.     select case (ta%kn)      case ("real32")      if (allocated (ta%vc)) la = .true.      if (allocated (tb%vc)) lb = .true.      if (allocated (tc%vc)) lc = .true.      if (la .and. lb .and. lc)         ta%vc = tb%vc + tc%vc      end if     end select   end subroutine vsum  end module vtest   program test use vtest type (v) :: z     z% kn = "real32"    allocate (z% vc_real32(3))    write (*,*) "z = u + v"    write (*,*) "z% kn: ", z% kn    z = u + v    write (*,*) "z% kn: ", z% kn    write (*,*) "z: ", z% vc_real32  end program vtest

let's @ subroutine vsum. work done in subroutine depends on definition of components of 3 objects ta, tb , tc. summation proceed expect necessary allocatable components allocated , match of ta%kn 'real32'.

there's no complete working example, things work when there call like

call vsum(z, u, w)  ! v name of type, call variable w

z%kn set 'real32'.

however, defined operation

z = u + w

there function reference

z = vadd(u, w)

in there call

call vsum (a, b, c)  ! actual arguments "z", u, w.

now, a function result of vadd. consequent of that, dummy argument intent(out) a undefined.

that means, essentially, when passed through vadd z makes vsum undefined component kn , unallocated component vc. so, requirements summation not met (and indeed select case forms invalid reference).

i suppose can comment on how fix issue, also.

in vsum, you're taking, in effect, type of operation "result". instead, like

subroutine vsum (ta, tb, tc)   type (v), intent (out)  :: ta   type (v), intent (in) :: tb, tc    logical :: lb, lc    lb = .false.;  lc = .false.    ! in here want consistency checking...   select case (tb%kn)       case ("real32")     if (allocated (tb%vc)) lb = .true.     if (allocated (tc%vc)) lc = .true.     if (lb .and. lc)        ta%vc = tb%vc + tc%vc   ! ta%vc being automatically allocated     end if   end select   end subroutine vsum

this, naturally, doesn't address seeming desire (from title) take type of operation result. could original subroutine approach, except in 1 instance isn't philosophy of fortran functions.

if want use defined operations, can't have operation defined left-hand side of assignment. in expression u+w there isn't left-hand side, defined operation still expected behave. say, in statement

z = u + w

you aren't saying "z result of operation + applied u , w". but: "the operation + applied u , w evaluated , result (through defined or implicit assignment) assigned z". there assignment before "type of result" reached. why earlier on put z actual argument in quotes.


-

Comments

Post a Comment

Popular posts from this blog

scala - 'wrong top statement declaration' when using slick in IntelliJ -

Image
i'm porting project scala 2.9 2.10, therefore have use slick instead of scalaquery. i'm using slick 2.1.0 since supports ms access. according this tutorial , upgrade guide changed robs object class , added val robs : why error message wrong top statement declaration , how rid of it? edit: i'm new scala... according this question seems can't put val outside of methods or classes, right? code above directly in packge. right approach slick then? should move code in class or trait? according this changed val robs = tablequery[robs] to object robs extends tablequery(new robs(_)) {} . no warnings or errors anymore. :)
Read more

c# - DevExpress.Wpf.Grid.InfiniteGridSizeException was unhandled -

Image
when using devexpress, see error: devexpress.wpf.grid.infinitegridsizeexception unhandled message="by default, infinite grid height not allowed since grid rows rendered , hence grid work slowly. fix issue, should place grid container give finite height grid, or should manually specify grid's height or maxheight. note can avoid exception setting gridcontrol.allowinfinitegridsize static property true, in case grid run slowly." the problem dxgrid has infinite height. to fix, set height non-infinite. snoop absolutely invaluable this: if "height" xaml element infinite (i.e. 0 or nan ), can set using one of following: option 1: height="{binding path=actualheight, relativesource={relativesource mode=findancestor, ancestortype=uielement}}" option 2: verticalalignment="stretch" option 3: height="auto" hint: use verticalalignment="stretch" if child of grid <rowdefinition height=
Read more

PySide and Qt Properties: Connecting signals from Python to QML -

i trying run simple pyside appliccation design made in qml. call function in qml python in separated thread. i have following qml view: import qt 4.7 rectangle { visible: true width: 100; height: 100 color: "lightgrey" text { id: datestr objectname: "datestr" x: 10 y: 10 text: "init text" font.pixelsize: 12 function updatedata(text) { datestr.text = qstr(text) console.log("you said: " + text) } } } here application: some imports.... #!/usr/bin/python import sys import time pyside import qtdeclarative, qtcore pyside.qtgui import qapplication pyside.qtdeclarative import qdeclarativeview the qthread worker.... #worker thread class worker(qtcore.qthread): updateprogress = qtcore.signal(int) def __init__(self, child): qtcore.qthread.__init__(self) self.child = child def run(self):
Read more