       Re: CompiledFunction for matrices ??

• To: mathgroup at smc.vnet.net
• Subject: [mg16089] Re: CompiledFunction for matrices ??
• From: dreiss at !SPAMscientificarts.com (David Reiss)
• Date: Thu, 25 Feb 1999 08:24:52 -0500
• References: <7atr4s\$7m9@smc.vnet.net>
• Sender: owner-wri-mathgroup at wolfram.com

```In article <7atr4s\$7m9 at smc.vnet.net>, "Ersek, Ted R"
<ErsekTR at navair.navy.mil> wrote:

> While writing some code to answer a question Peter Klamser sent in I found
> myself wanting to write a CompiledFunction that takes a length (n) vector of
> real numbers.
> The best solution I could find is to explicitly write out (n) variables as
> in:
> foo=Compile[{x1,x2,x3, ... , xn}, expr]
>
> Then I can evaluate something like:
> Apply[foo, {6.2,4.1,2.5,8.6,7.7,9.1,2.2,1.4}]
>
> If I want to write a CompiledFunction that takes an (m) by (n) matrix of
> real numbers it's also a real chore by any method I know of.
>
> The documentation says:
> Compile[{{x1, t1, n1}, ... }, expr] assumes that xi is a
> rank ni array of objects each of a type which matches ti.
>
> I used to think this allowed the sort of thing I am looking for, but  I
> didn't think about it long enough.  It seems
> Compile[{{m,_Real,n}}, expr]   (for n>2) takes a tensor.
>
> __________________
>
> So how can we write a CompiledFunction that takes a large vector or large
> matrix?
>
>
> Thanks,
> Ted Ersek

Hi Ted,

If I understand your question correctly you are confusing the "n"
in Compile[{{m,_Real,n}}, expr] with the size of (number of elements in)
the object in question.  I had this confusion when I first encountered
this new feature in Version 3.  Actually n=1 corresponds to a vector
of arbitrary length, n=2 corresponds to a 2-dimensional
rectangular matrix, etc...

For example, this computes the dot product of a (Real)
vector with itself:

In:=
testCompile1=Compile[{{x, _Real,1 }}, x.x]

Out=
CompiledFunction[{x},x.x,"-CompiledCode-"]

Here is a length 1000 vector:

In:=
testVector=Table[Random[],{j,1000}];

In:=
testCompile1[testVector]

Out=
332.942

Check the result:

In:=
testVector.testVector

Out=
332.942

This computes the product of an p x q matrix with its transpose:

In:=
testCompile2=Compile[{{x, _Real,2 }}, x.Transpose[x]]

Out=
CompiledFunction[{x},x.Transpose[x],"-CompiledCode-"]

Here is a 2 x 1000 matrix:

In:=
testMatrix=Table[Random[],{i,2},{j,1000}];

In:=
testCompile2[testMatrix]

Out=
{{326.89,246.928},{246.928,334.346}}

Check the result:

In:=
testMatrix.Transpose[testMatrix]

Out=
{{326.89,246.928},{246.928,334.346}}

Regards,

David

--

----------------------------------------

Scientific Arts:
Creative Services and Consultation
for the Applied and Pure Sciences

David Reiss
Email: dreiss at !SPAMscientificarts.com

----------------------------------------

Remove the !SPAM to send email

```

• Prev by Date: Re: Strange results from Integrate
• Next by Date: Re: How can I plot for array variable?
• Previous by thread: Re: CompiledFunction for matrices ??
• Next by thread: Bug using Options[]