simon987/wavelib
1
0
Fork 0
mirror of https://github.com/simon987/wavelib.git synced 2025-04-19 18:16:44 +00:00
Code Issues Projects Releases Wiki Activity

COMMIT : ICWT available only for wt->pow = 2

Browse Source
This commit is contained in:
Rafat Hussain 2016-10-13 14:39:00 +05:30
parent 35c426bb12
commit eeae494901
1 changed files with 2 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
src/wavelib.c
View File

@ -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->scale[i] = s0*pow((double) power, (double)(i)*dj);
} }
wt->sflag = 1; wt->sflag = 1;
wt->pow = power;
} }
else if (!strcmp(wt->type, "lin") || !strcmp(wt->type, "linear")) { 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); cdel = cdelta(wt->mother, wt->m, psi);
//printf("\n PSI %g CDEL %g param %g mother %d \n", psi, cdel,wt->m,wt->mother); //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); icwavelet(wt->params, N, wt->params+nj2, wt->J, wt->dt, wt->dj, cdel, psi, cwtop);
} else { } else {
printf("Inverse CWT is only available for power of 2.0 scales \n"); printf("Inverse CWT is only available for power of 2.0 scales \n");