be more fussy when fixing the Wigner 3j signs (#180)

LSSTDESC · Oct 20, 2023 · d9af8b8 · d9af8b8
1 parent b45317f
commit d9af8b8
Showing 1 changed file with 1 addition and 2 deletions.
diff --git a/src/utils.c b/src/utils.c
@@ -368,8 +368,7 @@ int drc3jj(int il2,int il3,int im2, int im3,int *l1min_out,
 	sumuni=sum1;
 
       cnorm=1./sqrt(sumuni);
-      if(thrcof[nfin-1]<0) sign1=-1;
-      else sign1=1;
+      sign1 = copysign(1., thrcof[nfin-1]);
 
       if(sign1*sign2<=0)
 	cnorm=-cnorm;