Dave Keenan wrote: ↑Tue Jun 30, 2020 1:07 pm
"d" is problematic, as I've used it for denominator many times above. Can we call it "w"? My earlier use of "w" for what you're now calling "c", was short lived ... I'm happy to go back and edit my earlier copfr coefficients from "w" to "c".
You're right, and I should have noticed that. Yes, let's make it w. I can go back and edit my previous post from "d" to "w" too.
Dave Keenan wrote: ↑Tue Jun 30, 2020 1:32 pm
I'm afraid this ... is really pushing our idiosyncratic function-naming scheme beyond the breaking point.
...we could always fall back to the real name
Dave Keenan wrote: ↑Tue Jun 30, 2020 12:07 pm
BTW, changing that variable name from α to e is a really bad idea, Mexican food notwithstanding, because using e as a variable log base is exactly the place where it would be confused with the constant log base e ≈ 2.718.
[Edit: Oops! You suggested epsilon ε, not e. That wouldn't be so bad, except it is traditionally used to represent an infinitesimal.]
I'm not sure what "real" name you're referring to. As soon as we went beyond sopfr and sopf, aren't these all new things without established/"real" names? Or do you mean something else by that?
These unwieldy names are helpful during the development process. Once we reach the final step of naming these functions nicely for the outside world we can rely a lot more on the descriptions of what exactly the functions do, and reduce the name to something optimized for pronounceability, e.g. "soapfar" for "sum of adjusted prime factors (with) adjusted repetition", or something like that.
It might be cool if we added a plugin to the forum for LaTeX or MathJax. That might help get these formulas across in a less disgusting and/or intimidating way
Dave Keenan wrote: ↑Tue Jun 30, 2020 12:07 pm
I'm keen to know how low you can get the SoS in the vicinity of these results.
Do you still want me to check those, even though I've since found one with 0.004250806?
Actually, scratch that. I think your method will need to fulfill the role of the "home stretch": finding the exact values down to the millionths place or whatnot. The way I'm doing things, it's not really tractable to look deeper than the thousandths place. So if you're already working with SoS-billionths, I'm not going be able to help you get any more precise.
I could at least double-check them, though, if you want.
Dave Keenan wrote: ↑Tue Jun 30, 2020 12:31 pm
cmloegcmluin wrote: ↑Tue Jun 30, 2020 11:15 am
The simple copfr experiment was not fruitful. I couldn't find any example of a weight > 0 on it helping. Cool idea, but it doesn't look like its our winner.
No, but the fact that a negative weight on it helped, was interesting.
I actually didn't check it with negative values, FYI. I should have tried that, though, and I should have been more clear.
Dave Keenan wrote: ↑Tue Jun 30, 2020 12:31 pm
I realised that if instead of the simple c × copfr(n/d) you split it into
w×copfr
y(n) + k×w×copfr
y(d), then when you add it to the existing
sop
αfr
y(n) + k× sop
αfr
y(d), where p
α = log
α(p), you get
sop
αwfr
y(n) + k×sop
αwfr
y(d), where p
αw = log
α(p)+w.
it was simpler to roll copfr into so{..p..}fr, whereupon your "c" was shown to be almost equivalent to your "d". I'm happy to go back and edit my earlier copfr coefficients from "w" to "c".
The almostness is due to the fact that copfr was not previously split into separate calculations for numerator and denominator, and was not previously applied to prime exponents raised to the power y.
I can't figure out how that works. Is it using a logarithmic identity I'm not familiar with? I'm interested, certainly, since it seems like you found a way to consolidate the count of primes into their sum.
Dave Keenan wrote: ↑Tue Jun 30, 2020 12:31 pm
I don't think it makes as much sense to try that on the term.
Not sure what you mean by that. But I should mention that I also tried weighting the prime exponent
before raising it to the power y, i.e. ((log
α(p)+w)×r)
y, but I couldn't get the SoS quite as low with that as with (log
α(p)+w)×r
y.
Dang, I was afraid that might not be clear enough, but I was in a rush. What I meant was that while it feels right to adjust the primes by some constant – if only because we're starting with 5, just kind of floating out there in Fiveland
– adjusting the terms of the monzos (AKA the repetitions, or "r") by a constant doesn't make as much sense (whereas I
do think adjusting them by a power or log does make enough sense).
Dave Keenan wrote: ↑Tue Jun 30, 2020 1:44 pm
cmloegcmluin wrote: ↑Tue Jun 30, 2020 12:36 pm
d = -2.08 (so basically subtract 2 from each prime)
Do you mean "so basically subtract 2 from the
log of each prime"?
Yes, what I meant was subtracting 2 outside the quotes, like log
α(p) + w. I did
not mean log
α(p + w). In my defense, I'm pretty sure logarithms come before subtraction in order of operations, so my "log
2-p" was accurate (and I think the only reason I left off the parens was because I was trying to have it double as part of something like a name... the problem which we've already covered above)
Dave Keenan wrote: ↑Tue Jun 30, 2020 12:31 pm
I should mention that I also tried weighting the prime exponent
before raising it to the power y, i.e. ((log
α(p)+w)×r)
y, but I couldn't get the SoS quite as low with that as with (log
α(p)+w)×r
y.
Ah. I'm glad for your "i.e." clause because I would have interpreted what precedes it differently otherwise. I thought by "weighting the prime exponent before raising it to the power" you meant log
α(p + w). I thought that because we're still entertaining both sublinear exponents and logarithms (though logarithms seem to be winning now) and so when you said "before raising it to the power" I thought that was standing in for either raising to a power or putting to a logarithmic base. I had not even yet considered the prospect of raising the whole thing to some power (or to some base), even after the repetition count had been raised to some power (or put to some base). All of these things are real and distinct possibilities though! We could have:
((log
α(p + v)+w)×r
x)
y or log
y(((p + v)
α+w)×log
xr)
or any other combination of logs and exps, where v is just some other constant and x is just some other exponent...
Yeah... let's not go down that path...
Dave Keenan wrote:
We're definitely in trunk-waggling territory now, with 5 or 6 parameters. We should try to cut that down if possible without going over SoS = .0055.
I agree with the sentiment here.
Honestly, I think that while "max(so{log2p-2}f√r(num), so{log2p-2}f√r(den)) + ⅗copfr(num/den)" looks kind of gross, it wouldn't actually be that bad once you got in LaTeX. Especially if you have devised a way to fold in copfr. And I think the log base 2 on the prime is kind of pretty, even, in its recognition of human base-2 pitch perception. The square root bit of r makes less immediate musical sense, but it offends me less than the ⅗ on the copfr, insofar as I sense there might be an underlying truth to something like a square root...
Dave Keenan wrote:I was hoping you could set c=0 without too much damage because d (or w?) should pick up the slack. But then I realised that k=0 means the denominator is irrelevant, except in the count of its primes, including repetitions. So if you set c=0 then you will need k≠0.
I will investigate this k=0, c≠0 regime.
I look forward to your results. I'm not really sure what more I can produce tonight until I hear back from you in more detail about your consolidation of copafry into sopafry, so that I can understand it well enough to implement it myself.
Dave Keenan wrote: ↑Tue Jun 30, 2020 2:59 pm
What do you get when you set k=0, c≠0, and what when you set c=0, k≠0?
Welp, you snuck one more in before I could get my post out.
with k = 0, I can get 0.004609100.
that's with c = 0.723, a = 1.753, y = 0.473, w = -2.620
those parameters are disconcertingly different from the c = 0.577, a = 1.994, y = 0.455, w = -2.08 which got us 0.004250806 with k = 0.038. If you just take those parameters and set k = 0 then you get 0.004749566 which is also quite close.
with c = 0, I can only get 0.006251296.
that's with k = 0.635, a = 1.430, y = 0.850, w = -2.770
So it would seem that the right path forward is to obliterate the smaller of the num and den, and use the count of primes to account for the presence of harmonic information on the other side.