some ideas stolen from the cloud

orbcorrsnippet.py

+quads = ring0.get_refpts(at.Quadrupole])                                                                                               
+bpms = ring0.get_refpts(at.Monitor)
+
+rm = LinoptResponseMatrix(ring0, bpms, quads, ['beta', 'mu', 'dispersion'], [[0, 1], [0, 1], 0], fit_tune=True, obsw=[1.0, 1.0, 1000.0])
+rm.compute_fullrm(ring0, use_mp=True)
+target, _ = rm.get_vals(ring0)
+dk, exp = rm.correct(ring_err, target, apply_correction=True)
+
+#Warning!!! for this to work ring and ring_err MUST have the same element sequence because variables and observables are based on their refpts
+
+