Cleanup

demo/cygnusA.py

@@ -16,7 +16,6 @@ def main():
     nthreads = 8
     rve.set_nthreads(nthreads)
     rve.set_wstacking(False)
-    # FIXME Add script that can open measurement sets and given some statistics
     obs = rve.ms2observations(ms, 'DATA')[0].restrict_to_stokes_i()
 
     rve.set_epsilon(1/10/obs.max_snr())
@@ -29,10 +28,6 @@ def main():
     sky = rve.vla_beam(dom, np.mean(obs.freq)) @ sky
     invcovop = ift.InverseGammaOperator(obs.vis.domain, 1, 1/obs.weight).reciprocal().ducktape('invcov')
 
-    plotter = rve.Plotter()
-    # plotter.add(invcovop)
-    plotter.add(sky, norm=LogNorm())
-
     # imaging_lh = rve.ImagingLikelihood(obs, sky)
     # ham = ift.StandardHamiltonian(imaging_lh, ift.AbsDeltaEnergyController(0.5, 3, 2000))
     # ham = ift.EnergyAdapter(0.1*ift.from_random(ham.domain), ham, want_metric=True)

resolve/__init__.py

@@ -6,6 +6,6 @@ from .likelihood import (CalibrationLikelihood, ImagingCalibrationLikelihood,
                          ImagingLikelihoodVariableCovariance)
 from .ms_import import ms2observations
 from .observation import Observation, tmin_tmax, unique_antennas, unique_times
-from .plotter import Plotter
+# from .plotter import Plotter
 from .points import PointInserter
 from .primary_beam import vla_beam

resolve/plotter.py

@@ -9,6 +9,7 @@ from .util import my_assert_isinstance
 
 class Plotter:
     def __init__(self):
+        raise NotImplementedError
         self._nifty = []
         self._hists = []