Adding type checking helps

Preben Alsholm — Tue, 22 Nov 2011 01:39:02 Z

If you add type checking Compile works for me (both in Maple 12 and 15).

ggg:=(y::float,a::float,b::posint)->
(2*(-a+6*y^2*a^2+2*b*y^2*a^2+9*b*y^3*a .........

g:=Compiler:-Compile(ggg);

g(1.23,.89,2);

#That resulted in 6.39854925121287010

typed arguments

pagan — Tue, 22 Nov 2011 01:49:05 Z

You should declare the arguments to your procedure with typechecking. You might want to change the type checks for `a` and `b` below from float to integer (and then pass them as such, accordingly, of course).

restart:
ggg_expression:=(2*(-a+6*y^2*a^2+2*b*y^2*a^2+9*b*y^3*a+3*b*y^2*a
+b*y^6*a^2+4*b*y^5*a^2+5*b*y^4*a^2+2*b*y^5*a+8*b*y^4*a+3*b*y^3*a^2
+4*b*y^2-y^3*a-3*y^2*a+y^6*a^3+5*y^5*a^3+8*y^4*a^3+4*y^3*a^3+y^5*a^2
+5*y^4*a^2+9*y^3*a^2+b*y^4+4*b*y^3+y*a^2+2*a*b*y+y^3*a*(y+1)^(-2*b)
+a*(y+1)^(-2*b)+3*y^2*a*(y+1)^(-2*b)+3*y*a*(y+1)^(-2*b)
-2*y^2*a^2*(y+1)^(-2*b)-y*a^2*(y+1)^(-2*b)-3*y*a-y^3*a^2*(y+1)^(-2*b))
/(y^2*(a+b)*(y+1)*(y+2)^2*(y*a+1)^2)+(-2-y^2+2*y^2*a^2-2*y+4*b*y^3*a
+8*b*y^2*a+2*b*y^4*a^2+4*b*y^3*a^2+2*b*y^2+4*b*y+2*y^2*a+y^4*a^2
+4*y^3*a^2+(y+1)^(-2*b)-2*y^2*a*(y+1)^(-2*b)-2*y*a*(y+1)^(-2*b)
+y^2*a^2*(y+1)^(-2*b)+(y+1)^(-2*b)*y^2+2*(y+1)^(-2*b)*y+(y+1)^(-2*b-2)*y^2
+2*(y+1)^(-2*b-2)*y+(y+1)^(-2*b-2)*y^4*a^2+2*(y+1)^(-2*b-2)*y^3*a
+2*(y+1)^(-2*b-2)*y^3*a^2+4*(y+1)^(-2*b-2)*y^2*a+(y+1)^(-2*b-2)*y^2*a^2
+2*(y+1)^(-2*b-2)*y*a+(y+1)^(-2*b-2))/(y^2*(a+b)*(y+2)^2*(y*a+1)^2))
/((-a+6*y^2*a^2+2*b*y^2*a^2+9*b*y^3*a+3*b*y^2*a+b*y^6*a^2+4*b*y^5*a^2
+5*b*y^4*a^2+2*b*y^5*a+8*b*y^4*a+3*b*y^3*a^2+4*b*y^2-y^3*a-3*y^2*a
+y^6*a^3+5*y^5*a^3+8*y^4*a^3+4*y^3*a^3+y^5*a^2+5*y^4*a^2+9*y^3*a^2
+b*y^4+4*b*y^3+y*a^2+2*a*b*y+y^3*a*(y+1)^(-2*b)+a*(y+1)^(-2*b)
+3*y^2*a*(y+1)^(-2*b)+3*y*a*(y+1)^(-2*b)-2*y^2*a^2*(y+1)^(-2*b)
-y*a^2*(y+1)^(-2*b)-3*y*a-y^3*a^2*(y+1)^(-2*b))^2
/(y^4*(a+b)^2*(y+1)^2*(y+2)^4*(y*a+1)^4)-(1/4)*(-2-y^2+2*y^2*a^2-2*y
+4*b*y^3*a+8*b*y^2*a+2*b*y^4*a^2+4*b*y^3*a^2+2*b*y^2+4*b*y+2*y^2*a
+y^4*a^2+4*y^3*a^2+(y+1)^(-2*b)-2*y^2*a*(y+1)^(-2*b)-2*y*a*(y+1)^(-2*b)
+y^2*a^2*(y+1)^(-2*b)+(y+1)^(-2*b)*y^2+2*(y+1)^(-2*b)*y+(y+1)^(-2*b-2)*y^2
+2*(y+1)^(-2*b-2)*y+(y+1)^(-2*b-2)*y^4*a^2+2*(y+1)^(-2*b-2)*y^3*a
+2*(y+1)^(-2*b-2)*y^3*a^2+4*(y+1)^(-2*b-2)*y^2*a+(y+1)^(-2*b-2)*y^2*a^2
+2*(y+1)^(-2*b-2)*y*a+(y+1)^(-2*b-2))^2/(y^4*(a+b)^2*(y+2)^4*(y*a+1)^4)):

ggg:=unapply(ggg_expression,[a::float,b::float,y::float]):

cggg:=Compiler:-Compile(ggg):

ggg(1.0,1.0,3.2);

                          9.924398803

cggg(1.0,1.0,3.2);

                      9.92439880079097314

N:=10000:

V1:=CodeTools:-Usage( Vector(N,(i)->ggg(1.0,1.0,2.0+1.0/i),datatype=float[8]) ):
memory used=0.88GiB, alloc change=39.37MiB, cpu time=20.44s, real time=20.62s

V2:=CodeTools:-Usage( Vector(N,(i)->evalhf(ggg(1.0,1.0,2.0+1.0/i)),datatype=float[8]) ):
memory used=1.90MiB, alloc change=0 bytes, cpu time=891.00ms, real time=899.00ms

V3:=CodeTools:-Usage( Vector(N,(i)->cggg(1.0,1.0,2.0+1.0/i),datatype=float[8]) ):
memory used=4.35MiB, alloc change=0 bytes, cpu time=187.00ms, real time=205.00ms

LinearAlgebra:-Norm(V1-V2), LinearAlgebra:-Norm(V2-V3);

                            -8                        -10
      3.23561790693815965 10  , 9.92976367797382409 10

good

icegood — Tue, 22 Nov 2011 02:27:01 Z

completly good! One problem less!

MaplePrimes - answers and comments on Question, doesnt compiled

Adding type checking helps

typed arguments

good