Luigi Tibaldo / gammalib

Commit 681214c59112da9e080e3a98e5a17cdfc9f4e141

Authored by Jürgen Knödlseder
1 parent 5e154566

Initialise frequency and frequency derivatives with zero

Inline Side-by-side
Showing 1 changed file with 6 additions and 6 deletions
src/model/GModelTemporalPhaseCurve.cpp
  View file @681214c
1 1 /***************************************************************************
2 2 * GModelTemporalPhaseCurve.cpp - Temporal phase curve model class *
3 3 * ----------------------------------------------------------------------- *
4   - * copyright (C) 2017-2018 by Juergen Knoedlseder *
  4 + * copyright (C) 2017-2021 by Juergen Knoedlseder *
5 5 * ----------------------------------------------------------------------- *
6 6 * *
7 7 * This program is free software: you can redistribute it and/or modify *
... ... @@ -506,7 +506,7 @@ double GModelTemporalPhaseCurve::phase(const GTime& time) const
506 506 // Set constants
507 507 const double c1 = 0.5;
508 508 const double c2 = 1.0 / 6.0;
509   -
  509 +
510 510 // Compute time since reference time in seconds
511 511 double t = time - mjd();
512 512  
... ... @@ -645,7 +645,7 @@ void GModelTemporalPhaseCurve::init_members(void)
645 645 m_f0.name("F0");
646 646 m_f0.unit("Hz");
647 647 m_f0.scale(1.0);
648   - m_f0.value(1.0);
  648 + m_f0.value(0.0);
649 649 m_f0.range(0.0,1000.0);
650 650 m_f0.free();
651 651 m_f0.gradient(0.0);
... ... @@ -656,7 +656,7 @@ void GModelTemporalPhaseCurve::init_members(void)
656 656 m_f1.name("F1");
657 657 m_f1.unit("s^-2");
658 658 m_f1.scale(1.0);
659   - m_f1.value(0.1);
  659 + m_f1.value(0.0);
660 660 m_f1.range(-1000.0,1000.0);
661 661 m_f1.free();
662 662 m_f1.gradient(0.0);
... ... @@ -667,7 +667,7 @@ void GModelTemporalPhaseCurve::init_members(void)
667 667 m_f2.name("F2");
668 668 m_f2.unit("s^-3");
669 669 m_f2.scale(1.0);
670   - m_f2.value(0.01);
  670 + m_f2.value(0.0);
671 671 m_f2.range(-1000.0,1000.0);
672 672 m_f2.free();
673 673 m_f2.gradient(0.0);
... ... @@ -831,7 +831,7 @@ void GModelTemporalPhaseCurve::load_nodes(const GFilename& filename)
831 831 "with normalizations not exceeding 1.";
832 832 throw GException::invalid_value(G_LOAD_NODES, msg);
833 833 }
834   -
  834 +
835 835 } // endfor: looped over nodes
836 836  
837 837 // If first phase is larger than zero then add last node with phase-1 as
... ...