MaplePrimes Posts

MaplePrimes Posts are for sharing your experiences, techniques and opinions about Maple, MapleSim and related products, as well as general interests in math and computing.

Latest Post
  • Latest Posts Feed
  • This post will explain how to configure the compiler and other tools that will be necessary for you to build the External Calling examples that will come in later posts.  This is an advanced topic and so this post is fairly complex.

    First, I am going to be using the compilers via the command line, so you will need to familarize yourself with the terminal program on your particular OS.  You'll have to do this for yourself, but here are a few starting points:

    Windows

    Apple

    I am going to assume that Linux and Solaris users are familar with using the terminal.

    For Linux, Apple and Solaris, I am going to use gcc as the compiler.  For Linux you should use your distribution's package management system to get it, for Apple you need to install Xcode and for Solaris, well, gcc is probably already installed or you'll want to talk to you sys admin to have it installed (or if you are your own sys admin, you probably know how to install gcc for yourself).  For Windows, you need to install the Windows Software Development Kit.  If you already have a copy of Visual Studio C++ (Express or Professional) installed, then you already have these tools.

    I am also going to use the "make" program to manage the building of the examples, thus you will need to install a version of make as well (you won't need to learn how make works unless you want to modify the examples).  I will be using gnu make, which should be easy to install on Linux and Solaris (similar to how you installed gcc) and it is included in Xcode for Apple.  For Windows, use this:

    http://gnuwin32.sourceforge.net/packages/make.htm

    Installing 32 bit make on 64 bit windows is fine.

    Now you'll need to launch a terminal.  For Linux, Apple and Solaris this should be easy, on Windows go to the Windows SDK folder (or Windows Visual Studio folder) on the Start menu, there should be an icon for Windows SDK Command Prompt.  Click that to launch the terminal.  This version of the terminal has the environment configured to run the compiler.

    On Windows you'll also have to add the location you installed make to your path, which can be done on 32 bit windows like this:

    path=%PATH%;C:\Program Files\GnuWin32\bin

    and on 64 bit Windows like this:

    path=%PATH%;C:\Program Files (x86)\GnuWin32\bin

    assuing you used the default install location for make.

    You can test this by running "make" in the terminal.  If everything is set up correctly, make should run but not find a Makefile and it will raise an error.  If the path is not set properly, make won't be found you'll get a message saying that.

    Path not set properly:

    C:\Program Files\Microsoft SDKs\Windows\v7.1>make
    'make' is not recognized as an internal or external command,
    operable program or batch file.

    Set the path:

    C:\Program Files\Microsoft SDKs\Windows\v7.1>path=%PATH%;"C:\Program Files (x86)
    \GnuWin32\bin"

    Make is now found, but there is no makefile in the current directory

    C:\Program Files\Microsoft SDKs\Windows\v7.1>make
    make: *** No targets specified and no makefile found.  Stop.

    As a final test, I've attached a small example (test.zip) that contains a Makefile and a simple source file.  If you extract the files to a new directory, go to that new directoy in the terminal and run make (with make added to the path as described above) it should build an executable (test or test.exe).  You can run the executable by executing "test" on the command line.

    By default the Makefile is configured for Windows, so Windows users won't need to change it, however other users will need to comment out the

    WINDOWS=true

    line in Makefile by changing it to

    #WINDOWS=true

    I know this is a little confusing, especially if you are not familar with the command line interface, therefore I encourage you post replies if you have problems.  Hopefully we will be able to answer your questions.  Once everyone has figured out how to get this simple example to compile and run on their system, the upcoming external calling examples will be (relatively) easy.

    Good Luck!

    Darin

    Maple T.A.  7 and Maple T.A. MAA Placement Test Suite 7 are now available.

    Maple T.A. 7 provides new options for analyzing grades and gaining a deep understanding of how students are performing, including:

    • Instructors can view all responses to an assignment question, to easily look for patterns
    • The grading scheme for the entire course can be defined inside Maple T.A., with appropriate weightings and flexible policies for dealing with missed, repeated, and worst assignments...

    For the last few days I have been getting as many a 4 email notifications for a single update to a thread in which I am participating.

    A prospective customer recently asked if we had a MapleSim model of a double pipe heat exchanger. Heat exchangers are a critical unit operation in the process industries, and accurate models are needed for process control studies.  I couldn't find an appropriate model so I decided to derive the dynamic equations, and implement them using MapleSim's custom component interface.  I'll outline my modeling strategy in this blog post.

    I've had a few request to provide some more information on External Calling, so I thought I would make a few posts about it. This first post will be a high level description of External Calling and how it works, with examples coming later. As External Calling is an advanced topic, I am going to assume you know how to compile a shared library and are generally familiar with the C language. Although this first post won't require any real programming knowledge.

    What is External Calling?

    External Calling is the name for Maple's ability to connect to and call functions from other programming languages. Maple uses this for various reasons. We have written our own libraries in C, C++ and Java to solve particular problems. We partner with various labs around the world who have developed code, often in languages like C or C++, so external calling is used to interface with their code. We also connect to high performance libraries like NAG and BLAS to provide those high performance routines in Maple. Of course, you can use External Calling to connect Maple to your code as well.

    Although Maple can call various programming languages, the most common languages we connect to are C and C++, and those are the languages I am going to focus on.

    How does it work?

    In Maple, you call ?define_external or use the ?ExternalCalling package. Both these methods take a description of the function that you want to call and returns a Maple procedure. Normally you would assign the procedure to a name and then call the externally defined function just like any other Maple procedure.

    There are a couple different ways to use define_external to connect to a shared library, the differences are mostly concerned with how the parameters given in Maple are converted to parameters used in the external function.

    • Wrapperless external calling. With wrapperless external calling, Maple calls a function implemented in the shared library by automatically converting the values given in Maple into valid types for the external function.
    • Generated wrappers: With generated wrappers, Maple automatically generates a small C library that handles conversions from Maple values to the values used in the external function. Using generated wrappers allows Maple to handle more data types, like call back procedures.
    • Custom wrappers: A custom wrapper is a C function that you write yourself. This function accepts arguments as Maple data structures and returns a Maple data structure. You are responsible for converting the Maple data structures into whatever forms you need and converting your computed value back into a Maple data structure. Maple provides the External Calling API to assist in working with Maple from the externally defined function.

    The first two forms of external calling are the easiest to do, however they are also the most limited. Internally we exclusively (I think) use the third, custom wrapper, form of external calling. That is the form I am going to talk about.

    Custom Wrapper

    The name "custom wrapper" is a bit of a misnomer. The function that you write does not need to "wrap" anything, it can implement anything you want. As long as you can convert the result into a Maple data structure, you can pass it back into Maple. In fact Maple also supports returning generic data, via the ?MaplePointer routines, but that is a more complex topic for a later blog post.

    Your external function is simply a C function with the following calling convention:

    ALGEB CustomWrapper( MKernelVector kv, ALGEB args )

    ALGEB is the C data type that represents a Maple data structure. The MKernelVector is a data structure that acts as an intermediary between your external calling routines and the Maple engine. You will need to pass this structure back into the External Calling API functions. Both of these types, plus the External Calling API functions are defined in a header, maplec.h, that needs to be included in your code. I will provide more details when I provide examples.

    The External Calling API

    The External Calling API is a set of functions that we make available for working with the Maple Engine from external code. Maple also allows third party applications to load the Maple engine as a shared library, we call this ?OpenMaple. The External Calling functions are also available in OpenMaple, so you will often see OpenMaple used in the Maple help pages. Most functions can be used in both OpenMaple and External Calling, except for a few that are OpenMaple specific and involve starting and stopping the Maple Engine.

    Maple's help system documents all the External Calling functions so you can see what is available. There is an overview of the external calling functions on this page, ?ExternalCalling,C,API. Briefly, however there are functions for converting Maple types to C and back, creating and interacting with Maple data structures (list, set, rtable, table, string, etc), creating and interacting with Maple language elements (names, procedures, etc), printing to the Maple interface, memory allocation, evaluating Maple statements and raising exceptions. There is even a C interface to the Task Programming Model.

    Next Time...

    In my next post I will provide some examples of using the External Calling API to actually do stuff in an externally defined procedure.  However, I am going to spend some time trying to figure out the easiest way for you to get the tools you'll need to be able to develop externally defined functions yourself, so my next post might take a bit of time.

    Darin

    since i can't find the answer in F.A.Q or by searching so i started this thread on my curiosity on BLOGGING ...

    i want to know is there anyway that common members can promote to a degree that allow them become a BLOG Contributor in mapleprimes ?

    what are Qualifications and Requirements ?

     

    thanks a lot for your guides and comments !

    restart:
    with(LinearAlgebra):
    with(ArrayTools):
    with(Statistics):
    randomize():
    with(plots):
    with(combinat):

    n := 100:
    nstock := 7:
    corr := .8:

    R := Matrix(nstock, nstock, proc (i, j) options operator, arrow; `if`(i = j, 1, corr) end proc):
    CD := Matrix(LUDecomposition(evalf(R), 'method' = 'Cholesky'), datatype = float[8]):

    ev := [seq((1/5)*(rand(-3 .. 4))(), i = 1 .. nstock)]:
    st := [seq((rand(1 .. 2))(), i = 1 .. nstock)]:

    just a comment in passing, for the record,

    it appears that the numpoints plot option can be used this way:

    numpoints=round(1e+4)

    but not in the simpler form:

    numpoints=1e+4

    that's because 1e+4 is "1000." with a dot while round(1e+4) is "1000" with no dot.

    There probably are good reasons why 1e+4 is not treated as a natural number, but I haven't investigated them.

    Coon's patch defines a patch from 4 curves forming a chain. 4-examples to Coon's patch are given following worksheet. Boundary curves in the examples are constructed by parametric cubic curves (in Bezier and Hermite form). First example is in 3D and comparable to parametric bicubic surface. Rest of examples are in 2D (abstract surface) which can be described shortly as:

    • 2D Example#1: 4 circular-arc like curves are combined by Coon's patch to form a circular-like region

    I have always preferred the notation  for the derivative of

    My father’s first car when our family moved to North America was a 1970 Buick Skylark sedan, and the color was a majestic deep green. I was seven years old at the time and this was my first experience with a green car. It’s ironic that my life in North America started with a green car and has come full circle with green cars all over the place as far as my job is concerned. But of course, today’s green car is really about highly fuel efficient cars...

    http://www.tiny9.com/u/Math5543

    Fun interactive paradox about the MRB constant.

    Comes with free Mathematica player.

    I notice that a Seach on this site for the word "cipher" fails to show this result as a hit. Why is that?

    It is by no means the first time that a keyword Search on Mapleprimes has failed to return a hit for an old comment/post/answer of mine. I'd really like to know why there are so many such misses.

    Though Maple provides routines to express integers or floating in different
    integer bases the output looks not as nice as I want to have it.

    And for example it does not work as expected using the 'simple ways' (while
    for hex there is no simple way for floats):

      Digits:=16; # for all the following

      x:=frac(evalf(Pi)):
      x:= x* SFloat(1,-14);
                  ...

    This post gives equation of stability region of a Runge-Kutta ODE solver (Dormand-Prince45) as implicit domain equation and ODE-IVP.

    Contents are:

    • Runge-Kutta solver is implemented
    • Solver is applied to test equation and growth factor is obtained
    • Stability regions are defined by separating hyperplane (GF=growth factor=1)
    • Stability regions are plotted by bruteforce approach (implicitplot)
    • Boundary of stability regions...
    First 117 118 119 120 121 122 123 Last Page 119 of 297