Macaulay2 / Macaulay2/M2

minimalBetti sometimes doesn't get cached by Complexes

Open
#4,391 3 comments 0 reactions 0 assignees View on GitHub

Nobody has claimed this yet.

package issue
Dominant language
Macaulay2
Stars
435
Forks
297
Avg merge
4d 20h
Merged PRs (30d)
11

Description

S = ZZ/101[s_0..s_11]
I = ideal{s_8*s_9-s_7*s_10+42*s_9*s_10-42*s_8*s_11, s_6*s_9+42*s_7*s_10-40*s_9*s_10-42*s_6*s_11+40*s_8*s_11,
      s_5*s_9-40*s_7*s_10-28*s_9*s_10-42*s_5*s_11+40*s_6*s_11+28*s_8*s_11, s_4*s_9-28*s_7*s_10+42*s_9*s_10-42*s_4*s_11+40*s_5*s_11+28*s_6*s_11-42*s_8*s_11,
      s_3*s_9+42*s_7*s_10+48*s_9*s_10-42*s_3*s_11+40*s_4*s_11+28*s_5*s_11-42*s_6*s_11-48*s_8*s_11,
      s_2*s_9+48*s_7*s_10+13*s_9*s_10-42*s_2*s_11+40*s_3*s_11+28*s_4*s_11-42*s_5*s_11-48*s_6*s_11-13*s_8*s_11,
      s_1*s_9+13*s_7*s_10+27*s_9*s_10-42*s_1*s_11+40*s_2*s_11+28*s_3*s_11-42*s_4*s_11-48*s_5*s_11-13*s_6*s_11-27*s_8*s_11,
      s_0*s_9+27*s_7*s_10-13*s_9*s_10-42*s_0*s_11+40*s_1*s_11+28*s_2*s_11-42*s_3*s_11-48*s_4*s_11-13*s_5*s_11-27*s_6*s_11+13*s_8*s_11, s_8^2-s_6*s_10, s_7*s_8+7*s_9*s_10-7*s_8*s_11, s_6*s_8-s_5*s_10,
      s_5*s_8-s_4*s_10, s_4*s_8-s_3*s_10, s_3*s_8-s_2*s_10, s_2*s_8-s_1*s_10, s_1*s_8-s_0*s_10, s_7^2-42*s_7*s_9+7*s_9^2-7*s_7*s_11, s_6*s_7+7*s_7*s_10+9*s_9*s_10-7*s_6*s_11-9*s_8*s_11,
      s_5*s_7+9*s_7*s_10-23*s_9*s_10-7*s_5*s_11-9*s_6*s_11+23*s_8*s_11, s_4*s_7-23*s_7*s_10-6*s_9*s_10-7*s_4*s_11-9*s_5*s_11+23*s_6*s_11+6*s_8*s_11,
      s_3*s_7-6*s_7*s_10+9*s_9*s_10-7*s_3*s_11-9*s_4*s_11+23*s_5*s_11+6*s_6*s_11-9*s_8*s_11, s_2*s_7+9*s_7*s_10-33*s_9*s_10-7*s_2*s_11-9*s_3*s_11+23*s_4*s_11+6*s_5*s_11-9*s_6*s_11+33*s_8*s_11,
      s_1*s_7-33*s_7*s_10+10*s_9*s_10-7*s_1*s_11-9*s_2*s_11+23*s_3*s_11+6*s_4*s_11-9*s_5*s_11+33*s_6*s_11-10*s_8*s_11,
      s_0*s_7+10*s_7*s_10+13*s_9*s_10-7*s_0*s_11-9*s_1*s_11+23*s_2*s_11+6*s_3*s_11-9*s_4*s_11+33*s_5*s_11-10*s_6*s_11-13*s_8*s_11, s_6^2-s_4*s_10, s_5*s_6-s_3*s_10, s_4*s_6-s_2*s_10, s_3*s_6-s_1*s_10,
      s_2*s_6-s_0*s_10, s_1*s_6-s_0*s_8, s_5^2-s_2*s_10, s_4*s_5-s_1*s_10, s_3*s_5-s_0*s_10, s_2*s_5-s_0*s_8, s_1*s_5-s_0*s_6, s_4^2-s_0*s_10, s_3*s_4-s_0*s_8, s_2*s_4-s_0*s_6, s_1*s_4-s_0*s_5,
      s_0*s_4-49*s_0*s_5-37*s_0*s_6+33*s_0*s_8-31*s_7*s_9+50*s_9^2-20*s_0*s_10+46*s_1*s_10+30*s_2*s_10+s_3*s_10-36*s_4*s_10-3*s_5*s_10+46*s_6*s_10-45*s_7*s_10+43*s_8*s_10-24*s_9*s_10+13*s_10^2-24*s_4*
      s_11+17*s_5*s_11-4*s_6*s_11+32*s_7*s_11+50*s_8*s_11+49*s_9*s_11+19*s_10*s_11+38*s_11^2, s_3^2-s_0*s_6, s_2*s_3-s_0*s_5,
      s_1*s_3-49*s_0*s_5-37*s_0*s_6+33*s_0*s_8-31*s_7*s_9+50*s_9^2-20*s_0*s_10+46*s_1*s_10+30*s_2*s_10+s_3*s_10-36*s_4*s_10-3*s_5*s_10+46*s_6*s_10-45*s_7*s_10+43*s_8*s_10-24*s_9*s_10+13*s_10^2-24*s_4*
      s_11+17*s_5*s_11-4*s_6*s_11+32*s_7*s_11+50*s_8*s_11+49*s_9*s_11+19*s_10*s_11+38*s_11^2,
      s_0*s_3-14*s_0*s_5+38*s_0*s_6-19*s_0*s_8+46*s_7*s_9-39*s_9^2-25*s_0*s_10-39*s_1*s_10-44*s_2*s_10+13*s_3*s_10-50*s_4*s_10-26*s_6*s_10-7*s_7*s_10-s_8*s_10+46*s_9*s_10+31*s_10^2-24*s_3*s_11-48*s_4*
      s_11+21*s_5*s_11-45*s_6*s_11-35*s_7*s_11+35*s_8*s_11+15*s_9*s_11+22*s_10*s_11+44*s_11^2,
      s_2^2-49*s_0*s_5-37*s_0*s_6+33*s_0*s_8-31*s_7*s_9+50*s_9^2-20*s_0*s_10+46*s_1*s_10+30*s_2*s_10+s_3*s_10-36*s_4*s_10-3*s_5*s_10+46*s_6*s_10-45*s_7*s_10+43*s_8*s_10-24*s_9*s_10+13*s_10^2-24*s_4*s_
      11+17*s_5*s_11-4*s_6*s_11+32*s_7*s_11+50*s_8*s_11+49*s_9*s_11+19*s_10*s_11+38*s_11^2,
      s_1*s_2-14*s_0*s_5+38*s_0*s_6-19*s_0*s_8+46*s_7*s_9-39*s_9^2-25*s_0*s_10-39*s_1*s_10-44*s_2*s_10+13*s_3*s_10-50*s_4*s_10-26*s_6*s_10-7*s_7*s_10-s_8*s_10+46*s_9*s_10+31*s_10^2-24*s_3*s_11-48*s_4*
      s_11+21*s_5*s_11-45*s_6*s_11-35*s_7*s_11+35*s_8*s_11+15*s_9*s_11+22*s_10*s_11+44*s_11^2,
      s_0*s_2-42*s_0*s_5-32*s_0*s_6+33*s_0*s_8+32*s_7*s_9-4*s_9^2-16*s_0*s_10-6*s_1*s_10+29*s_2*s_10-36*s_3*s_10+s_4*s_10+33*s_5*s_10+37*s_6*s_10+22*s_7*s_10+27*s_8*s_10+3*s_9*s_10-20*s_10^2-24*s_2*s_
      11-48*s_3*s_11-12*s_4*s_11-9*s_5*s_11-21*s_6*s_11-45*s_7*s_11-21*s_8*s_11+23*s_9*s_11-37*s_10*s_11+27*s_11^2,
      s_1^2-42*s_0*s_5-32*s_0*s_6+33*s_0*s_8+32*s_7*s_9-4*s_9^2-16*s_0*s_10-6*s_1*s_10+29*s_2*s_10-36*s_3*s_10+s_4*s_10+33*s_5*s_10+37*s_6*s_10+22*s_7*s_10+27*s_8*s_10+3*s_9*s_10-20*s_10^2-24*s_2*s_11
      -48*s_3*s_11-12*s_4*s_11-9*s_5*s_11-21*s_6*s_11-45*s_7*s_11-21*s_8*s_11+23*s_9*s_11-37*s_10*s_11+27*s_11^2,
      s_0*s_1+31*s_0*s_5-6*s_0*s_6-44*s_0*s_8+7*s_7*s_9+24*s_9^2-38*s_0*s_10+42*s_1*s_10+12*s_2*s_10+43*s_3*s_10+36*s_4*s_10+12*s_5*s_10+40*s_6*s_10+32*s_7*s_10-32*s_8*s_10+25*s_9*s_10+41*s_10^2-24*s_
      1*s_11-48*s_2*s_11-12*s_3*s_11-7*s_4*s_11-14*s_5*s_11+13*s_6*s_11+9*s_7*s_11+20*s_8*s_11-36*s_9*s_11-10*s_10*s_11-20*s_11^2,
      s_0^2-2*s_0*s_5-8*s_0*s_6+50*s_0*s_8-25*s_7*s_9+19*s_9^2-45*s_0*s_10+22*s_2*s_10+5*s_3*s_10+17*s_4*s_10+32*s_5*s_10-44*s_6*s_10+6*s_7*s_10+21*s_8*s_10-25*s_9*s_10+s_10^2-24*s_0*s_11-48*s_1*s_11-
      12*s_2*s_11-7*s_3*s_11+23*s_4*s_11-9*s_5*s_11+43*s_6*s_11+29*s_7*s_11+17*s_8*s_11-24*s_9*s_11+17*s_10*s_11+34*s_11^2,
      s_9^3-22*s_7*s_9*s_10-6*s_9^2*s_10-31*s_7*s_10^2-3*s_9*s_10^2+2*s_0*s_5*s_11+18*s_0*s_6*s_11-39*s_0*s_8*s_11+19*s_7*s_9*s_11+50*s_9^2*s_11+36*s_0*s_10*s_11-32*s_1*s_10*s_11+12*s_2*s_10*s_11-9*s_
      3*s_10*s_11+34*s_4*s_10*s_11-25*s_5*s_10*s_11-31*s_6*s_10*s_11+10*s_7*s_10*s_11-41*s_8*s_10*s_11-3*s_9*s_10*s_11-47*s_10^2*s_11-48*s_5*s_11^2-23*s_6*s_11^2-30*s_7*s_11^2+48*s_8*s_11^2-39*s_9*s_
      11^2+9*s_10*s_11^2+18*s_11^3, s_7*s_9^2-21*s_7*s_9*s_10-48*s_9^2*s_10+8*s_7*s_10^2+15*s_9*s_10^2-2*s_0*s_5*s_11-17*s_0*s_6*s_11+49*s_0*s_8*s_11+43*s_7*s_9*s_11+23*s_9^2*s_11+5*s_0*s_10*s_11-10*s
      _1*s_10*s_11+31*s_2*s_10*s_11-43*s_3*s_10*s_11-40*s_4*s_10*s_11+36*s_5*s_10*s_11-21*s_6*s_10*s_11+24*s_8*s_10*s_11-39*s_9*s_10*s_11+20*s_10^2*s_11+48*s_5*s_11^2-s_6*s_11^2+8*s_7*s_11^2-33*s_8*s_
      11^2-s_9*s_11^2+37*s_10*s_11^2-27*s_11^3};
i3 : elapsedTime minimalBetti I;
 -- 3.35017s elapsed

i4 : elapsedTime minimalBetti I;
 -- 3.23379s elapsed

Previously:

i4 : elapsedTime minimalBetti I;
 -- 3.33406s elapsed

i5 : elapsedTime minimalBetti I;
 -- .000518093s elapsed

Also, isn't minimalBetti supposed to give an error if the input is over QQ?

Contributor guide

No contributing guide indexed for this repository

First steps

  1. Read the whole issue, then the project's contributing guide.
  2. Comment on the issue to say you are picking it up — it saves two people doing the same work.
  3. Fork the repository and make your change on a branch.
  4. Open a pull request that references the issue number.

Research direction

Start by reproducing the two consecutive minimalBetti I calls in Macaulay2 using the provided ring and ideal, then compare the elapsed times. Investigate minimalBetti's caching behavior and its response to input over QQ; done means the caching result is consistent and the QQ behavior matches the documented expectation.

Written by the indexing model from the issue text.

Assessment

Domain
tooling
Issue type
Bug
Difficulty
4/5
Estimated time
3-5 days
Activity status
Quiet
Clarity
Mostly clear
Newbie friendliness
45/100

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.