Hi,
I tried to run the software per supplementary.pdf
, but I get this on the most recent M2 1.23:
i4 : load("example-2111_1-numerical.m2")
--loading configuration for package "NumericalAlgebraicGeometry" from file /Users/rfabbri/Library/Application Support/Macaulay2/init-NumericalAlgebraicGeometry.m2
--loading configuration for package "Bertini" from file /Users/rfabbri/Library/Application Support/Macaulay2/init-Bertini.m2
-- warning: symbol "old'random'Type" in User#"private dictionary" is shadowed by a symbol in MonodromySolver#"private dictionary"
-- use the synonym old'random'Type$0
warning: 'rowSelector' redefined
/usr/local/share/Macaulay2/NumericalAlgebraicGeometry/systems.m2:21:0-21:11: here is the first use of 'rowSelector'
: debug with expression debug 4808 or with command line option --debug 4808
warning: 'pCompose' redefined
/usr/local/share/Macaulay2/MonodromySolver/galois-group.m2:2:0-2:8: here is the first use of 'pCompose'
: debug with expression debug 8688 or with command line option --debug 8688
warning: 'writePermutations' redefined
/usr/local/share/Macaulay2/MonodromySolver/galois-group.m2:100:8-100:25: here is the first use of 'writePermutations'
: debug with expression debug 1619 or with command line option --debug 1619
currentString:1:5:(3):[12]: error: no method for assignment to binary operator _ applied to objects:
-- z (of class Rz)
-- _ 2 (of class ZZ)
currentString:1:5:(3):[12]: --entering debugger (type help to see debugger commands)
currentString:1:0-1:18: --source code:
z_2 = inputGate z_2
The problems seems to be when running the following code similar to the initial part of the file:
restart
setRandomSeed 0
m=3 -- number of cameras
D = (4,0,{{0,1},{0,2},{1,2}}) -- encoding for CLEVELAND, a problem of degree 216
needsPackage "NumericalAlgebraicGeometry"
needs "numerical-problem-builder.m2" -- <-- This is what gives the problem.
The remaining example files of section 2.1 run without problem:
Macaulay2, version 1.23
--loading configuration for package "FourTiTwo" from file /Users/rfabbri/Library/Application Support/Macaulay2/init-FourTiTwo.m2
with packages: ConwayPolynomials, Elimination, IntegralClosure, InverseSystems, Isomorphism, LLLBases, MinimalPrimes, OnlineLookup, PrimaryDecomposition, ReesAlgebra, Saturation,
TangentCone, Truncations, Varieties
i1 : load("example-2111_1-unrolled.m2")
i2 : load("example-2111_1.m2")
--building: 3 points, 1 pins, 1 free lines.
-- dependent points:
MutableHashTable{2 => {0, 1}}
-- commonPoint: #equations = 15
-- commonPoint: #equations = 15
-- commonPoint: #equations = 15
-- seeLine: #equations = 4
-- seeLine: #equations = 4
-- seeLine: #equations = 4
i3 : load("example-2111_1-jacobian.m2")
--building: 3 points, 1 pins, 1 free lines.
-- dependent points:
MutableHashTable{2 => {0, 1}}
Tested on MacOS X 12.6.6
Trying M2 1.12 gives:
Macaulay2, version 1.12
with packages: ConwayPolynomials, Elimination, IntegralClosure, InverseSystems, LLLBases, PrimaryDecomposition, ReesAlgebra, TangentCone
i1 : load("example-2111_1-numerical.m2")
--warning: symbol "Jacobian" in ReesAlgebra.Dictionary is shadowed by a symbol in NAGtypes.Dictionary
-- use the synonym ReesAlgebra$Jacobian
--loading configuration for package "NumericalAlgebraicGeometry" from file /Users/rfabbri/Library/Application Support/Macaulay2/init-NumericalAlgebraicGeometry.m2
--loading configuration for package "PHCpack" from file /Users/rfabbri/Library/Application Support/Macaulay2/init-PHCpack.m2
--loading configuration for package "Bertini" from file /Users/rfabbri/Library/Application Support/Macaulay2/init-Bertini.m2
--warning: symbol "sampleInitFile" in User#"private dictionary" is shadowed by a symbol in Core#"private dictionary"
-- use the synonym sampleInitFile$0
-- 0.019499 seconds elapsed
number of polynomials is 60
{0, 1, 2}
common.m2:327:13:(3):[11]: error: no method found for applying evaluate to:
argument 1 : SLProgram ( (of class RawSLProgram)
inputs: 17
#mNodes: 0
outputs: 12
.............
argument 2 : | -.726792+.686857i .080013+.99679. (of class MutableMatrix)
argument 3 : 0 (of class MutableMatrix)
common.m2:327:13:(3):[11]: --entering debugger (type help to see debugger commands)
common.m2:327:13-327:66: --source code:
evaluate(cam, mutableMatrix{sampleCameraParameters}, M);
Trying M2 1.15 gives:
Macaulay2, version 1.15
--loading configuration for package "FourTiTwo" from file /Users/rfabbri/Library/Application Support/Macaulay2/init-FourTiTwo.m2
--loading configuration for package "Topcom" from file /Users/rfabbri/Library/Application Support/Macaulay2/init-Topcom.m2
with packages: ConwayPolynomials, Elimination, IntegralClosure, InverseSystems, LLLBases, PrimaryDecomposition, ReesAlgebra, TangentCone, Truncations
+
i1 : load("example-2111_1-numerical.m2")
--loading configuration for package "NumericalAlgebraicGeometry" from file /Users/rfabbri/Library/Application Support/Macaulay2/init-NumericalAlgebraicGeometry.m2
--loading configuration for package "PHCpack" from file /Users/rfabbri/Library/Application Support/Macaulay2/init-PHCpack.m2
--loading configuration for package "Bertini" from file /Users/rfabbri/Library/Application Support/Macaulay2/init-Bertini.m2
--warning: symbol "sampleInitFile" in User#"private dictionary" is shadowed by a symbol in Core#"private dictionary"
-- use the synonym sampleInitFile$0
-- 0.01712 seconds elapsed
number of polynomials is 60
{0, 1, 2}
-- 0.00331732 seconds elapsed
{0, 1, 2}
example-2111_1-numerical.m2:13:22:(3):[6]: error: no method found for applying evaluate to:
argument 1 : {{(((((l * 1) + (l * 0) + (. (of class GateMatrix)
1,0,1 1,0,2 .
argument 2 : | -.726792+.686857ii | (of class Matrix)
| .0800129+.996794ii |
| .742918+.669382ii |
| -.939129+.343564ii |
.......................
example-2111_1-numerical.m2:13:22:(3):[6]: --entering debugger (type help to see debugger commands)
example-2111_1-numerical.m2:13:22-13:36: --source code:
max(cStarts/(c->norm evaluate(F,c||yStart)))
Thanks!