Services & Resources / Wolfram Forums
-----
 /
MathGroup Archive
2001
*January
*February
*March
*April
*May
*June
*July
*August
*September
*October
*November
*December
*Archive Index
*Ask about this page
*Print this page
*Give us feedback
*Sign up for the Wolfram Insider

MathGroup Archive 2001

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Pattern matching "on the fly"

  • To: mathgroup at smc.vnet.net
  • Subject: [mg30112] Re: Pattern matching "on the fly"
  • From: "Allan Hayes" <hay at haystack.demon.co.uk>
  • Date: Sat, 28 Jul 2001 22:08:50 -0400 (EDT)
  • References: <9jtk77$87r$1@smc.vnet.net>
  • Sender: owner-wri-mathgroup at wolfram.com

Tony,
Not a general solution, but for your example we get


(1+x)^1000000 +O[x]^3//Normal//Timing

{0.27*Second, 1 + 1000000*x + 499999500000*x^2}
--
Allan
---------------------
Allan Hayes
Mathematica Training and Consulting
Leicester UK
www.haystack.demon.co.uk
hay at haystack.demon.co.uk
Voice: +44 (0)116 271 4198
Fax: +44 (0)870 164 0565

"Tony MacKenzie" <mackenzi at usq.edu.au> wrote in message
news:9jtk77$87r$1 at smc.vnet.net...
> Hello from down under,
>
> I am trying to use Mathematica's pattern matching abilities to achieve
> something I can do using the computer algebra software REDUCE.  I want to
do
> this because I need to take advantage of the powerful numerical algorithms
> of Mathematica.
>
> I want to set all powers of some variable greater than a certain value to
> zero, but I want mathematica to apply this pattern matching "on the fly"
> (while an expression is being evaluated, not after the expression has been
> evaluated).
>
> I'll give a quick example. I want to set all powers of x greater than 2 to
> zero i.e. x^3=>0, x^4=>0 and so on. In REDUCE I can use the following
> statement
>
> let x^3=>0;   (* This replaces x^3 with 0, x^4 with 0 and so on*)
> Then if I evaluate (1+x)^1000000, I quickly (a few seconds) find this
> expression gives 499999500000*x^2  + 1000000*x + 1.  The pattern matching
is
> done as the expression is being expanded.
>
> Now in Mathematica I have tried  x^n_ ^:=0 /; n> 2; This works, but it
only
> appears to be applied after an expression has been expanded.  For example
if
> I try
>
> x^n_ ^:=0 /; n> 2;
> Expand[(1+x)^1000000];
>
> in Mathematica, the evaluation is very slow (which I think is because the
> pattern matching is applied after the expansion and not on the fly).
>
> Any help would be greatly appreciated as I am very new to Mathematica.
>
> Tony MacKenzie
> University of Southern Queensland, Australia
>
>
>
>
>




  • Prev by Date: A Penny for Thought
  • Next by Date: Re: Programming Problem
  • Previous by thread: Re: Pattern matching "on the fly"
  • Next by thread: Re: Pattern matching "on the fly"