PulseWave.hh

 
//
// Copyright (C) 2011-2022 Yokokura, Yuki
// This program is free software;
// you can redistribute it and/or modify it under the terms of the FreeBSD License.
// For details, see the License.txt file.
 
#ifndef PULSEWAVE
#define PULSEWAVE
 
#include "Matrix.hh"
 
namespace ARCS {    // ARCS名前空間
double PulseWave(const double freq, const double phase, const double time); 
double PulseWave(const double freq, const double phase, const double time, const double starttime); 
 
template <size_t M>
void PulseWave(const double freq, const double phase, const double time, Matrix<1,M>& output){
    double r = sin(2.0*M_PI*freq*time + phase);
    if(0 < r){
        output = Matrix<1,M>::ones();
    }else{
        output = Matrix<1,M>::zeros();
    }
}
 
template <size_t M>
Matrix<1,M> PulseWaveVec(const double freq, const double phase, const double time){
    Matrix<1,M> ret;
    PulseWave(freq, phase, time, ret);
    return ret;
}
 
}
 
#endif