The speed of gravitational waves for a single observation can be measured by the time delay among gravitational-wave detectors with Bayesian inference. Then multiple measurements can be combined to produce a more accurate result. From the near simultaneous detection of gravitational waves and gamma rays originating from GW170817/GRB 170817A, the speed of gravitational wave signal was found to be the same as the the speed of the gamma rays to approximately one part in $10^{15}$. Here we present a different method of measuring the speed of gravitational waves, not based on an associated electromagnetic signal but instead by the measured transit time across a geographically separated network of detectors. While this method is far less precise, it provides an independent measurement of the speed of gravitational waves. For GW170817 a binary neutron star inspiral observed by Advanced LIGO and Advanced Virgo, by fixing sky localization of the source at the electromagnetic counterpart the speed of gravitational waves is constrained to 90% confidence interval (0.97c, 1.02c), where c is the speed of light in a vacuum. By combing seven BBH events and the BNS event from the second observing run of Advanced LIGO and Advanced Virgo, the 90% confidence interval is narrowed down to (0.97c, 1.01c). The accurate measurement of the speed of gravitational waves allows us to test the general theory of relativity. We further interpret these results within the test framework provided by the gravitational Standard-Model Extension (SME). In doing so, we obtain simultaneous constraints on 4 of the 9 nonbirefringent, nondispersive coefficients for Lorentz violation in the gravity sector of the SME and place limits on the anisotropy of the speed of gravity.

The ever-increasing sensitivity of the network of gravitational-wave detectors has resulted in the accelerated rate of detections from compact binary coalescence systems in the third observing run of Advanced LIGO and Advanced Virgo. Not only has the event rate increased, but also the distances to which phenomena can be detected, leading to a rise in the required sky volume coverage to search for counterparts. Additionally, the improvement of the detectors has resulted in the discovery of more compact binary mergers involving neutron stars, revitalizing dedicated follow-up campaigns. While significant effort has been made by the community to optimize single telescope observations, using both synoptic and galaxy-targeting methods, less effort has been paid to coordinated observations in a network. This is becoming crucial, as the advent of gravitational-wave astronomy has garnered interest around the globe, resulting in abundant networks of telescopes available to search for counterparts. In this paper, we extend some of the techniques developed for single telescopes to a telescope network. We describe simple modifications to these algorithms and demonstrate them on existing network examples. These algorithms are implemented in the open-source software \texttt{gwemopt}, used by some follow-up teams, for ease of use by the broader community.