Added functions to handle matrix, vector operations
Added some thoughts to the generic handling of IR
Started reworking on the driving function for WFS