example-2111_1-numerical.m2 fails on M2 1.23

[rfabbri]

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!

[timduff35]

Thanks, I’ll take a look at this within the next week.

[rfabbri]

The errors for M2 1.12 and 1.15 I now understand, these are changes to NumcericalAlgebraicGeometry packages such as reordering arguments and renaming of methods. I think the focus should be to make it work for the most recent M2,
and from what I see the errors in older versions are to help trace back to whatever M2 version you were using. I'll have a look at it.