%% demo_HearingBehindWalls
% Demonstrating "split room localization" (mics and srcs separated by a
% soundproof obstacle)
%
%     Copyright (C) Srdan Kitic, Inria Rennes
% 
%     This program is free software: you can redistribute it and/or modify
%     it under the terms of the GNU Affero General Public License as
%     published by the Free Software Foundation, either version 3 of the
%     License, or (at your option) any later version.
% 
%     This program is distributed in the hope that it will be useful,
%     but WITHOUT ANY WARRANTY; without even the implied warranty of
%     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
%     GNU Affero General Public License for more details.
% 
%     You should have received a copy of the GNU Affero General Public License
%     along with this program.  If not, see <http://www.gnu.org/licenses/>.
% 
% Contact: nancy.bertin@irisa.fr or remi.gribonval@inria.fr
% PANAMA Research Group,  IRISA  -  (Inria & CNRS)
% Campus de Beaulieu, 35042 Rennes cedex, France
%
%%
clc
clear all

deleteOutputFile = true;

% Setup waveeq/tools search path (but NOT Optimization!)
setupwave

%% 

% Solves min_x ||Ax||_{2,1} + epsilon||Mx-y||^2 + sigma||A_{c}x||^2
% (A_{c} corresponds to pressumed zeros in the spacetime domain)
%
% Set fileName='2_2D.mat' to simulate 2D setting, or fileName='2_3D.mat' to simulate 3D setting
% Simulating for k=2 sources and m=10 microphones [INVERSE CRIME!]
%
% Finer parameterization possible by calling "ProblemSpecifications.m" script:
% it will generate an appropriate .mat file containing problem parameters
%

fileName = 'demo_HBW_2_2D.mat';
nExperiments = 1;
epsilon = 100;
sigma = 100;

[wf,xHat,zHat,outputFile] = SplitRoom(fileName, nExperiments, epsilon, sigma);

%%

if deleteOutputFile
    delete(outputFile)
    disp('Output .mat file deleted.')
else
    disp(['Output .mat file recorded (' outputFile '). Delete to restart the demo.'])
end
    
% Restores previous search path
setupwave(0)