From 2f994276ae594f769c1c4b3b63cd42652885dd8f Mon Sep 17 00:00:00 2001
From: Paul Brossier <piem@piem.org>
Date: Fri, 10 Mar 2017 22:10:17 +0100
Subject: [PATCH] src/cvec.h: add cvec_logmag

---
 src/cvec.c | 7 +++++++
 src/cvec.h | 8 ++++++++
 2 files changed, 15 insertions(+)

diff --git a/src/cvec.c b/src/cvec.c
index 2e98cf5f..4ed2a43f 100644
--- a/src/cvec.c
+++ b/src/cvec.c
@@ -139,3 +139,10 @@ void cvec_zeros(cvec_t *s) {
   cvec_norm_zeros(s);
   cvec_phas_zeros(s);
 }
+
+void cvec_logmag(cvec_t *s, smpl_t lambda) {
+  uint_t j;
+  for (j=0; j< s->length; j++) {
+    s->norm[j] = LOG(lambda * s->norm[j] + 1);
+  }
+}
diff --git a/src/cvec.h b/src/cvec.h
index d0a9c2ba..96cc04e5 100644
--- a/src/cvec.h
+++ b/src/cvec.h
@@ -230,6 +230,14 @@ void cvec_phas_ones(cvec_t *s);
 */
 void cvec_zeros(cvec_t *s);
 
+/** take logarithmic magnitude
+
+  \param fftgrain cvec to compress
+  \param lambda value to use for normalisation
+
+*/
+void cvec_logmag(cvec_t *s, smpl_t lambda);
+
 #ifdef __cplusplus
 }
 #endif
-- 
2.11.0