Activity
From 02 Nov 2015 to 01 Dec 2015
01 Dec 2015
-
22:39 Bug #831: primary.cocoa5
- I did remember doing something like this.
This was because of a problem in cocoa 4 about 0 in the first row.
I th... -
19:00 Bug #831: primary.cocoa5
- 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 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:...
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...
09 Nov 2015
-
- John Abbott wrote:
> Do we want this function to handle symbol ranges?
>
> It probably makes sense for @symbols... -
- Do we want this function to handle symbol ranges?
_e.g._ @symbols("x[1..3,2..5]")@ would produce the same result as ... -
- Currently this operation is implemented in a hidden way inside @RingBase::myNew(symbol)@.
The implementation there... -
- Should there be a function which tests whether a given @symbol@ is in given @ring@?
I am thinking of "recursive" s... -
- I had assumed that this issue already existed, but a quick search did not find it.
There is already a partial impl... -
- Implement a PPMonoid which uses a sparse repr for the PPs.
This is useful when there are many indets.
-
- The vector of indets is used in the (mem) function @mySymbolValue@, which currently returns a reference to the image ...
-
- Is it truly useful to have each PPMonoid keep a vector of its own indets?
This is not a problem if the indets are ...
07 Nov 2015
-
- The speed gains varies form platform to platform -- surprise, surprise!
The old @DUPFF@ code was (significantly) f... -
- After comparing @gcd@ for @DUPFp@ and @DUPFF@, I saw that the latter was significantly faster (_e.g._ factor 2 or mor...
05 Nov 2015
-
- This problem is generally easier than @SqFreeFactor@, indeed in @SqFreeFactorPosDerChar0@ it is just computed by the ...
-
- Anna Maria Bigatti wrote:
> John Abbott wrote:
> > I'm not sure; I suspect it will not work in every case.
>
> I... -
- Good point! So if @gcd@ works then we can expect that squarefree decomposition works.
-
- John Abbott wrote:
> I'm not sure; I suspect it will not work in every case.
I thought it would: in char 0 the squar... -
- The *@Fatal error@* messages appeared once again when I tried building CoCoALib from a new unpackaged tar file. The ...
-
- Implemented, documented. No separate test; it is effectively tested inside @test-NumTheory1@
-
- For efficiency it would be nice to have public access to the "fast" function for computing inverses modulo an integer.
-
- [oops cut-and-paste did the wrong thing]
Changing to feedback; implemented the fn mentioned in comments 18-20. See #795 -
- John Abbott wrote:
> I am thinking of the name @InvModNoCheck@; it is pretty clear; it is a bit long, so that ought ... -
- I am thinking of the name @InvModNoCheck@; it is pretty clear; it is a bit long, so that ought to discourage "casual ...
-
- In @NumTheory.H@ there is a fn called @InvMod@ which expects two @MachineInt@; internally it calls the fn @InvMod_pri...
-
- I have resolved the mystery of the two identical functions running at different speeds: my test program was giving (v...
04 Nov 2015
-
- I have implemented @NonRedValue@.
Checked everything last night; revised examples, doc, etc.
Checked in today.
T...
03 Nov 2015
-
- I'm not sure; I suspect it will not work in every case.
What coeff ring did you have in mind? Were you thinking o... -
- In @factor.C@ there is the line:...
-
- I'm compiling some other programs; compilation time varies from 2mins to almost 5mins. I'll ask the technical guy if...
-
- I'm recompiling (with full debugging), and the mysterious @Fatal error@ messages have not reappeared; and anyway desp...
-
- John Abbott wrote:
> The new computer took *100mins* to compile and run the examples. Just running the examples take... -
- The new computer took *100mins* to compile and run the examples. Just running the examples takes only 30secs.
-
- I shall try to investigate this, but it may take a long time. The new computer seems to be quite fast (_e.g._ more t...
-
- The archive creator @ar@ gives warnings about the use of @u@ flag with the @r@ flag (because the default is @D@).
... -
- I quite like @NonRedValue@; it is quite clear, and a bit long (to discourage "casual use"). I'll try changing the co...
02 Nov 2015
-
- John Abbott wrote:
> There is a question about names:
> * *@SmallFpImpl::value@* is the class for reduced values
>... -
- Compilation on the new computer (with fedora 23, due for release tomorrow ?!!? according to Wikipedia) produced sever...
-
- Since I believe that the @find@ command is the culprit, the solution would be to eliminate it. This could be done by...
-
- The situation was the following. The new computer did not have GMP installed, so the attempt to use @-lgmp@ failed, ...
-
- I have new computer (here in Kassel). I am trying to build CoCoALib/CoCoA-5. The @configure@ script took *ages* to ...
Also available in: Atom