Fortraneries/Modules/wavmetrics.f90

86 lines
2.4 KiB
Fortran
Raw Normal View History

2024-02-07 00:50:27 +01:00
module wavmetrics
2024-02-07 02:36:08 +01:00
! new: Thu Jan 4 00:08:04 UTC 2024
2024-02-07 00:50:27 +01:00
use mathstuff2
implicit none
type wavmetric
2024-02-07 02:36:08 +01:00
integer :: nbre ! number of slices/samples
real :: freql, freqr ! zero-crossing estimation
integer :: maxl, maxr ! maximum of abs values
real :: meanl, meanr
2024-02-07 00:50:27 +01:00
end type
type intsample
integer :: left, right
end type
contains
!-------------------------------------------------------------
2024-02-07 02:36:08 +01:00
!-
! main computation routine, still full buggy
!-
subroutine compute_wavmetric(samples, size, metrics)
2024-02-07 00:50:27 +01:00
type(intsample), intent(in) :: samples(:)
2024-02-07 02:36:08 +01:00
integer, intent(in) :: size
2024-02-07 00:50:27 +01:00
type(wavmetric), intent(out) :: metrics
2024-02-07 02:36:08 +01:00
integer :: Lmax, Rmax
integer :: Lval, Rval
integer :: idx
integer :: Lfreq, Rfreq
real :: Lsum, Rsum
Lmax = 0 ; Rmax = 0
Lfreq = 1 ; Rfreq = 1
Lsum = 0.0 ; Rsum = 0.0
do idx=1, size
Lval = samples(idx)%left
Rval = samples(idx)%right
! print *, Rval, Lval
if (abs(Lval) .GT. Lmax) Lmax = abs(Lval)
if (abs(Rval) .GT. Rmax) Rmax = abs(Rval)
if (idx .GT. 1) then
if (diff_sign(samples(idx-1)%left, Lval)) Lfreq = Lfreq + 1
if (diff_sign(samples(idx-1)%right, Lval)) Rfreq = Rfreq + 1
endif
Lsum = Lsum + Lval
Rsum = Rsum + Rval
enddo
metrics%num = size
metrics%maxl = Lmax ; metrics%maxr = Rmax
metrics%freql = 1.0 / (Lfreq / real(size))
metrics%freqr = 1.0 / (Rfreq / real(size))
metrics%meanl = Lsum / real(size)
metrics%meanr = Rsum / real(size)
2024-02-07 00:50:27 +01:00
end subroutine
!-------------------------------------------------------------
subroutine display_wavmetrics(metrics)
type(wavmetric), intent(in) :: metrics
2024-02-07 02:36:08 +01:00
! print '(1X, "metrics are :")'
print '(1X, " | num ", I0)', metrics%num
print '(1X, " | freq ", 2F12.2)', metrics%freql, metrics%freqr
print '(1X, " | mean ", 2F12.2)', metrics%meanl, metrics%meanr
print '(1X, " | maxi ", 2I8)', metrics%maxl, metrics%maxr
2024-02-07 00:50:27 +01:00
end subroutine
!-------------------------------------------------------------
!-------------------------------------------------------------
end module