From eeae494901d957766166350b48d51a09bc4f3131 Mon Sep 17 00:00:00 2001
From: Rafat Hussain <rafat.hsn@gmail.com>
Date: Thu, 13 Oct 2016 14:39:00 +0530
Subject: [PATCH] COMMIT : ICWT available only for wt->pow = 2

---
 src/wavelib.c | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/src/wavelib.c b/src/wavelib.c
index d33671c..a2157d9 100644
--- a/src/wavelib.c
+++ b/src/wavelib.c
@@ -1334,6 +1334,7 @@ void setCWTScales(cwt_object wt, double s0, double dj,char *type,int power) {
 			wt->scale[i] = s0*pow((double) power, (double)(i)*dj);
 		}
 		wt->sflag = 1;
+		wt->pow = power;
 		
 	}
 	else if (!strcmp(wt->type, "lin") || !strcmp(wt->type, "linear")) {
@@ -1418,7 +1419,7 @@ void icwt(cwt_object wt, double *cwtop) {
 	cdel = cdelta(wt->mother, wt->m, psi);
 
 	//printf("\n PSI %g CDEL %g param %g mother %d \n", psi, cdel,wt->m,wt->mother);
-	if (!strcmp(wt->type, "pow") || !strcmp(wt->type, "power")) {
+	if ((!strcmp(wt->type, "pow") || !strcmp(wt->type, "power")) && wt->pow == 2) {
 		icwavelet(wt->params, N, wt->params+nj2, wt->J, wt->dt, wt->dj, cdel, psi, cwtop);
 	} else {
 		printf("Inverse CWT is only available for power of 2.0 scales \n");