MathGroup Archive 2006

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

Search the Archive

Re: Problems solving using Solve

  • To: mathgroup at smc.vnet.net
  • Subject: [mg68828] Re: Problems solving using Solve
  • From: Jean-Marc Gulliet <jeanmarc.gulliet at gmail.com>
  • Date: Sun, 20 Aug 2006 04:43:47 -0400 (EDT)
  • Organization: The Open University, Milton Keynes, UK
  • References: <ec65bn$1f8$1@smc.vnet.net>
  • Sender: owner-wri-mathgroup at wolfram.com

akil wrote:
> I have a problem suing the following code:
> 
> Cone[beta_, v_,rcAB_,rcCD_,ax_,ay_,cx_,cy_,centerOfMassx_,centerOfMassy_] := 
> Module[
> {betaUsed, vtemp, fx, fy, ex, ey, 
> waarde,criticalRC,criticalHeight,xVerplaatsing,ePx2Py2},
>       betaUsed = beta ;
>       vtemp = v;
>       fline[x_] := Tan[betaUsed]x - Tan[betaUsed]vtemp;
>       fx = x /. Solve[fline[x] == cd[x,rcCD,cx,cy], x][[1]];
>       fy = cd[fx,rcCD,cx,cy];
>       criticalRC = (centerOfMassy - fy)/(centerOfMassx - fx);
>       criticalHeight = -criticalRC*fx + fy;
>       ex = (-ay + rcAB*ax + criticalHeight)/(rcAB + (-criticalRC));
>       ey = ab[ex,rcAB,ax,ay];
>       xVerplaatsing = ex - vtemp;
>       ePx2Py2 = Sqrt[xVerplaatsing^2 + ey^2];
>       waarde = xVerplaatsing/ePx2Py2;
>       ArcCos[-waarde]
> ];
> 
> afg = 
> Simplify[D[Cone[beta,v,rcAB,rcCD,ax,ay,cx,cy,centerOfMassx,centerOfMassy],v], 
> TimeConstraint -> 1 ];
> Solve[afg== 0, v];
> 
> eg.
> {ax,ay}= {-44.9376, 36.1678} and {cx,cy}= {17.8329, 122.477} and 
> {centerOfMassx,centerOfMassy}={-18.6338,72.3356} and rcAB=-0.166667 and 
> rcCD=-0.37931
> 
> Doing the Solve[afg== 0, v] takes ages (after an hour I quit it). Is there 
> any way I can make it work faster? Is any code in Cone such that it 
> unnecessarily makes it run faster? So are there any ways I can try to make 
> everything run faster? Using integers it runs, but it takes too long when 
> using real cones (made up of four points A B C D and a centerOfMass).
> 
> Akil 
> 
> 
Hi Akil,

I tried to run your code on my system but the definition of the function 
'cd' is missing.

Anyway, since your are interested in numerical solutions, why don't you 
use NSolve in place of Solve? Solve is looking for an analytic solution 
-- a closed form expression that can be quite time consuming to find, 
assuming that such a solution exists at all -- whereas NSolve uses 
directly numerical algorithms.

HTH,
Jean-Marc


  • Prev by Date: Re: too many special linear matrices
  • Next by Date: Re: too many special linear matrices
  • Previous by thread: Problems solving using Solve
  • Next by thread: RE: Problems solving using Solve