src/cvec.h: add cvec_logmag

diff --git a/src/cvec.c b/src/cvec.c
index 2e98cf5..4ed2a43 100644 (file)
--- 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 d0a9c2b..96cc04e 100644 (file)
--- 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