From eaa4eaf0c0c936645b48e5ce0fc877064fd92089 Mon Sep 17 00:00:00 2001 From: Martin Reinecke Date: Wed, 19 Dec 2018 14:19:53 +0100 Subject: [PATCH] safety margin --- libsharp/sharp_ylmgen_c.c | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/libsharp/sharp_ylmgen_c.c b/libsharp/sharp_ylmgen_c.c index 622fc30..f0672b0 100644 --- a/libsharp/sharp_ylmgen_c.c +++ b/libsharp/sharp_ylmgen_c.c @@ -73,9 +73,9 @@ void sharp_Ylmgen_init (sharp_Ylmgen_C *gen, int l_max, int m_max, int spin) gen->mfac[0] = inv_sqrt4pi; for (int m=1; m<=gen->mmax; ++m) gen->mfac[m] = gen->mfac[m-1]*sqrt((2*m+1.)/(2*m)); - gen->root = RALLOC(double,2*gen->lmax+6); - gen->iroot = RALLOC(double,2*gen->lmax+6); - for (int m=0; m<2*gen->lmax+6; ++m) + gen->root = RALLOC(double,2*gen->lmax+7); + gen->iroot = RALLOC(double,2*gen->lmax+7); + for (int m=0; m<2*gen->lmax+7; ++m) { gen->root[m] = sqrt(m); gen->iroot[m] = (m==0) ? 0. : 1./gen->root[m];