
We present first results from the third GRavitational lEnsing Accuracy
Testing (GREAT3) challenge, the third in a sequence of challenges for testing
methods of inferring weak gravitational lensing shear distortions from
simulated galaxy images. GREAT3 was divided into experiments to test three
specific questions, and included simulated space and groundbased data with
constant or cosmologicallyvarying shear fields. The simplest (control)
experiment included parametric galaxies with a realistic distribution of
signaltonoise, size, and ellipticity, and a complex point spread function
(PSF). The other experiments tested the additional impact of realistic galaxy
morphology, multiple exposure imaging, and the uncertainty about a
spatiallyvarying PSF; the last two questions will be explored in Paper II. The
24 participating teams competed to estimate lensing shears to within systematic
error tolerances for upcoming StageIV dark energy surveys, making 1525
submissions overall. GREAT3 saw considerable variety and innovation in the
types of methods applied. Several teams now meet or exceed the targets in many
of the tests conducted (to within the statistical errors). We conclude that the
presence of realistic galaxy morphology in simulations changes shear
calibration biases by $\sim 1$ per cent for a wide range of methods. Other
effects such as truncation biases due to finite galaxy postage stamps, and the
impact of galaxy type as measured by the S\'{e}rsic index, are quantified for
the first time. Our results generalize previous studies regarding sensitivities
to galaxy size and signaltonoise, and to PSF properties such as seeing and
defocus. Almost all methods' results support the simple model in which additive
shear biases depend linearly on PSF ellipticity.

We present the results of the first strong lens time delay challenge. The
motivation, experimental design, and entry level challenge are described in a
companion paper. This paper presents the main challenge, TDC1, which consisted
of analyzing thousands of simulated light curves blindly. The observational
properties of the light curves cover the range in quality obtained for current
targeted efforts (e.g.,~COSMOGRAIL) and expected from future synoptic surveys
(e.g.,~LSST), and include simulated systematic errors. \nteamsA\ teams
participated in TDC1, submitting results from \nmethods\ different method
variants. After a describing each method, we compute and analyze basic
statistics measuring accuracy (or bias) $A$, goodness of fit $\chi^2$,
precision $P$, and success rate $f$. For some methods we identify outliers as
an important issue. Other methods show that outliers can be controlled via
visual inspection or conservative quality control. Several methods are
competitive, i.e., give $A<0.03$, $P<0.03$, and $\chi^2<1.5$, with some of
the methods already reaching subpercent accuracy. The fraction of light curves
yielding a time delay measurement is typically in the range $f = $2040\%. It
depends strongly on the quality of the data: COSMOGRAILquality cadence and
light curve lengths yield significantly higher $f$ than does sparser sampling.
Taking the results of TDC1 at face value, we estimate that LSST should provide
around 400 robust timedelay measurements, each with $P<0.03$ and $A<0.01$,
comparable to current lens modeling uncertainties. In terms of observing
strategies, we find that $A$ and $f$ depend mostly on season length, while P
depends mostly on cadence and campaign duration.

Weak gravitational lensing has the potential to place tight constraints on
the equation of the state of dark energy. However, this will only be possible
if shear measurement methods can reach the required level of accuracy. We
present a new method to measure the ellipticity of galaxies used in weak
lensing surveys. The method makes use of direct deconvolution of the data by
the total Point Spread Function (PSF). We adopt a linear algebra formalism that
represents the PSF as a Toeplitz matrix. This allows us to solve the
convolution equation by applying the Hopfield Neural Network iterative scheme.
The ellipticity of galaxies in the deconvolved images are then measured using
second order moments of the autocorrelation function of the images. To our
knowledge, it is the first time full image deconvolution is used to measure
weak lensing shear. We apply our method to the simulated weak lensing data
proposed in the GREAT10 challenge and obtain a quality factor of Q=87. This
result is obtained after applying image denoising to the data, prior to the
deconvolution. The additive and multiplicative biases on the shear power
spectrum are then +0.000009 and +0.0357, respectively.

We present four new seasons of optical monitoring data and six epochs of
Xray photometry for the doublyimaged lensed quasar Q J01584325. The
highamplitude, shortperiod microlensing variability for which this system is
known has historically precluded a time delay measurement by conventional
methods. We attempt to circumvent this limitation by application of a Monte
Carlo microlensing analysis technique, but we are only able to prove that the
delay must have the expected sign (image A leads image B). Despite our failure
to robustly measure the time delay, we successfully model the microlensing at
optical and Xray wavelengths to find a half light radius for soft Xray
emission log(r_{1/2,X,soft}/cm) = 14.3^{+0.4}_{0.5}, an upper limit on the
halflight radius for hard Xray emission log(r_{1/2,X,hard}/cm) <= 14.6 and a
refined estimate of the inclinationcorrected scale radius of the optical
Rband (rest frame 3100 Angstrom) continuum emission region of log(r_s/cm) =
15.6+0.3.