Fortran passing statement function to function -



Fortran passing statement function to function -

i looking pass statement function evaluated in function, possible? trying have function can numerically integrate general function. there other values in real problem don't want have external function.

edit 2014/11/1 02:22 gmt

there couple of bugs in code @ianh set in. code works after take care of bugs.

program trapezoidz implicit none ! main programme global variables , constants here. phone call main contains subroutine main real :: u, integral1 integer :: integer, parameter :: n = 5000 ! delete declaration of trapezoid_integration , f. integral1 = trapezoid_integration(trigsin, n, 0., 3.14159 / 2.0) write (*,*) "i - trap = ",1 - integral1 end subroutine main function my_f(u) real :: my_f real, intent(in) :: u my_f = (u**4)*exp(u)/((exp(u)-1.0)**2) end function my_f function trigsin(u) result(my_f) real :: my_f real, intent(in) :: u my_f = sin(u) end function trigsin function trapezoid_integration (f, n, start_val, end_val) result (integral) ! delete declaration of integrand (assuming have it). procedure(my_f) :: f integer :: n real :: start_val, end_val real :: integral,u,h integer :: integral = 0.0 i=0,n u = start_val + ((end_val - start_val)*i)/n ! implement eqn g.4 if ((i.eq.0).or.(i.eq.n)) integral = integral+integrand(f, u) else integral = integral+(2.0*integrand(f, u)) end if end h=(end_val - start_val)/n integral = (h/2.0)*integral end function trapezoid_integration function integrand(f, x) result (value) implicit none real :: x real :: value real :: f if (x .lt. 0.00001) x = 0.00001 end if value = f(x) end function integrand end programme trapezoidz

below original code. not use

program trapezoidz implicit none integer, parameter :: n = 10 real :: u, integral integer :: real :: f, trapezoid_integration f(u) = (u**4)*exp(u)/((exp(u)-1.0)**2) integral = trapezoid_integration(f, n, 0.2, 3.0) write (*,*) '#trapezoidal integration = ',integral end programme trapezoidz function trapezoid_integration(f, n, start_val, end_val) implicit none integer :: n real :: start_val, end_val real :: f real :: integral,u,h integer :: integral = 0.0 i=0,n u = start_val + ((end_val - start_val)*i)/n ! implement eqn g.4 if ((i.eq.0).or.(i.eq.n)) integral = integral+integrand(f, u) else integral = integral+(2.0*integrand(f, u)) end if end h=(end_val - start_val)/n integral = (h/2.0)*integral end subroutine trapezoid_integration function integrand(f, x) result (value) implicit none real :: x real :: value real :: f if (x .lt. 0.00001) x = 0.00001 end if value = f(x) end function integrand

it not permitted in standard fortran - statement function not in list of permitted procedure types in constraint c1235 in fortran 2008.

however, constraint indicates, not restricted external functions. since fortran 90 actual argument can module procedure. of fortran 2008 can internal procedure.

that offers possibility of compact , much more robust source:

program trapezoidz implicit none ! main programme global variables , constants here. phone call main contains subroutine main real :: u, integral integer :: integer, parameter :: n = 10 ! delete declaration of trapezoid_integration , f. integral = trapezoid_integration(my_f, 0.2, 3.0) write (*,*) '#trapezoidal integration = ',integral end subroutine main function my_f(u) real :: my_f real, intent(in) :: u my_f = (u**4)*exp(u)/((exp(u)-1.0)**2) end function my_f function trapezoid_integration(f, n, start_val, end_val) ! delete declaration of integrand (assuming have it). procedure(my_f) :: f ... end function trapezoid_integration function integrand(f, x) result (value) ... end function integrand end programme trapezoidz

function fortran

Comments

Post a Comment

Popular posts from this blog

xslt - DocBook 5 to PDF transform failing with error: "fo:flow" is missing child elements. Required content model: marker* -

xslt - DocBook 5 to PDF transform failing with error: "fo:flow" is missing child elements. Required content model: marker* - i've inherited doc publishing process takes docbook , builds html , pdf output, using apache fop both. project had been started, never completed. lots of tweaking, i've able of doc sets build (20+), except handful, , fail severe code 1 of these 2 conditions: "fo:flow" missing kid elements. required content model: marker* "fo:block" not valid kid of "fo:root" the xml output comes out of docbook , docbook-fop checks out well-formed. the fop error in 1 book, "endecagloss", example, points end of string: <fo:flow flow-name="xsl-region-body" start-indent="4pc" end-indent="0pt"/> here’s log output in case: [java] severe: javax.xml.transform.transformerexception: file:/scratch/publishing/hudson/jobs/build-main-endeca-documentation/workspace/serverdoc/w...
Read more

mediawiki - How do I insert tables inside infoboxes on Wikia pages? -

mediawiki - How do I insert tables inside infoboxes on Wikia pages? - i want insert tables in infoboxes. example, have proteinbox i've developed using infobox template (http://health-and-medicine.wikia.com/wiki/template:proteinbox), , fasta field big fit comfortably in proteinbox. hence collapsible table within box. this tried: {{infobox | row 1 = {| class = "mw-collapsible mw-collapsed wikitable" |- | <!-- table content --> |} }} your problem is, putting nested table within template. pipe characters ( | ) in table syntax collides pipe usage in templates. the commonly used hack around create template called template:! , containing pipe character, , utilize when need set tables, parser functions or other stuff using pipe characters, within templates. table (with every | replaced {{!}} ): {{{!}} {{!}}- {{!}} a1 {{!}} b1 {{!}}- {{!}} a2 {{!}} b2 {{!}}} ...the equivalent of {| |- | a1 | b1 |- | a2 | b2 |} furthermore, hav...
Read more

Local Service User Logged into Windows -

Local Service User Logged into Windows - i created windows service tracks app server user connected through load balancer. it retrieves user logged windows on machine, cross references active directory (ad) business relationship homecoming first name, lastly name , e-mail (for notification through website created project) there 1 scheme on network users remote system, vista system, reason returning values "local service" logged machine, it's weird , security guys impossible "local service" business relationship logged windows. here how user logged windows managementscope ms = new managementscope("\\\\.\\root\\cimv2"); objectquery query = new objectquery("select * win32_computersystem"); managementobjectsearcher searcher = new managementobjectsearcher(ms, query); foreach (managementobject mo in searcher.get()) { _username = mo["username"].tostring(); } // remove domain part username string[] usernameparts = _...
Read more