Activity
From 10 Nov 2015 to 09 Dec 2015
09 Dec 2015
-
22:08 CoCoALib Bug #834 (In Progress): Fix test failures (after revising MatrixForOrdering)
- First bug trapped:
IsPositiveOrd wants @rk(M) == NumRows(M)@.
This function is called when making a term ordering o... -
13:41 CoCoALib Bug #834: Fix test failures (after revising MatrixForOrdering)
- JAA believes the root cause is the check that a grading is strictly positive when weakly positive is sufficient.
-
- The current CVS produces 3 failed CoCoALib tests:
@test-MatrixForOrdering1, test-RingWeyl1, test-SparsePolyRing1@
a... -
- JAA will try to implement this before xmas (2015!)
But first some Gluehwein ;-) -
- I wonder whether @CoCoA_THREADSAFE@ is to be just a temporary measure. Presumably at some point practically all comp...
-
- I spoke to Mario this morning about threadsafety in CoCoALib (also thinking about Mario's new code which has just bee...
-
- Bearing in mind the suggestion in issue #142 about using a compile-time switch to indicate whether we want threadsafe...
-
- See the examples in issue #784; I think this should cover many cases.
-
- Adapt Mario's new code so that it is threadsafe (when @CoCoA_THEADSAFE@ is set).
-
- Currently a @MachineInt@ is converted to a @long@ using the function @AsSignedLong@ which includes a check for overfl...
08 Dec 2015
-
- The fact that UIBC needs to be visible in @RingWeyl.C@ is "unfortunate". John and Mario think this may be a symptom ...
-
- With Mario's new code the file @RingWeyl.C@ must include @TmpUIBC.H@ otherwise it does not compile.
This is likely... -
- I have renamed the fns and changed all calls -- well, everything compiles.
Will check doc, then check-in.
-
- A simple interface is: a function which forces computation of a GBasis (without copying the result), and a function w...
-
- Here is a brief "one-way discussion" about some function interfaces: I shall consider just @IsMonomial@, but it is a ...
-
- With Mario's help we have integrated his code into (my copy of) the CoCoA sources.
Some questions have arisen whic...
04 Dec 2015
-
- After skype with Anna: we prefer the nicer sounding names.
So I'll change the names in CoCoALib so that *@Mat@* is a ... -
- John Abbott wrote:
> Perhaps @MakeTermOrd@ or @MakePositiveTermOrdering@, or something similar?
@MakeTermOrd@! -
- I do not much like the names @CompleteToOrd@ and @NewMatCompleteOrd@.
Perhaps @MakeTermOrd@ or @MakePositiveTermO... -
- Currently matrices for gradings must be full rank and have no negative entries.
In princicple we could accept negati...
03 Dec 2015
-
- The CoCoA-5 tests pass as they did before; and the speed test with 20000 indets seemed to go OK. So marking as resol...
-
- I didn't sleep so well this last night (ate something which didn't agree with me), so today's a good day to hack away...
-
- One possible approach is to devise a hash for @symbol@ then just comparing hashes will probably show quickly that the...
-
- Given the cost of converting between a "compressed" format and a "vector" format, it is probably better to represent ...
02 Dec 2015
-
- I had a quick look on Wikipedia for information about sparse matrices. There they distinguished between representati...
-
- A very simple speed test: compute the order of the field/group generator by naive repeated multiplication.
Prototype ... -
- Here are some options for implementing extensions of (prime) finite fields:
* *(A)* use a CoCoALib quotient ring
* ... -
- I have just tried a simplistic speed test (on the "new" computer here in Kassel). The aim was to see how costly it i...
-
- Certainly an explicit quotient ring which is a finite field will be isomorphic to one of "my" finite fields. It migh...
-
- Some (more or less) obvious comments:
* the advantage of having log/exp tables is that some computations can be noti... -
- @IdentityMat@ and @ZeroMat@ are now in @matrix.H@ -- Anna's suggestion.
The args are now @MachineInt@ instead of jus...
01 Dec 2015
-
- I did remember doing something like this.
This was because of a problem in cocoa 4 about 0 in the first row.
I th... -
- I have put this under CoCoALib even though the example comes from the C5 tests, because I think that te real problem ...
-
- The C5 test @primary.cocoa5@ fails because it calls @CompleteToOrd@ on the following matrix:...
-
- Refering to comment 4: I have now removed the signature with 3 args (after speaking to Anna about it).
-
- I shall shortly try to implement the semantics given in comment 4 above.
After speaking to Anna, the idea is that ... -
- I'm mostly interested in a solution which will make most code easy to read (without significant run-time overhead).
... -
- I've just spoken to Mario about this issue (even though he wasn't on the watchers list).
He made some suggestions:
*... -
- It might make sense to allow other types of "condition":
@RangeCheck(arg, PositiveLong)@ or perhaps simply @RangeCh... -
- Mmm, I guess Scott Meyers would not be impressed by a fn @IsInRange@ which takes 3 integer args -- one could guess th...
-
- Checked in the updated code.
Still need: decide fn names, write proper doc!
-
- The documentation for @ElimMat@ and @HomogElimMat@ is largely missing; and the names given in the doc are different f...
-
- I suppose strictly the fns should be called @ElimOrdMat@ and @HomogElimOrdMat@.
Is that a good idea or a not-so-good... -
- Anna suggested change the arg order for @ElimMat@ and @HomogElimMat@ so that the list of indets to elim is the first ...
-
- There are two ways to create a "lex" ordering for @N@ indets:
@NewLexOrdering(N)@ or @lex(N)@.
I wonder whether... -
- I've made the change (plus some very minor cleaning), incl changing the arg names!
Also made the analogous chang... -
- Implemented. No need to change doc (I believe).
A few minor consequential changes.
Will check in shortly.
-
- The same applies to @SetEntry@, @SwapRows@, @SwapCols@, @DeleteRow@, @DeleteCol@.
-
- The pseudo-ctors @ZeroMat@, @IdentityMat@, @MatByCols@, @MatByRows@ have params of type @long@; these should be of ty...
-
- Since the fn @AssignZero@ (for matrices) was no longer used anywhere, and especially since its semantics were rather ...
-
- Here are some other names.
With *SUFFIX* @Mat@:
@ZeroMat@, @IdentityMat@, @BlockMat@, @DiagMat@, @ColMat@, @RowMa...
30 Nov 2015
-
- @IsGreaterThan@?
Or just overwrite @operator>@? -
- Ok, for the Check function (with the function name).
But I 'd rather have (n, lo, hi)... Am I too late to notice ... -
- In prima battuta farei solo @MaxLong@ perche' voglio "scoraggiare" l'uso di @unsigned long@.
Non ti viene i mente ... -
- John Abbott wrote:
> Perhaps the existing @IsInRange@ is not too far from what I want; the real hitch is that an upp... -
- Perhaps the existing @IsInRange@ is not too far from what I want; the real hitch is that an upperbound has to be spec...
-
- Many fns accepting a @MachineInt@ actually want a positive (or non-negative) value.
There is a fn called @IsInRang... -
- Ok, good for operator()!
Still I 'd call the variable @OrdCtor@ ;-) -
- Yes, it should be possible to write @lex(NumIndets)@, but normally it would be enough to write just @lex@ and let the...
-
- Would it be possible to write @PPOrdering O = lex(NumIndets);@?
If not, I fear it might be more confusing.
In fact... -
- We could also have both, but I do not usually like to offer two ways of achieving the same thing... a user might be p...
-
- Currently @PPOrderingCtor@ has a mem fn called @myCtor@ which accepts an integer arg and returns the ordering for tha...
-
- John Abbott wrote:
> @GradingMat@ is used in @hp.cpkg5@ but just for testing whether the grading is "standard"; Ann... -
- John Abbott wrote:
> Does the 0x0 matrix define a term-ordering?
> Do we want to allow users to make a term-orderin... -
- Does the 0x0 matrix define a term-ordering?
Do we want to allow users to make a term-ordering on 0 indets?
JAA think...
27 Nov 2015
-
- Maybe I can do this next week?
-
- Would it be possible to have this by summer 2016?
-
- Just spoken to Anna (via Skype) about the idea of non specifying the ring.
It seems safer for @IdentityMat@ and @Z... -
- @GradingMat@ is used in @hp.cpkg5@ but just for testing whether the grading is "standard"; Anna can you fix this?
-
- I suggest that you try with @degree@ objects containing just a single component, that way the saving as a percentage ...
-
- I note that @test-RingWeyl1.C@ tries to create a an elim mat for no indets. Why? Is that really what it should do?
-
- What should the @ElimMat@ function do in these two limit cases?
* elim no indets
* elim all indets
Note that the...
26 Nov 2015
-
- If there truly is a need to convert a matrix of rationals into a matrix of integers defining the same term ordering, ...
-
- Also converting a (non-neg?) matrix of rationals into a matrix of integers defining an equivalent ordering allows for...
-
- Perhaps the first question should be: do people who build order matrices really start from matrices with rational (no...
-
- In @test-MatrixForOrdering1.C@ the function @NewPositiveMat@ is called with a matrix containing rational (non-integer...
-
- I am pretty convinced that @ZZ@ is the right ring for order matrices in CoCoALib.
That said, I am making sure that... -
- John Abbott wrote:
> The CoCoALib name is @NewMatCompleteToOrd@ while in CoCoA-5 it is @CompleteToOrd@.
> Is this w... -
- John Abbott wrote:
> I suggest that the semantics of @NewMatCompleteToOrd@ be the following:
yes, I agree.
I'm j... -
- I suggest that the semantics of @NewMatCompleteToOrd@ be the following:
* expect as input a matrix @M@ (over ring of... -
- The CoCoALib name is @NewMatCompleteToOrd@ while in CoCoA-5 it is @CompleteToOrd@.
Is this what we want?
-
- What exactly should @NewMatCompleteOrd@ do?
What does its expect as input? And what does it give as output? When d... -
- I have hacked the code: assuming all matrices are over @ZZ@ makes the code a lot simpler. I have added a fn which cr...
-
- Putting into Feedback now that Anna is happy with the C5 manual entry.
I'll add @ILogBase@ to @obsolescent.cpkg5@.
-
- John Abbott wrote:
> I'm still happy with the single combined page for @FloorLog2@, @FloorLog10@, and @FloorLogBase@... -
- I think it is useful. I'll investigate and document it
-
- John Abbott wrote:
> In some ways it would be nice to have a default impl of @AssignZero@; presumably as the mem fn ... -
18:04 CoCoALib Design #825: IsPositiveGrading -- really need 2 signatures?
- John Abbott wrote:
> Just to clarify the semantics:
> * @IsPositiveGrading@ accepts a matrix of integers (error if ... -
- Just to clarify the semantics:
* @IsPositiveGrading@ accepts a matrix of integers (error if entries are not integers... -
- John Abbott wrote:
> What is the difference between @IsTermOrdering@ and @IsPositiveGrading@?
>
> One allows a ze... -
- What is the difference between @IsTermOrdering@ and @IsPositiveGrading@?
One allows a zero-column, the other not?
-
- The fn @IsPositiveGrading@ has 2 signatures:
@IsPositiveGrading(M)@ and
@IsPositiveGrading(M,GrDim)@
Do we reall... -
- Mario now thinks he wants sparse matrices.
We need a clear description of the requirements. -
- JAA thinks that names such as @StdDegRevLexMat@ read more easily than @MatStdDegRevLex@, but the names @MatXYZ@ all a...
-
- There is an incompatibility between the names in CoCoALib and in CoCoA-5 for the pseudo-ctors for matrices for orderi...
-
- Following the principle of "no nasty surprises", I think that @IsTermOrdering@ should do a complete check (incl. corr...
-
- JAA is now convinced that having @OrdMat@ return a "reference" (namely a @ConstMatrixView@) is correct: it was relati...
-
- Could it be useful to have a pseudo-ctor which accepts just the "grading matrix"? The pseudo-ctor would then "auto-c...
-
- Now I wonder why there is the first version of the pseudo-ctor:...
-
- An alternative would be to accept a @record@ as an argument, but this does make the input longer: for instance...
-
- According to the C5 manual, if we want to create a new poly ring with given weights then we must supply a complete or...
-
- The original report did not contain an example of input which produces the unhelpful error message. I think the foll...
25 Nov 2015
-
- I really feel I have opened Pandora's box :-(
I'll keep plodding on... this is all very frustrating :-(
-
- In my new implementation the fixed order matrices return @ConstMatrix@.
At the moment @ElimMat@ returns a @matrix@; ... -
- These matrices for ordering really seem to be a minefield :-(
Who uses @NewIntegerOrdMat@?
Unfortunately I need... -
- @MatrixForOrdering@ offers a fn called @IsTermOrdering@.
The doc says its tests whether the matrix defines a term or... -
- JAA really thinks that order matrices should be over @ZZ@ rather than @QQ@ since we insist that the entries be intege...
-
- There are 2 pseudo-ctors for @NewMatrixOrdering@:...
-
- I have added the check to the ctor for @MatrixOrderingImpl@, and also added a new test @test-bug5.C@. I'll check the...
-
- If the input matrix does not have integer entries then the pseudo-ctor should give error (it may delegate the task to...
-
- The following example compiles and runs, and prints out a matrix with rational entries....
-
- I am getting ready to check in the changes I have made, and a few oversights are coming up :-(
Now each concrete @... -
- The function @GradingMat@ bothers me a bit. Is it actually useful?
(yes)
Can it be removed?
(no)
24 Nov 2015
-
- I have another doubt/question.
The pseudo-ctors for matrices for orderings do not accept a ring as an arg; the mat... -
- I am unsure where the put the new defns of @ZeroMat@ and @IdentityMat@.
They should not be in @MatrixView.C@ since... -
- I now think that the old code was "safe" (_i.e._ no risk of slicing etc). However it was definitely confusing that t...
-
- I have found a reasonable way to circumvent the problem of passing a @ZeroMat@ (of type @ConstMatrix@) to @BlockMat2x...
-
- I think we should implement *(B)* with the "LCM" meaning from comment 4.
I can do this after checking in the curre... -
- John Abbott wrote:
> A (minor) disadvantage of having to pass in a vector of @symbol@ when creating a @PPMonoid@ is ... -
- A (minor) disadvantage of having to pass in a vector of @symbol@ when creating a @PPMonoid@ is that the monoid is obl...
-
- Should the @NewPolyRing@ fns automatically choose @PPMonoidSparse@ if there are "lots of" indets? If so, how many ar...
-
- Cleaned the @PPMonoidSparse@ code; there were still a couple of incomplete fns.
-
- I have a first (complete) version which seems to work OK -- just simple tests.
I even managed to create a sparse pol... -
- What is the difference between @BackwardCompatible.cpkg5@ and @obsolescent.cpkg5@?
-
- Anna, do you mean "yes" to both proposals?
-
- John Abbott wrote:
> Should there be a function @isqrt@ in the obsolescent package?
> I note that @Isqrt@ is in @Ba... -
- Should there be a function @isqrt@ in the obsolescent package?
I note that @Isqrt@ is in @BackwardCompatible.cpkg5@;... -
- I have now checked in the changes. Anna, could you verify? Thanks!
-
- Not too hard. Renamed @isqrt@ to @FloorSqrt@ (in both CoCoALib and CoCoA-5); changed tests, examples, doc.
Status... -
- I'm still happy with the single combined page for @FloorLog2@, @FloorLog10@, and @FloorLogBase@. But if you can conv...
-
- The message "Imminent disaster" is entertaining but not so useful to the user (especially if the message has not been...
-
- I awoke this morning with a plausible explanation of the problems with @test-F5@.
The problem is that the dtor for...
23 Nov 2015
-
- Arri's life is getting shorter by the minute -- grrr!
There is a 3 arg definition in @TmpF5Mat.C:867@, but the 3rd a... -
- With 3rd formal arg commented out, it compiles and runs fine, but produces a different answer from before... :-/ At ...
-
- So why did this problem not appear earlier? And why did it appear now?
-
- An easy workaround is just to comment out the 3rd arg in the declaration.
Since no one ever calls it with a 3rd arg ... -
- If Arri comes my way I'll shorten his life by at least 2 hours -- the time I've wasted tracking this bug down. Grrr!...
-
- In @TmpF5.H@ two fns called @F5@ are declared: one (inline) with 2 args, and one with 3 args (3rd arg has a default v...
-
- Checked in the renaming @ILogBase@ to @FloorLogBase@.
Currently the C5 manual page is single; it's what I'd done, an... -
- Same should be for @isqrt@! @FloorSqRt@? @FloorSqrt@?
-
- John Abbott wrote:
> Should the CoCoA-5 manual pages for @FloorLog2@ and @FloorLog10@ be separate from that for @Flo... -
- I've almost finished making the changes. I think @FloorLog2@ and @FloorLog10@ do improve readability of the code whe...
-
- John Abbott wrote:
> Personally I think that *@FloorLogBase@* is too long. Anyway, the fn should expect two args (o... -
- At the Dagstuhl meeting last week some people were interested in having @GenRepr@ available in CoCoALib. I said it w...
-
- Rename isqrt to FloorSqrt (following ilog renaming to FloorLog)
-
- John Abbott wrote:
> What about @NewPPMonoid("x", 10)@ as a shorthand for @NewPPMonoid(SymbolRange("x",0,10-1), StdD... -
- What about @NewPPMonoid("x", 10)@ as a shorthand for @NewPPMonoid(SymbolRange("x",0,10-1), StdDegRevLex)@?
Or should... -
- John Abbott wrote:
> We could slightly modify Mario's request:
> have new pseudo-ctors which do not need the orderi... -
- I note that in issue #320 the main reason for rejection is that J and A disagreed about what symbol names should be u...
-
- Mario asked about having an easy way to create a @PPMonoid@ without specifying explicitly the ordering (or even the n...
-
- I have put the lines into a separate (static member) function. What should it be called?
Currently I have called ... -
- The real problem with @terminate@ is that it gives the program no chance to handle the situation; in fact the only "b...
-
- The real problem is that the check for still-existing @RingElem@ values occurs inside the dtor for @GlobalManager@, a...
-
- I've just made some big changes to matrices for orderings: everything compiles, but @test-F5@ produces "imminent disa...
-
- During destruction @GlobalManager@ checks the ref counts of @RingQQ@ and @RingZZ@; if they are too large then "Immine...
-
- I notice only now that these pseudo-ctors take only the dimension as arg; *what should the ring of the matrix be?*
... -
- I think that @NewMatrixOrdering@ does not check the ring of the matrix handed to it.
Check if I am right; anyway, ... -
- In CoCoALib we have member function @myElim@.
Choose and implement the "proper" function.
-
- John Abbott wrote:
> Checked in. Should be OK, but I do have some other files "pending", so it is possible I missed... -
- John Abbott wrote:
> Did you also try tests which create lots of ring elements? Maybe a non-trivial GBasis?
I tr... -
- We could simply merge the code in @PPMonoidEvZZ@ into @PPMonoidEv@.
What do you think?
Are they similar enough th...
21 Nov 2015
-
- Checked in. Should be OK, but I do have some other files "pending", so it is possible I missed something -- can you ...
-
- Implemented; documented (also C5); added test (NumTheory3)
-
- In some applications it is helpful to use "binary rationals" (denom is a power of 2).
Add a new fn *@SimplestBinaryR... -
- After the SC^2 meeting it is clear that several people would like @RealRoots@ to be ported into C++. If the project ...
-
- Did you also try tests which create lots of ring elements? Maybe a non-trivial GBasis?
20 Nov 2015
-
- John Abbott wrote:
> Change @MinLoafBytes@ to 4*1024. Maybe you could also try 1024.
> Change @MaxLoafBytes@ to 64... -
- John Abbott wrote:
> Are the problems my fault? Or some problem with the interface of @MemPool@? :-(
no, my pro... -
- Are the problems my fault? Or some problem with the interface of @MemPool@? :-(
-
- Personally I think that *@FloorLogBase@* is too long. Anyway, the fn should expect two args (of which the second is ...
-
- I was wondering about using a better name for @ILogbase@.
I used it recently in a new function, and couldn't remembe...
19 Nov 2015
-
- about to try..
(some compilation problems...) -
- Is the function @FastCmp@ for @degree@ actually useful?
Check whether it is perceptibly faster than @cmp@; if so, ...
18 Nov 2015
-
- The following variant works OK:...
-
- The following produces an NYI error for the reciprocal:...
14 Nov 2015
-
- I've still got problems with the new (internal) design for @BlockMat2x2@. While investigating I was quite surprised ...
-
- As a general rule we say that mem fns do not perform arg sanity checks because such checks are the responsibility of ...
-
- The mem fn @mySetEntry@ for @DiagMat@ checks that the indices indicate a diagonal element only if @CoCoA_DEBUG@ has b...
-
- In some ways it would be nice to have a default impl of @AssignZero@; presumably as the mem fn @MatrixView::myAssignZ...
-
- Before spending too much time deciding the exact semantics, we should also decide if the procedure is even useful -- ...
-
- There is a procedure @AssignZero@ for a matrix; it calls the mem fn @myAssignZero@.
What should these fns do? -
- I've implemented the other defn of @StdDegRevLexMat@; we can decide later which impl will remain "active".
13 Nov 2015
-
- As predicted in comment 1, there is a problem with @BlockMat2x2@. Not yet sure how to deal with it :-/
-
- In principle the concrete classes for @IdentityMat@ and so on could derive from @matrix@ but then they have to impl l...
-
- I believe that @ConstMatrixView@ is not the right type, because a matrix "view" is a way of looking at another object...
-
- I believe that implementing @IdentityMat@, @StdDegRevLexMat@ and so on as objects of type @ConstMatrixView@ is not co...
-
- What do we want the @StdDegRevLexMat@ to be?
There are two obvious candidiates: here are 3x3 instances...
12 Nov 2015
-
- Rather to my surprise I have failed to create an example where @ZeroMat@ (or @IdentityMat@) misbehaves -- no idea why...
-
- I think that @ZeroMat@ and @IdentityMat@ should produce a result of type @matrix@ and not one of type @ConstMatrixVie...
11 Nov 2015
-
- After some hasty cut-and-paste work I have some code which compiles :-)
I'm not so foolish as to try running it... t... -
- JAA is convinced that this is a good idea.
To get it finished quickly I shall ignore issue #592; it probably shoul... -
- JAA now thinks that there is an example (see issue #803) where it could be useful to have the order matrix(view) quic...
-
- How important is it that @OrdMat@ is fast? Is there a good example where we would want to compute repeatedly the @Or...
-
- The main task I want to simplify is that of implementing @myWDeg@, @myCmpWDeg@ and @myCmpWDegPartial@ for @PPMonoidSp...
-
- Should objects of type @PPOrdering@ be able to compute WDeg vectors and/or individual components of WDeg vectors?
... -
- I have lowered the priority of this issue because it will take some time to complete, and the resulting improvement w...
-
- As far as I can tell @PPMonoidBase::myIndets@ is used only once, in the (inefficient) fn @PPMonoidBase::mySymbolValue...
-
- John Abbott wrote:
> I think updating would mean replacing the div-mask with that for the LCM of the old value and t...
10 Nov 2015
-
- John Abbott wrote:
> Might it be important to store the @StdDeg@ along with the list of (index, exp) pairs? If it i... -
- For @Wdeg@ I am inclined to compute it "on the fly" every time it is needed. If someone reports poor execution speed...
-
- The currently missing functions are: @myCmp@, @myWDeg@ and @myCmpWDeg@ (& a partial version). Also @myComputeDivMask...
-
- There was a very incomplete impl in @PPMonoidSparse.C@, which was already recognized by @Makefile@.
Apparently the... -
- I think updating would mean replacing the div-mask with that for the LCM of the old value and the new indet-power. T...
-
- John Abbott wrote:
> If we do extend the interface of @DivMask@ to accommodate the sparse repr, what should the exte... -
- If we do extend the interface of @DivMask@ to accommodate the sparse repr, what should the extended interface expect:...
-
- Currently the only way to make a (non-trivial) @DivMask@ is to use the function @myAssignFromExpv@ which clearly suit...
-
- The new @PPMonoidSparse@ impl needs to supply a @myComputeDivMask@ mem fn.
The current interface is well suited to... -
- That's fine by me -- less work for me :-)
Indeed @symbols("x[1..3,2..5]")@ is not really that much more readable tha...
Also available in: Atom