CoCoA » CoCoALib

24 Oct 2016

17:51 Feature #951 (In Progress): New function: IsSqFree

A function which tests whether a value is square-free is probably going to be quite similar to a function which deter... John Abbott

17:29 Feature #951: New function: IsSqFree

Currently there is a function called @IsRadical(PP)@ for a PPMonoidElem.
Should it be called @IsSquareFree@?
"A n... Anna Maria Bigatti

17:19 Feature #951 (Closed): New function: IsSqFree

Write function IsSquareFree (or IsSqFr, or IsSqFree) for a RingElem Anna Maria Bigatti

17:14 Feature #796: CoCoALib function for radical (or SqFree) of a polynomial

I like the idea of a fn @IsSqFree@; it certainly could be faster than testing equality with the radical (and also cle... John Abbott

16:54 Feature #796: CoCoALib function for radical (or SqFree) of a polynomial

After spending some time thinking and writing, I realized that we would need a
@IsSqFree(f)@ which could be consider... Anna Maria Bigatti

15:32 Feature #838: Differential algebra

After some more reading, it seems that the usual structure is @QQ(x_1, x_2, ..., x_n)[u1, u2,...]@ where the expressi... John Abbott

21 Oct 2016

14:24 Feature #947: IsRadical for ideal?

I think there could be problems with removing indets by giving them values: if there is an embedded 0-dim component t... John Abbott

14:14 Design #950 (In Progress): factor and SmoothFactor for integers --> FactorINT, FactorINT_TrialDiv, FactorINT_PollardRho

Does that mean that you are in favour of offering the user the choice between a fast (and probably correct) answer, o... John Abbott

07:46 Design #950: factor and SmoothFactor for integers --> FactorINT, FactorINT_TrialDiv, FactorINT_PollardRho

@ProbSmoothFactor@?
or @SmoothFactor(n)@ and @SmoothFactor(n, "NoProb")@? Anna Maria Bigatti

20 Oct 2016

20:10 Design #950 (Closed): factor and SmoothFactor for integers --> FactorINT, FactorINT_TrialDiv, FactorINT_PollardRho

I can make @factor@ and @SmoothFactor@ for integers faster in some cases by calling @IsProbPrime@.
The risk is that ... John Abbott

20:08 Feature #796 (Resolved): CoCoALib function for radical (or SqFree) of a polynomial

Checked in the code. Also two new tests @test-NumTheory6@ and @test-SparsePolyRing3@.
No documentation yet.
Name... John Abbott

16:10 Feature #796: CoCoALib function for radical (or SqFree) of a polynomial

John Abbott wrote:
> Currently I have used *@radical@* in CoCoALib, and *@rad@* in CoCoA-5 (to avoid a clash with th... Anna Maria Bigatti

15:44 Feature #796: CoCoALib function for radical (or SqFree) of a polynomial

Before checking in which name should I use?
Currently I have used *@radical@* in CoCoALib, and *@rad@* in CoCoA-5 ... John Abbott

13:30 Feature #796 (In Progress): CoCoALib function for radical (or SqFree) of a polynomial

I have an implementation in CoCoALib. It took so long because of a "mysterious bug" in @ContentFreeFactor@.
My cu... John Abbott

16:32 Feature #838 (In Progress): Differential algebra

The ctor args for a @DifferentialAlgebra@ are the following:
* *@n@* number of "independent variables"
* *@k@* numb... John Abbott

19 Oct 2016

10:26 Feature #440: Port RealRoots to C++

Some of Erika's people would like this to be in CoCoALib sooner rather than later!
John Abbott

18 Oct 2016

17:49 Feature #949 (New): New function: IdealWithIncrementalGBasis(I, f)

(from SC^2)
We could (easily?) implement a function doing an incremental GBasis: if I has GBasis GB, and we add a po... Anna Maria Bigatti

16:06 Feature #947: IsRadical for ideal?

John Abbott wrote:
> I wonder if a probabilistic approach could work. If the ideal is 0-dim, the speed is adequate.... Anna Maria Bigatti

16:02 Feature #947: IsRadical for ideal?

I wonder if a probabilistic approach could work. If the ideal is 0-dim, the speed is adequate. If not, perhaps we c... John Abbott

15:51 Feature #947 (In Progress): IsRadical for ideal?

The crude idea was just to run the usual code for computing the radical, and if at some point we determine with certa... John Abbott

15:19 Feature #947: IsRadical for ideal?

Apart from 0-dim (already implemented), and monomial, how can you determine it?
Anna Maria Bigatti

15:03 Feature #947: IsRadical for ideal?

There is now "IsRadical" for 0-dimensional ideals. (and it's in our paper in progress ;-)
I've just realized I forg... Anna Maria Bigatti

14:56 Feature #947 (In Progress): IsRadical for ideal?

Werner and Mario ask whether it could make sense to have an @IsRadical@ function for ideals; can this be determined s... John Abbott

15:12 Feature #796: CoCoALib function for radical (or SqFree) of a polynomial

This would be quite useful. For the time being should we just add... Anna Maria Bigatti

13 Oct 2016

11:31 Support #944: Release: CoCoALib-0.99550 (together with CoCoA-5.2.0)

Thanks John! Now we have only the "16bytes"... Anna Maria Bigatti

12 Oct 2016

18:08 Support #944: Release: CoCoALib-0.99550 (together with CoCoA-5.2.0)

I'll do @test-OrderedRing2@.... there is a leak in twinfloat.
I know where it is; now I just have to find a nice way... John Abbott

18:02 Support #944: Release: CoCoALib-0.99550 (together with CoCoA-5.2.0)

Anna Maria Bigatti wrote:
> First run of valgrind (I had to reinstall valgrind on my mac...)
> [...]
The "16-byt... Anna Maria Bigatti

17:49 Support #944: Release: CoCoALib-0.99550 (together with CoCoA-5.2.0)

First run of valgrind (I had to reinstall valgrind on my mac...)... Anna Maria Bigatti

16:10 Support #944 (In Progress): Release: CoCoALib-0.99550 (together with CoCoA-5.2.0)

Anna Maria Bigatti

16:06 Support #944: Release: CoCoALib-0.99550 (together with CoCoA-5.2.0)

John Abbott wrote:
> Add also to check whether the CoCoA-5 tests pass with @BackwardCompatible.cpkg5@
OK, moved i... Anna Maria Bigatti

15:42 Support #944: Release: CoCoALib-0.99550 (together with CoCoA-5.2.0)

Add also to check whether the CoCoA-5 tests pass with @BackwardCompatible.cpkg5@
John Abbott

15:04 Support #944 (Closed): Release: CoCoALib-0.99550 (together with CoCoA-5.2.0)

Everything related to making a CoCoALib release:
# Redmine Roadmap: close or postpone issue
# Redmine Release iss... Anna Maria Bigatti

10 Oct 2016

21:20 Support #942: Which names to use? Intersection/saturation vs intersect/saturate

I have put this issue under CoCoALib, but it is just as relevant to CoCoA-5.
In general my personal preference is ... John Abbott

21:18 Support #942 (In Progress): Which names to use? Intersection/saturation vs intersect/saturate

I notice that we have the following 2 names for functions: *@intersection@* and *@saturate@*.
In CoCoA-3 and CoCoA... John Abbott

08 Oct 2016

22:24 Feature #910: BigRat: read from a string in "decimal" format?

Unfortunately the code for reading "decimal numbers" appears twice: once in *@operator>>@* for @BigRat@ and once agai... John Abbott

22:10 Bug #936 (Feedback): ReadExpr: cannot handle exponent which is not integer literal

I have now implemented the following:
after @^@ there may be an unsigned integer literal or @(@ optional-sign unsign... John Abbott

21:55 Bug #938 (Feedback): ReadExpr: misreads 2/3^4

I have fixed it now, and added a new @test-RingElem5@ for the function @ReadExpr@ (for rational numbers).
Trickier... John Abbott

07 Oct 2016

15:15 Bug #938 (In Progress): ReadExpr: misreads 2/3^4

I see no choice but to make @2/3^4@ read as an expression @div(2, power(3,4))@, and this implies that all rationals e... John Abbott

06 Oct 2016

11:38 Bug #938: ReadExpr: misreads 2/3^4

It is a slightly subtle bug: at first sight it seems reasonable that @<constant>^<integer>@ should yield the obvious ... John Abbott

11:37 Bug #938 (Closed): ReadExpr: misreads 2/3^4

After my change to @ReadExpr@ so that it can read rational numbers, the fn now misreads @"2/3^4"@ and produces 16/81 ... John Abbott

08:30 Feature #937 (In Progress): LinDepMill: Mill for linear dependencies

Currently in MinPoly.C there is a class IncrLinDep.
It is small and neat: I think this is the perfect paradigm for a... Anna Maria Bigatti

07:47 Bug #936: ReadExpr: cannot handle exponent which is not integer literal

John Abbott wrote:
> What exactly do we want to allow in the exponents?
> The simplest solution would be to allow e... Anna Maria Bigatti

05 Oct 2016

17:12 Bug #936: ReadExpr: cannot handle exponent which is not integer literal

The relevant source is around line 100 of @RingElemInput.C@
John Abbott

17:12 Bug #936: ReadExpr: cannot handle exponent which is not integer literal

What exactly do we want to allow in the exponents?
The simplest solution would be to allow either an integer liter... John Abbott

17:08 Bug #936 (Closed): ReadExpr: cannot handle exponent which is not integer literal

The fn @ReadExpr@ fails if input contains an exponent which is not an integer literal: for instance... John Abbott

16:40 Feature #910: BigRat: read from a string in "decimal" format?

I'm not sure how related this is: a twin-float value can be printed out with a trailing @<dot>@ if the value is very ... John Abbott

03 Oct 2016

17:29 Bug #935: ElimMat, ElimHomogMat: complain about zero or negative weights

It seems that the weights must be integer. This should be documented!
John Abbott

17:28 Bug #935 (Closed): ElimMat, ElimHomogMat: complain about zero or negative weights

The fns @ElimMat@ and @ElimHomogMat@ should give error if the user supplies weights are not strictly positive (intege... John Abbott

30 Sep 2016

22:11 Design #934: MachineInt: change semantics?

I tried on my old MacBook with g++-4.2, and the results were disappointing. The new class @INDEX@ was significantly ... John Abbott

15:01 Design #934: MachineInt: change semantics?

In the example code above, initially I used the name @index@ for the class (rather than @INDEX@) but that produced st... John Abbott

14:21 Design #934: MachineInt: change semantics?

Here is the test code I have been using: (functional rather than elegant)... John Abbott

14:12 Design #934: MachineInt: change semantics?

For the record, here are the timings I obtained (on the Linux tower in my office):
built-in @IdentityMat@ 1.66s
@Id... John Abbott

13:46 Design #934: MachineInt: change semantics?

I implemented a prototype, and ran a speed test. The results are mixed: usefully faster than current @MachineInt@, b... John Abbott

11:38 Design #934: MachineInt: change semantics?

Advantages:
* One hope is that the simpler implementation may allow the compiler to be clever (and optimise away alm... John Abbott

11:28 Design #934 (In Progress): MachineInt: change semantics?

The implementation of @MachineInt@ is a bit messy as it tries to guarantee that no information will be lost (_e.g._ a... John Abbott

13:31 Feature #813: Implement "elim" in CoCoALib

This should be easy once elim orderings are properly implemented, right?
John Abbott

11:20 Design #933 (In Progress): Separate configure scripts for CoCoALib and CoCoA-5

Currently there is a single @configure@ script for both CoCoALib and CoCoA-5. Strictly they should be separate.
I... John Abbott

11:06 Design #932: CoCoALib configuration: BOOST dependency

A reasonable solution might simply to be to remove the BOOST dependent code in @ThreadsafeCounter@. This would mean ... John Abbott

11:01 Design #932: CoCoALib configuration: BOOST dependency

As far as I can see, the code in @ThreadsafeCounter@ will probably become much simpler when we switch to C++11 (since... John Abbott

10:55 Design #932 (Closed): CoCoALib configuration: BOOST dependency

Currently we have a single configuration process (script) for both CoCoALib and CoCoA-5. This is not ideal!
Gereo... John Abbott

10:46 Feature #319: BOOST -- how it could help in CoCoALib

An aspect to bear mind if we do opt to become BOOST dependent is that linking would then require linking (first) with... John Abbott

29 Sep 2016

19:01 Design #822: Should ElimMat return a ConstMatrix

John Abbott wrote:
> Should @ElimMat([1,2,3], 4)@ and @ElimMat([3,2,1], 4)@ give the same result? Why or why not?
... Anna Maria Bigatti

17:15 Design #822: Should ElimMat return a ConstMatrix

Should @ElimMat([1,2,3])@ and @ElimMat([3,2,1])@ give the same result? Why or why not?
In other words, does the o... John Abbott

16:04 Design #805 (Closed): New type for "constant" matrices?

John Abbott

16:03 Design #592 (Closed): Review design of ConstMatrixView

I'm no longer sure what worried me so much about @ConstMatrixView@; I suspect that the new class @ConstMatrix@ has pr... John Abbott