Receiver.h

// This file is part of openPSTD.                                       //
//                                                                      //
// openPSTD is free software: you can redistribute it and/or modify     //
// it under the terms of the GNU General Public License as published by //
// the Free Software Foundation, either version 3 of the License, or    //
// (at your option) any later version.                                  //
//                                                                      //
// openPSTD 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 General Public License for more details.                         //
//                                                                      //
// You should have received a copy of the GNU General Public License    //
// along with openPSTD.  If not, see <http://www.gnu.org/licenses/>.    //
//                                                                      //

//
// Date: 17-9-15
//
//
// Authors: Omar Richardson
//
//
// Purpose: This class models a receiver; a way to measure the sound
// pressure and velocity on a fixed location within the domain.
//
//
#ifndef OPENPSTD_RECEIVER_H
#define OPENPSTD_RECEIVER_H

#include <vector>
#include <iostream>
#include <memory>
#include "../KernelInterface.h"
#include "Geometry.h"
#include "Domain.h"
#include <string>
#include "kernel_functions.h"
#include "WisdomCache.h"

namespace OpenPSTD {
    namespace Kernel {
        class Receiver {

        public:
            const float x;
            const float y;
            const float z;
            int id;

            std::vector<float> location;

            Point grid_location;

            std::vector<float> grid_offset;

            std::shared_ptr<PSTDSettings> config;

            std::shared_ptr<Domain> container_domain;

            std::vector<float> received_values;

            Receiver(std::vector<float> location, std::shared_ptr<PSTDSettings> config, int id,
                     std::shared_ptr<Domain> container);

            float compute_local_pressure();

        private:
            Eigen::ArrayXcf get_fft_factors(Point size, CalcDirection bt);

            float compute_with_nn();

            float compute_with_si();

            Eigen::ArrayXXf calc_domain_fields(std::shared_ptr<Domain> container, CalcDirection bt);

        };

    }
}

#endif //OPENPSTD_RECEIVER_H