FhSim  3.1.0
Marine systems simulation
Loading...
Searching...
No Matches
Samplifier.h
1#ifndef SAMPLIFIER_H
2#define SAMPLIFIER_H
3
4#include <sfh/util/CircularBuffer.h>
5
6#include <fhsim/simobject/SimObject.h>
7#include <memory>
8#include <mutex>
9#include <random>
10#include <vector>
11
83class Samplifier : public SimObject
84{
85 public:
86 Samplifier(const std::string& simObjectName, ISimObjectCreator* const creator);
87 ~Samplifier() { };
88
89 virtual void OdeFcn(const double T, const double* const X, double* const XDot) const override;
90 virtual void AcceptedStep(const double T, const double* const X) override;
91 virtual void InitialConditionSetup(const double T, const double* const currentIC, double* const updatedIC, ISimObjectCreator* const creator);
92 virtual void FinalSetup(const double T, const double* const X, ISimObjectCreator* const creator) { m_IsSettingUp = false; };
93 const double* Output(const double T, const double* const X, const int Index);
94
95#ifdef FH_VISUALIZATION
96 // The rendering is not used.
97 virtual void RenderInit(Ogre::Root* const ogreRoot, ISimObjectCreator* const creator) { }
98
99 // The rendering is not used.
100 virtual void RenderUpdate(const double T, const double* const X) { }
101#endif
102
103 protected:
104 std::unique_ptr<ISignalPort*[]> m_Input_i;
105 void SetOutputPortValues(const double T, const double* const X) const;
106 ICommonComputation* m_SetOutputPortValues;
107
108 private:
109 bool isOutlier(double probability);
110 bool fiftyFifty();
111 double getNoise(const int PortNo, const int PortIndex);
112
113 bool m_IsSettingUp;
114 int m_PortWidth, m_seed, m_DummyIndex;
115 std::unique_ptr<int[]> m_PortSize_i;
116
117 std::vector<std::unique_ptr<double[]>> m_noiseBias;
118 std::vector<std::unique_ptr<double[]>> m_noiseStdDev;
119 std::vector<std::unique_ptr<double[]>> m_samplingPeriod;
120 std::vector<std::unique_ptr<double[]>> m_signalSampleTime;
121 std::vector<std::unique_ptr<double[]>> m_firstSampleTime;
122 std::vector<std::unique_ptr<double[]>> m_outlierProbability;
123 std::vector<std::unique_ptr<double[]>> m_lossProbability;
124 std::vector<std::unique_ptr<int[]>> m_lostSamples;
125 std::vector<std::unique_ptr<int[]>> m_avgLostSamples;
126 std::vector<std::unique_ptr<double[]>> m_timeDelay;
127 mutable std::vector<std::unique_ptr<double[]>> m_WriteBuffer;
128 mutable std::vector<std::unique_ptr<double[]>> m_ReadBuffer;
129
130 mutable std::mutex m_SignalResource;
131 std::vector<std::unique_ptr<CircularBuffer<std::pair<double, double>, 1024>[]>> m_measurementBuffer;
132
133 std::mt19937 m_generator;
134 std::uniform_real_distribution<double> m_uniform;
135 std::vector<std::unique_ptr<std::normal_distribution<double>[]>> m_normals;
136};
137
138
139#endif // include guard
Definition Samplifier.h:84