Items tagged with compile compile Tagged Items Feed

Hi all,

I'm currently running a parameter space search via Maple scripting. The pseudo-code goes something like this:

fopen text file to write data to;

initialize parameters (variable) x, y, and z;

for x number of iterations

   A:-Simulate(model.msim, with [parameters]);

   Update parameters with new values;

   Write to text file;

end do;

fclose(text file)

Wide set of expressions can convert to compiled functions. Expressions can even include definite integrals.
Hope it helps for others who want really speedup calculations in maple as much as possible for now.

ex.mw

Checked under 15.01 version

How to do subj? Straightforward way doesn't compile:

restart;

myint:= proc(u::float)::float;
return int(exp(x^6), x=0.0..u);
end proc:


createcomp:=proc (expr::evaln, exprfunc::evaln);
exprfunc:=(v::float)->unapply(evaln(expr),a)(v);
return Compiler:-Compile(exprfunc,optimize);
end proc:
global_expr:=myint(a)+a^a;
s:=createcomp(global_expr, global_exprfunc);

"Error, (in Print) cannot translate value of lexically...

ggg:=(y,a,b)->
(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...

tryed to compile samples\OpenMaple\cmaple\omexample.c in MS Visual C++ Express 2010. Everything is OK but after running is writes "The program can't start because maplec.dll is missing from your computer..." Why? I've linked maplec.lib statically! Besides \Maple\bin.win\ in my library path where both maplec.dll and maplec.lib are. System - Win7 Home Premium +SP1

And so with this provocative title, "pushing dsolve to its limits" I want to share some difficulties I've been having in doing just that. I'm looking at a dynamic system of 3 ODEs. The system has a continuum of stationary points along a line. For each point on the line, there exist a stable (center) manifold, also a line, such that the point may be approached from both directions. However, simulating the converging trajectory has proven difficult.

I have simulated as...

Using the Open Maple application program interface (API), the ExternalCalling module, a C compiler, and some basic understanding of the Maple data structure format, you can create Maple utilities that run nearly as fast as Maple builtin procedures.

To show how this is done, I will start with a simple example and then proceed to a slightly more complex procedure that improves ?ListTools[SearchAll]

The Maple ?Compiler can compile a limited subset of Maple commands to native code. The result is substantially faster than running interpreted Maple code. This article shows how you can save a compiled procedure for reuse.

When a Maple procedure is compiled, a shared object library (dynamic-linked library for Windows) is created on the file system, as a temporary file.  When the Maple...

Hello there,  How can I compile this? I used   boundarycompiled := Compiler:-Compile(findboundary)    but this didnt work. Thanks in advance

> findboundary:=proc(a,b,c,d,e,f,g,h)  ;
> global A,B,lowerbound,upperbound;
> A, B := solve(a = b, x), solve(c = b, x);
> x1a, x1b := evalf(solve((x-e)^2+(a-f)^2-g^2)); if x1a < x1b then minx1, maxx1 := x1a, x1b else minx1, maxx1 := x1b, x1a end if;
> x2a, x2b := evalf(solve((x-e...

For each of the definitions

p := proc( x :: float ) :: float; x^(2/3) end proc:
p := proc( x :: float ) :: float; x^evalf((2/3)) end proc:

i get an error by calling Compiler:-Compile( p ):

"Error, (in printtab[CodeGeneration:-Names:-Rational])
rational numbers and arithmetic are not yet supported"

Any suggestions what to do (except manually converting
all exponents to floats)?

Edited to add: for example the following does not work:

c := proc (S::float,...