From 4fc2fa90574865a43156c030c9b5e976952045c5 Mon Sep 17 00:00:00 2001
From: Franz Zotter <fzotter@users.sourceforge.net>
Date: Wed, 13 Mar 2013 14:33:53 +0000
Subject: took a too large difference from unsigned int in legendre_a.c.
 corrected m<=n-2 to m+2<=n

svn path=/trunk/externals/iem/iemmatrix/; revision=17061
---
 src/mtx_spherical_harmonics/legendre_a.c | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/mtx_spherical_harmonics/legendre_a.c b/src/mtx_spherical_harmonics/legendre_a.c
index aa42949..602869c 100644
--- a/src/mtx_spherical_harmonics/legendre_a.c
+++ b/src/mtx_spherical_harmonics/legendre_a.c
@@ -45,7 +45,7 @@ static void legendre_second_recurrence (double *costheta, LegendreWorkSpace *wl)
    // P_n^m with m=n
    for (n=1; n<=wl->nmax; n++) {
       for (m=0; m<n; m++) {
-	 if (m<=n-2) {
+	 if (m+2<=n) {
 	    for (l=0,l0=0; l<wl->l; l++,l0+=incr) {
 	       wl->p[l0+n0+m] = (
 		     (2*n-1) * costheta[l] * wl->p[l0+nmo0+m]
-- 
cgit v1.2.1