/*
**  CUBIC.C - Solve a cubic polynomial
**  public domain by Ross Cottrell
*/

#include <math.h>
#include <stdlib.h>

void SolveCubic(double  a,
                double  b,
                double  c,
                double  d,
                int    *solutions,
                double *x)
{
      long double    a1 = b/a, a2 = c/a, a3 = d/a;
      long double    Q = (a1*a1 - 3.0*a2)/9.0;
      long double R = (2.0*a1*a1*a1 - 9.0*a1*a2 + 27.0*a3)/54.0;
      double    R2_Q3 = R*R - Q*Q*Q;

      double    theta;

      if (R2_Q3 <= 0)
      {
            *solutions = 3;
            theta = acos(R/sqrt(Q*Q*Q));
            x[0] = -2.0*sqrt(Q)*cos(theta/3.0) - a1/3.0;
            x[1] = -2.0*sqrt(Q)*cos((theta+2.0*M_PI)/3.0) - a1/3.0;
            x[2] = -2.0*sqrt(Q)*cos((theta+4.0*M_PI)/3.0) - a1/3.0;
      }
      else
      {
            *solutions = 1;
            x[0] = pow(sqrt(R2_Q3)+fabs(R), 1/3.0);
            x[0] += Q/x[0];
            x[0] *= (R < 0.0) ? 1 : -1;
            x[0] -= a1/3.0;
      }
}

#ifdef TEST

int main(void)
{
      double  a1 = 1.0, b1 = -10.5, c1 = 32.0, d1 = -30.0;
      double  a2 = 1.0, b2 = -4.5, c2 = 17.0, d2 = -30.0;
      double  x[3];
      int     solutions;

      SolveCubic(a1, b1, c1, d1, &solutions, x);

      /* should get 3 solutions: 2, 6 & 2.5   */

      SolveCubic(a2, b2, c2, d2, &solutions, x);

      /* should get 1 solution: 2.5           */

      return 0;
}

#endif /* TEST */
ef='/cgit.cgi/linux/net-next.git/?h=nds-private-remove'>summary</a><a href='/cgit.cgi/linux/net-next.git/refs/?h=nds-private-remove&amp;id=1254856c07e0e5707d3291f625db8ed69a866be2'>refs</a><a class='active' href='/cgit.cgi/linux/net-next.git/log/sound/soc/txx9/txx9aclc.c?h=nds-private-remove'>log</a><a href='/cgit.cgi/linux/net-next.git/tree/sound/soc/txx9/txx9aclc.c?h=nds-private-remove&amp;id=1254856c07e0e5707d3291f625db8ed69a866be2'>tree</a><a href='/cgit.cgi/linux/net-next.git/commit/sound/soc/txx9/txx9aclc.c?h=nds-private-remove&amp;id=1254856c07e0e5707d3291f625db8ed69a866be2'>commit</a><a href='/cgit.cgi/linux/net-next.git/diff/sound/soc/txx9/txx9aclc.c?h=nds-private-remove&amp;id=1254856c07e0e5707d3291f625db8ed69a866be2'>diff</a></td><td class='form'><form class='right' method='get' action='/cgit.cgi/linux/net-next.git/log/sound/soc/txx9/txx9aclc.c'>
<input type='hidden' name='h' value='nds-private-remove'/><input type='hidden' name='id' value='1254856c07e0e5707d3291f625db8ed69a866be2'/><select name='qt'>
<option value='grep'>log msg</option>
<option value='author'>author</option>
<option value='committer'>committer</option>
<option value='range'>range</option>
</select>
<input class='txt' type='search' size='10' name='q' value=''/>
<input type='submit' value='search'/>
</form>
</td></tr></table>
<div class='path'>path: <a href='/cgit.cgi/linux/net-next.git/log/?h=nds-private-remove&amp;id=1254856c07e0e5707d3291f625db8ed69a866be2'>root</a>/<a href='/cgit.cgi/linux/net-next.git/log/sound?h=nds-private-remove&amp;id=1254856c07e0e5707d3291f625db8ed69a866be2'>sound</a>/<a href='/cgit.cgi/linux/net-next.git/log/sound/soc?h=nds-private-remove&amp;id=1254856c07e0e5707d3291f625db8ed69a866be2'>soc</a>/<a href='/cgit.cgi/linux/net-next.git/log/sound/soc/txx9?h=nds-private-remove&amp;id=1254856c07e0e5707d3291f625db8ed69a866be2'>txx9</a>/<a href='/cgit.cgi/linux/net-next.git/log/sound/soc/txx9/txx9aclc.c?h=nds-private-remove&amp;id=1254856c07e0e5707d3291f625db8ed69a866be2'>txx9aclc.c</a></div><div class='content'><table class='list nowrap'><tr class='nohover'><th class='left'>Age</th><th class='left'>Commit message (<a href='/cgit.cgi/linux/net-next.git/log/sound/soc/txx9/txx9aclc.c?h=nds-private-remove&amp;id=1254856c07e0e5707d3291f625db8ed69a866be2&amp;showmsg=1'>Expand</a>)</th><th class='left'>Author</th><th class='left'>Files</th><th class='left'>Lines</th></tr>