well at the end of the day, the calculations are going to be done on the GPU. So most of these algorithms are going to be eventually hardcoded. I'm not sure how much that really impacts any of this discussion
doubling my memory requirements and killing my performance by switching to double precision was something I really didn't want to have to resort to... Like as in the performance hit is almost going to be so bad to kill the project
well, then try every choice at the prototyping phase (i.e. matlab and numpy)
try what happens if you compare full double vs partially double vs fully single with singulars removed
by "singulars removed" I mean looking at the spectral decomposition V^(-1)*D*V of your matrix (hopefully it exists), and finding the close-to-zero eigenvalues in D, and treating them separately with their corresponding eigenvectors
although 1/sqrt(M) should not exist for singular eigenvalues...
since sqrt(0)=0
as you probably know, the function f acting on the matrix A=V^(-1)*D*V is equal to f(A)=V^(-1)*f(D)*V if you have a diagonalizable matrix and D is a diagonal matrix containing the eigenvalues and V or its inverse contains the eigenvectors (I keep forgetting which)
I wouldn't keep your hopes up, though. The fact that your matrix is close to singular is a very bad starting point, generally speaking. Using singles on top of that doesn't help at all. But you'll only know by experimenting
So, if the matrix inverts under double but not single, your system is probably too stiff. If you don't need a very precise solution you can try using conjugate gradient (or similar) to do the inversion.
stackoverflow.com/questions/39714730/… here is one question which is up there from up about a week and also with a bounty but of mere 50 points. Thank you for looking at it.
@AnttiHaapala That one had a crappy pending edit. I fixed the edit so the question wouldn't go to the re-open queue when some robo-reviewer approved the pending edit.
As of Python 3.5, it supports matrix multiplication directly with the @ operator, versus C and Java, which implement these as library functions. Earlier versions of Python also used methods instead of an infix operator.[60][61]
@vaultah does it even make sense to talk about "atomic objects"?
@vaultah me neither, but it's usually preferred to back up edits with sources, otherwise it's an arbitrary change from one claim pulled out of one user's ass to another claim pulled out of another user's ass
Well, I would argue that the theories either helpfully describe a system behaviour, in which case the data is every time they are helpfully used and work. Or they don't. In which case they are disproven. Or there is no way to measure, in which case they live in a strange twilight world.
We were hoping that atomic-scale magnetic patterns might be safe, but then they developed spin-polarized scanning tunnelling microscopy, now they can measure magnetic moments on the atomic scale:(
@AndrasDeak That is pretty cool indeed (measuring atomic level magnetic moments). I don't understand the picture (except I assume that's an STM tip), which probably shows how far I am from understanding it!
@JRichardSnape b and c are the usual STM images, where the tip scans the surface and measure something that's close to the surface topology (where atoms sit)
it's a weird example of a skyrmion lattice on a surface with three-fold (120-degree-rotation) symmetry, with the skyrmion lattice having a four-fold (90-degree rotation) symmetry. Which is nuts.
(skyrmions are hedgehog-like localized spin patterns that happen to be topologically protected against local perturbations, if this says anything to you)
@JRichardSnape well eventually, sure. That's why anybody can get funding for research in magnetism:D
We already do that though, in HDD drives;)
but newer ideas include storing information in magnetic domain walls, or even skyrmions; these would all be MRAMs (magnetic RAMs)
I like the "hedehog like" I think that analogy makes sense to me. Comparing the spiral to the "hedgehog". Obviously, I am but an engineer of Very Small Brain, so there is no guarantee I'm interpreting them anywhere near correctly.
Sorry for my sloppy phrasing, I've been talking about magnetic skyrmions. Skyrmions generally are topologically protected meta-stable solutions of a field theoretical model. A kind of soliton, I think. Sorry for the confusion.
@khajvah it's the tiny ones that build up over time and end up breaking systems... as the major ones are normally caught as they're more obvious... good job though :)
Wiki says: "Often the term "jeans" refers to a particular style of pants, called "blue jeans," which were invented by Jacob W. Davis in partnership with Levi Strauss & Co. in 1871 and patented by Jacob W. Davis and Levi Strauss on May 20, 1873."