Bug in HeuristicDecomposition?

[AndrewVSutherland]

I am encountering the following error on many genus 2 curves over Q whose Jacobians are isogenous to products of elliptic curves over Q (I think these used to work):

> R<x>:=PolynomialRing(Rationals());
> C:=HyperellipticCurve(R![96,0,754,0,-118,0,-48]);
> HeuristicDecomposition(C);  

HeuristicDecomposition(
X: Hyperelliptic Curve defined by y^2 = -48*x^6 - 118*x^4 + 754...
)
DecompositionOverClosure(
X: Hyperelliptic Curve defined by y^2 = -48*x^6 - 118*x^4 + 754...
)
DecompositionOverClosure(
P: [-5.94849050897489773423212098305847989293559982752044076020...,
GeoEndoRep: [ [* [1 0] [0 1],  [1 0 0 0] [0 1 0 0] [0 0 1 0] [0 0 0 1] *...
)
ReconstructionFromFactor(
P: [-5.94849050897489773423212098305847989293559982752044076020...,
Q: [0.088028440588597724789372335904497769271430011621550501798...,
mor: [* [0 1],  [ 1  0  0  0] [ 0  1 -2  0] *]
)
ReconstructionFromFactorG1(
P: [-5.94849050897489773423212098305847989293559982752044076020...,
Q: [0.088028440588597724789372335904497769271430011621550501798...,
mor: [* [0 1],  [ 1  0  0  0] [ 0  1 -2  0] *]
)
In file "/home/drew/Dropbox/github/CHIMP/endomorphisms/endomorphisms/magma/heuristic/DecCurves.m", line 22, column 2:
>> E`period_matrix := Q;
    ^
Runtime error in :=: Bad LHS for indexed assign

[AndrewVSutherland]

Possibly the problem is in the curve_reconstruction package and not in endomorphisms.