cache the load status

Add a Struct.hash() method to return a unique hash over attributes.

Modify Noise.run() method to check for changes to the ifo and only run
the update method if the ifo has changed or other kwargs are supplied.

nb.precomp decorator can be used to wrap BudgetItem.calc methods with
functions that only need to be called once per budget run.

    nb.precomp(precomp_foo)
    def calc(self):
        ...

Precomp functions are all executed by the update() method (after attribute
update), and execution state is cached at the Budget level, to prevent
re-calculating the same functions multiple times per update.

The BudgetItem.update() method is no longer expected to be overloaded by
the user.

This provides a speed-up of about 50% for the Aplus model, but much less
for the CE models which are limited by other noise calculations.

API change: Budget/Noise calculations return BudgetTrace object

Closes #59

See merge request !98

This new object holds the PSD and frequency arrays, as well as the
list of all budget sub-noises, accessible via a dictionary interface
and attribute access.

All interfaces and tests are updated to process the new BudgetTrace
object.

A new HDF5 schemata (version 2) is also introduced that stores trace
data in a way that mirrors the BudgetTrace.  The io functions
(save_hdf5/load_hdf5) are updated to write/read this format.  The
load_hdf5 functions continues to support the version 1 SCHEMA.

Closes #59

remove reference to unused BeamRadius parameter in ifo.yamls

Closes #3

See merge request !101

We don't want to advertise these parameters yet, since they will
likely be changed to support different loss functions.

Adding Hong et al. calculations for Coating Brownian Noise

See merge request gwinc/pygwinc!52

Anchal Gupta authored 4 years ago and Jameson Rollins committed 4 years ago

Modified getCoatBrownian() to accept different Bulk and Shear Loss
Angles and calculate brownian noise using Hong et a. PRD 87, 082001
(2013).

New Function description:
This function calculates Coating ThermoOptic noise using
Hong et al . PRD 87, 082001 (2013).
All references to 'the paper', 'Eq' adn 'Sec' are to this paper.

***Important Note***
Inside this function phi is used for denoting the phase shift suffered
by light in one way propagation through a layer. This is in conflict
with present nomenclature everywhere else where it is used as loss
angle.

The layers are assumed to be alernating low-n high-n layers, with
low-n first.
Inputs:
         f = frequency vector in Hz
       ifo = parameter struct from IFOmodel.m
     wBeam = beam radius (at 1 / e^2 power)
      dOpt = the optical thickness, normalized by lambda, of each
             coating layer.
New required arguments:
       mTi = Mirror Transmittance
        Ic = Circulating Laser Power falling on the Mirror (W)
The following new arguments should be made available in the Materials
object and lines 108-110 and 116-118 should be uncommented:
 lossBlown = Coating Bulk Loss Angle of Low Refractive Index layer
 lossSlown = Coating Shear Loss Angle of Low Refractive Index layer
lossBhighn = Coating Bulk Loss Angle of High Refractive Index layer
lossShighn = Coating Shear Loss Angle of High Refractive Index layer
   PETlown = Relevant component of Photoelastic Tensor of High n layer*
  PEThighn = Relevant component of Photoelastic Tensor of Low n layer*
Returns:
SbrZ = Brownian noise spectra for one mirror in m^2 / Hz

*
Default values of PETlown and PEThighn are chosen from sec. A.1. to get
the longitudnal coefficent of photoelasticity as -0.5 as been asserted
by the paper there for Tantala and -0.27 for Silica.
These values also need to be added in Materials object.

API change: load_budget function return separate Budget and IFO

Closes #69

See merge request gwinc/pygwinc!97

This addresses issues about confusion over the IFO being set as a class
attribute of the Budget, and is required to close #69.

newlines '\n' in the string provided to ipshell.ex cause it to choke.

Resolve "Mistake in suspension thermal noise calculation"

Closes #68

See merge request gwinc/pygwinc!96

Christopher Wipf authored 4 years ago and Jameson Rollins committed 4 years ago

The joint order in hForce/vForce must remain consistent between
this code and the thermal noise calculation

Also remove a couple of unused variables

Fixes #68

Remove carrier density noise from pygwinc

Closes #50

See merge request gwinc/pygwinc!90

Christopher Wipf authored 4 years ago and Jameson Rollins committed 4 years ago

Based on the new calculation with Debye screening from
arXiv:2003.05345, there appears to be no plausible scenario where
this noise will affect the design of GW interferometers.  Koji
estimated the noise is below 1e-28/rtHz for Voyager, as did Florian
Bruns (lead author of arXiv:2003.05345).  It's also negligible for
ET, as discussed in the paper.

Fixes #50

Update coating parameters for aLIGO and Aplus

See merge request gwinc/pygwinc!93

Nicholas Demos authored 4 years ago and Jameson Rollins committed 4 years ago

From MIT coating thermal noise experiment results:
* https://dcc.ligo.org/LIGO-G2000880
and LMA 2020 paper:
* https://iopscience.iop.org/article/10.1088/1361-6382/ab77e9

convert more types to float for easier serialization