@acl I still enjoy gaming the site, at least a little. I'd like to get the Legendary badge even though I can't possible compete with R.M for badge count.
@Mr.Wizard yes, I was half-joking. I can't devote as much time to this nowadays as some time ago, and even if I could, I can't compete with most of you. The grapes of disappointment are always sour :)
(well I thought of it in Greek, but I think this gets close to what I mean)
by the way, nice to see you did not get eaten by pigs
Have some pride in the heritage my man! (re: pigs, I would have expected that kind of a story to originate in a place which has a serious feral pig problem, not here.)
@acl "A pathologist was unable to identify a cause or manner of death, the newspaper reported. The remains will be examined by a forensic anthropologist at the University of Oregon"
The first story I read suggested perhaps a heart-attack; that seems likely, as while I've heard that pigs will eat most anything, I don't think the domesticated ones attack people.
I've also heard of at least one murder case where the killer tried to use this to hide the evidence. As I recall he almost got away with it. I suppose this could be another.
Really? I haven't been on a working farm I'll admit, but I thought that was bred out of them. That's one of the reasons the feral hog problem in the southern US is such a concern.
@J.M. from profiling my code, it appears that I will need to compile Gamma. I'll get to that tomorrow as it's late. If you have something that would be much appreciated as I am sure it'll be better than my attempt.
@J.M. also thanks for the refs on symplectic integrators (and to @belisarius too); we've implemented a symplectic integrator and managed to get much larger system sizes than have been done in this context before, and all in pure mathematica!
proves that a better algorithm beats optimization most of the time
@acl Just to be sure: (modestly sized) nonnegative integers only, or are you expecting to evaluate real arguments as well? You don't need Lanczos if you just want to compute 10!.
@J.M. nonegative integers, and eventually the results will end up in a large matrix the eigenvectors of which will be calculated. thus, everything should fit in floating point. (but not just 10!). the application is enumerating states of a quantum system and obtaining the full set of eigenvectors of a hamiltonian (thus the floating point constraint)
but no, evaluating Gamma for reals is not necessary
@J.M. it might make sense if I precalculate the factorials and use a lookup table, but I am not sure how to do that efficiently in a compiled module. if I was doing everything in C that's what I'd do
maybe I can ask a question on how to inject precalculated lookup tables in compiled functions in an efficient way. but I'll have to think about this a bit first
@J.M. it does, but I'd need to work out how to efficiently inject that into compiled code. maybe it's small enough to literally inject it. I'll play with it (tomorrow, coding anything now is a recipe for disaster)
@J.M. Sorry, @J.M. i am not sure why it did not pass through. I would be happy to read it if you post it again. Maybe something just malfunctioned. Let me know what alias you posting under and when - and I will try to check right away.
@VitaliyKaurov Nah, it's okay. I had wanted to mention that Simon explained his use of that dreaded font in his answer here, but it seems another one of the more recent comments has taken care of it.
@J.M. I have quick question regarding MultinormalDistribution versus NormalDistribution. I seems to me mathematica is partially inconsistent because NormalDistribution expects as a second argument the standard deviation while MultinormalDistribution expects the matrix of variance covariances.
@J.M. as a consequence (I am not 100% sure this is the cause but I have strong suspicion) my code (mathematica.stackexchange.com/questions/10220/…) fails when the conditionnal probability only involves one variable.
as a simpler formal example, this fails to `Clear[pdf]; pdf = MultinormalDistribution[{Subscript[\[Mu], a], Subscript[\[Mu], b]}, {{Subscript[s, a], Subscript[s, c]}, {Subscript[s, c], Subscript[s, b]}}]
@J.M. uhm... it seems adding the Sqrt in `If[Length[\[Mu]2] == 1, NormalDistribution[\[Mu]2 // First, Sqrt[S2 // First // First]], MultinormalDistribution[\[Mu]2, S2]]` does the trick!
I have a text file with little endian UTF16 encoding (Windows, Scite). How on earth do I import this thing into Mathematica? I am about to look for Java code to do it, or does anybody have a better idea?
So what I do now is to manually save a copy of the text file from Scite to UTF8, and then things work fine (I only have russian characters, so right now there is no need for true UTF16), but why does Mathematica not have UTF16? I mean, there is MLGetUTF16String etc. in MathLink, so it seems to be used internally, or?
Yesterday I posted this question which, as Mr.Wizard noted, lacks constructivity. I have to agree with this, and when I started to think on editing I realized that there are actually three options which are of interest for me at least.
Animated process of construction from scratch, as it is des...
I need to implement a new "game of life" in Mathematica: I have three states — black, red and white. Black means go forward, red go back and white take no action.
I would like to implement the outcomes in this Excel file: http://cl.ly/0s2B2M3h3M2N
The code I have used is the following:
input13...
or if a have a central cell that is white, my code say that if the central cell is white and the black neighbors are equal to the grey one, the result should be a white cell
the meaning of my rule is the following: cell red(grey) interact with cell black. if the number of the cell in one color(black/red) is bigger than the cell of the other color (red/black), then the resulting cell will be the one that has the biggest number
for example if you look at a123:c125 you will see:
center cell black (2), one cell black (2), two red (1) and 5 white --> the result will be white
because the two red cell interact with the two black cell and they die, leaving one cell white
if you see A127:C129, there i have one black cell in the center, one black, three red and 4 white --> the two black cell interact with two of the red cell (so they die), and the only cell that doesn't interact is the red cell --> so the result will be the red cell
IF[(all the conditions to have as result 2 are correct; then Return[2]); if false the find all the conditions to have as result 0; then return [0];otherwise Return[1];Return[1]]
Yesterday I posted this question which, as Mr.Wizard noted, lacks constructivity. I have to agree with this, and when I started to think on editing I realized that there are actually three options which are of interest for me at least.
Animated process of construction from scratch, as it is des...
