Commit 8b8fa1b1 authored by Giovanni La Mura's avatar Giovanni La Mura
Browse files

Use built-in complex absolute function in cbf()

parent 48a977f2

src/libnptm/sph_subs.cpp

+2 −3
Original line number Diff line number Diff line
@@ -99,8 +99,7 @@ void cbf(int n, dcomplex z, int &nm, dcomplex *csj) { 
   *            point for backward recurrence

   * ==========================================================

   */

  double zz = real(z) * real(z) + imag(z) * imag(z);

  double a0 = sqrt(zz);

  double a0 = cabs(z);

  nm = n;

  if (a0 < 1.0e-60) {

    for (int k = 2; k <= n + 1; k++) {