FDTD_2D_Maxwell/incident.c

/*
    2D FDTD simulator
    Copyright (C) 2019 Emilia Blåsten

    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/>.
*/
/* incident.c: The file for defining the incident wave function.
 *   Input: current time (int)
 *   Output: Amplitude of ez-field on left side of domain (double)*/

#include <math.h>

double incidentWaveOnLeftSide(int time, double cdtds) {
  // Sine-wave on half-space
  return sin(time * cdtds / 2.5);

  // Delta-pulse-line
  //if(time==0)
  //  return 1/cdtds;
  //else
  //  return 0;
}