Solved

Can you simplify this mathematical algorithm

Posted on 2014-11-19
6
187 Views
Last Modified: 2014-11-20
Is it possible to refactor this mathematical algorithm to minimize the use of trigonometric functions?

input:
  z      :=  0 to ½PI
  Ma, Mb :=  0 to ∞
  Pa, Pb :=  0 to 2PI

algorithm:
  r = Ma * sin(Pa)
  i = Mb * sin(Pb)
  Sr = r + i
  Dr = r - i
  r = Ma * cos(Pa)
  i = Mb * cos(Pb)
  Si = r + i
  Di = r - i
  r = -cos(z)
  i = -sin(z)
  Tr = r * Si + i * Dr
  Ti = i * Si - r * Dr

output:
  Ra :=  Sr + Tr
  Ia :=  Ti + Di
  Rb :=  Sr - Tr
  RI :=  Ti - Di

Open in new window

0
Comment
Question by:JasonMewes
  • 3
  • 2
6 Comments
 
LVL 84

Accepted Solution

by:
ozo earned 500 total points
ID: 40451983
What do Ra,Ia,Rb,RI represent?
Unless you can establish some relationship between the inputs, there may not be much more you can do than to reduce half of the trigonometric functions to sqrt functions.
Or if you can transform your inputs so that instead of working with Ma,Mb,Pa,Pb,
you're working with inputs that correspond to  Ma*sin(Pa),Ma*cos(Pa),Mb*sin(Pb),Mb*cos(Pb),
But that may depend on refactoring other parts of the algorithm that use this part.
0
 

Author Comment

by:JasonMewes
ID: 40452008
Errors (for reference, no effect on question/solution):
Line 3 - bounds are actually -∞ to +∞
Line 24 - the name RI is actually Ib

Ra, Ia, Rb, Ib represent the four needed outputs of the algorithm.
There is no relationship between the inputs, only the bounds stated.

Other parts cannot be refactored, as inputs are external and outputs at this stage must be as specified.

SQRT is costlier than SIN and COS in this case.

Essentially you are saying that this is already as optimized as it can get?
Not possible to use identities or relationships (such as sin(A + B) = sin(A)*cos(B) + cos(A)*sin(B)) to reduce?
Given your (ozo) track record I trust your opinion implicitly!
0
 

Author Comment

by:JasonMewes
ID: 40452018
This probably does not help:

Ma, Pa = magnitude, phase pair A
Mb, Pb = magnitude, phase pair B

Ra, Ia = complex number out A
Rb, Ib = complex number out B
0
Netscaler Common Configuration How To guides

If you use NetScaler you will want to see these guides. The NetScaler How To Guides show administrators how to get NetScaler up and configured by providing instructions for common scenarios and some not so common ones.

 
LVL 84

Expert Comment

by:ozo
ID: 40452126
I'm not sure if this comes out simpler, but
Tr = -(Ma * cos(z-Pa) + Mb * cos(z+Pb))
Ti = -(Ma * sin(z-Pa) + Mb * sin(z+Pb))
0
 

Author Comment

by:JasonMewes
ID: 40452169
Since Sr and Di (requiring four sin/cos to calculate) are still required to calculate output - it seems calculating Tr and Ti this way would add additional processing.
0
 
LVL 16

Expert Comment

by:dhsindy
ID: 40452943
This page shows a history of how trig functions have been computed in the past and how complicated it becomes.  If you read this you will probably be convinced to just use the functions defined in program.  That is probably the simplest way.

http://www.clarku.edu/~djoyce/trig/compute.html
0

Featured Post

Is Your AD Toolbox Looking More Like a Toybox?

Managing Active Directory can get complicated.  Often, the native tools for managing AD are just not up to the task.  The largest Active Directory installations in the world have relied on one tool to manage their day-to-day administration tasks: Hyena. Start your trial today.

Question has a verified solution.

If you are experiencing a similar issue, please ask a related question

Suggested Solutions

Title # Comments Views Activity
Question of probability 6 101
Graph 5 75
algorithm 15 110
Properties of numbers 14 48
Article by: Nadia
Linear search (searching each index in an array one by one) works almost everywhere but it is not optimal in many cases. Let's assume, we have a book which has 42949672960 pages. We also have a table of contents. Now we want to read the content on p…
Have you ever thought of installing a power system that generates solar electricity to power your house? Some may say yes, while others may tell me no. But have you noticed that people around you are now considering installing such systems in their …
This is a video describing the growing solar energy use in Utah. This is a topic that greatly interests me and so I decided to produce a video about it.
Finds all prime numbers in a range requested and places them in a public primes() array. I've demostrated a template size of 30 (2 * 3 * 5) but larger templates can be built such 210  (2 * 3 * 5 * 7) or 2310  (2 * 3 * 5 * 7 * 11). The larger templa…

770 members asked questions and received personalized solutions in the past 7 days.

Join the community of 500,000 technology professionals and ask your questions.

Join & Ask a Question