function WCALCMRA, sig, LD, QMF
;+
; NAME:
;	WCALCMRA
;
; PURPOSE:
;	This procedure calculates a multi-resolution without the plot. 
;	It is used for testing purposes only as it only calculcates the
;	first dyad level.
;
; CATEGORY:
;	Wavelets.
;
; CALLING SEQUENCE:
;	mra = WCALCMRA(sig, LD, QMF)
;
; INPUTS:
;	sig:	1-d signal array
;	LD: 	lowest level resolution	
;	QMF: Quadrature Mirror Filter
;
; OUTPUTS:
;	mra: mra array.
;
; EXAMPLE:
;		Truncate sig array appropriately and get QMF
;		IDL> WINTWAVE, sig, WaveType, ParVal, QMF, LD
;		Call multi-resolution
;		IDL> mra = WCALCMRA(sig, LD, QMF)
;
;
; MODIFICATION HISTORY:
; 	Written by: Amara Graps February 1996
;-


;Calculate Periodized Orthononormal Forward Wavelet Transform
wc = WFWTPO(sig, LD, QMF )

;Initialize mra
mra = wc

wavecoef = WSHASROW(wc)

;Handle autoscaling
scal = 1. / MAX(ABS(wavecoef)) 
!c = 0

n = WDYADLNG(wavecoef,J)

jj=(J-1)
w = FLTARR(n)
w(WDYAD(jj)-1) = wavecoef(DYAD(jj)-1)
x = WIWTPO(w,jj,qmf)
mra(0: 2^(jj+1)-1) = scal*x

RETURN, mra

END