From c5fec9406d8d18da95a25a1d884cd394e926d156 Mon Sep 17 00:00:00 2001
From: Tomas Stolker <tomas.stolker@phys.ethz.ch>
Date: Sun, 3 May 2020 08:57:41 +0200
Subject: [PATCH] Type checks added to read_filter

---
 species/read/read_filter.py | 23 ++++++++++++++++-------
 1 file changed, 16 insertions(+), 7 deletions(-)

diff --git a/species/read/read_filter.py b/species/read/read_filter.py
index 5a2592c4..cff4e4ad 100644
--- a/species/read/read_filter.py
+++ b/species/read/read_filter.py
@@ -5,10 +5,13 @@
 import os
 import configparser
 
+from typing import Tuple
+
 import h5py
 import numpy as np
 
-from scipy.interpolate import interp1d, InterpolatedUnivariateSpline
+from typeguard import typechecked
+from scipy.interpolate import interp1d, InterpolatedUnivariateSpline, interpolate
 
 from species.data import database
 
@@ -18,8 +21,9 @@ class ReadFilter:
     Class for reading a filter profile from the database.
     """
 
+    @typechecked
     def __init__(self,
-                 filter_name):
+                 filter_name: str) -> None:
         """
         Parameters
         ----------
@@ -42,7 +46,8 @@ def __init__(self,
 
         self.database = config['species']['database']
 
-    def get_filter(self):
+    @typechecked
+    def get_filter(self) -> np.ndarray:
         """
         Function for selecting a filter profile from the database.
 
@@ -69,7 +74,8 @@ def get_filter(self):
 
         return data
 
-    def interpolate_filter(self):
+    @typechecked
+    def interpolate_filter(self) -> interpolate.interp1d:
         """
         Function for linearly interpolating a filter profile.
 
@@ -87,7 +93,8 @@ def interpolate_filter(self):
                         bounds_error=False,
                         fill_value=float('nan'))
 
-    def wavelength_range(self):
+    @typechecked
+    def wavelength_range(self) -> Tuple[float, float]:
         """
         Extract the wavelength range of the filter profile.
 
@@ -103,7 +110,8 @@ def wavelength_range(self):
 
         return np.amin(data[0, ]), np.amax(data[0, ])
 
-    def mean_wavelength(self):
+    @typechecked
+    def mean_wavelength(self) -> float:
         """
         Calculate the weighted mean wavelength of the filter profile.
 
@@ -117,7 +125,8 @@ def mean_wavelength(self):
 
         return np.trapz(data[0, ]*data[1, ], data[0, ]) / np.trapz(data[1, ], data[0, ])
 
-    def filter_fwhm(self):
+    @typechecked
+    def filter_fwhm(self) -> float:
         """
         Calculate the full width at half maximum (FWHM) of the filter profile.