News for hackers. Stuff that matters.
Re:Special Linear Systems and Cholesky Factorization for Pro... by makmanalp on Tuesday September 19, 2017 @12:22PM (#15286021)
For linear algebra basics, I also especially recommend 3Blue1Brown's youtube channel and his "Essence of" series:

https://www.youtube.com/watch?v=kjBOesZCoqc&list=PLZHQObOWTQ...

I've never seen the material presented more clearly - it's completely visual and focused on building intuition first.

[ Parent | on HN ]
Re:Special Linear Systems and Cholesky Factorization for Pro... by dragandj on Tuesday September 19, 2017 @11:39AM (#15285599)
And the source of the library used: https://github.com/uncomplicate/neanderthal
[ Parent | on HN ]
Re:Special Linear Systems and Cholesky Factorization for Pro... by tw1010 on Tuesday September 19, 2017 @11:37AM (#15285576)
This is awesome. I can't wait until the future when we'll get more informal blog posts from programmers about even more abstract subjects which thus far are dominated in exposition by academics (who very rarely innovate in how they present the material). The day I see a good engineer-style post on HN about sheaf theory or algebraic geometry will be amazing.
[ Parent | on HN ]
Re:Special Linear Systems and Cholesky Factorization for Pro... by dragandj on Tuesday September 19, 2017 @11:41AM (#15285615)
Author here. I hope I won't spoil the future, since I work at a university, but I am also a programmer, and these libraries and tutorials have been made to be 100% practical with the working programmers as the target audience.
[ Parent | on HN ]
Re:Special Linear Systems and Cholesky Factorization for Pro... by mlevental on Tuesday September 19, 2017 @01:03PM (#15286468)
personally I'm waiting for an engineer-style post on IUT :)
[ Parent | on HN ]
Re:Special Linear Systems and Cholesky Factorization for Pro... by santaclaus on Tuesday September 19, 2017 @12:09PM (#15285881)
Now the real tough question is how do you pronounce Cholesky? I had a Polish math prof who was mad at the American pronunciation, but Cholesky was actually French, so should we actually pronounce it in the francophone style?
[ Parent | on HN ]
Re:Special Linear Systems and Cholesky Factorization for Pro... by mturmon on Tuesday September 19, 2017 @02:28PM (#15287432)
The sophisticated linear algebra practitioners I learned from all said "ko-LESS-key". People who ask me naive questions about the decomposition ("wait, is it upper triangular or lower triangular?") say "choe-LESS-key", or "CHOLES-key", or mutter something even farther off the mark.

Thus, I'm sticking with ko-LESS-key.

Also seeking opinions on: Weiner (as in the process), and Jensen (as in the inequality), Fourier (as in the series).

[ Parent | on HN ]
Re:Special Linear Systems and Cholesky Factorization for Pro... by dmichulke on Tuesday September 19, 2017 @04:45PM (#15288650)
It's a Wiener (you had a typo there) process and the name refers to a sausage or an inhabitants of the Austrian capital in German. So I'd pronounce it VEEner in German. Of course, if the guy was American then idk. His parents were German and Polish though.

Jensen was Danish so it's probably pronounced similar to German (?) so it'd be YENsen.

Fourier is Foo-r(i)EH in regular French.

[ Parent | on HN ]
Re:Special Linear Systems and Cholesky Factorization for Pro... by mindcrime on Tuesday September 19, 2017 @03:34PM (#15288052)
Fourier

I've always heard this pronounced like "for-e-ay" (or "for-e-eh" where the last part is the Canadian "eh?").

[ Parent | on HN ]
Re:Special Linear Systems and Cholesky Factorization for Pro... by mturmon on Tuesday September 19, 2017 @04:03PM (#15288296)
I'm with you. You also hear "FOUR-e-err", and my favorite EE teacher said "FUR-ee-err", indistinguishable from the word for a fur coat-maker, which always produced a smile. (But, he's a fellow of the IEEE and I'm not, so who's laughing now?)

I say "WEE-ner" (since he was American, but, open to correction on that one) and "YEN-sen" (that's how I was taught).

[ Parent | on HN ]
Re:Special Linear Systems and Cholesky Factorization for Pro... by santaclaus on Tuesday September 19, 2017 @04:14PM (#15288382)
> FOUR-e-err

My favorites are 'you-ler' for Euler and 'goose-e-in' for Gaussian.

[ Parent | on HN ]
Re:Special Linear Systems and Cholesky Factorization for Pro... by llamaz on Wednesday September 20, 2017 @12:59AM (#15290963)
I've always heard Wiener pronounced wine-r, but that might just be to avoid saying ween-er
[ Parent | on HN ]
Re:Special Linear Systems and Cholesky Factorization for Pro... by jules on Tuesday September 19, 2017 @01:39PM (#15286831)
One can avoid this problem by calling it the LL^T-decomposition.
[ Parent | on HN ]
Re:Special Linear Systems and Cholesky Factorization for Pro... by CoreXtreme on Tuesday September 19, 2017 @02:29PM (#15287445)
Font size at 23px. Seriously?

Just skimmed and now my eyes are hurt.

Computation cost of calculating the Cholesky factorization of the original second derivative matrix is comparable to cost of calculating the inverse factorization.

[ Parent | on HN ]
Re:Special Linear Systems and Cholesky Factorization for Pro... by amelius on Tuesday September 19, 2017 @11:59AM (#15285774)
Since LAPACK is an excellent library for these matrix functions, why not use it, and translate it to Clojure instead of reinventing the wheel? (There's also a translator from Fortran to C).
[ Parent | on HN ]
Re:Special Linear Systems and Cholesky Factorization for Pro... by dragandj on Tuesday September 19, 2017 @12:05PM (#15285834)
Because JVM does not have access to the hardware features needed to execute this efficiently. Also, the power of lapack comes from existing highly tuned implementations, not the interface itself. Reference BLAS and LAPACK are actually quite slow.
[ Parent | on HN ]
Re:Special Linear Systems and Cholesky Factorization for Pro... by TurpIF on Tuesday September 19, 2017 @02:12PM (#15287249)
In fact, you can access the hardware through native calls via JNI (or JNA). Or course, then, you have to embedded multiplatform libraries and manage the associated issues. Also, the OpenBLAS implementation is very well optimized for several Intel and AMD processors (you can compile it so that it autodetects which one you're using). It can even reach the efficiency of the Intel's MKL implementation in mono-threaded mode.
[ Parent | on HN ]
Re:Special Linear Systems and Cholesky Factorization for Pro... by dragandj on Tuesday September 19, 2017 @02:27PM (#15287415)
We don't even have to guess, since that's exactly what Neanderthal does. Also, I micro-benchmarked lots of options and have yet to find one that fills similar use case that is faster than Neanderhtal+MKL on the CPU, regardless of the JNI overhead (minus the obvious direct use of MKL, but that is much more low-level code). Also, most higher level libraries have considerable overhead. Neanderthal's overhead is tiny.

OpenBLAS's huge drawback is that it only supports BLAS without LAPACK, sparse, tensors, FFT etc.

Anyway, regarding the OP's comment, I guess that they meant to suggest implementing all that in pure Java, not Java + FFI, since then the native code has to be written in non-Java.

[ Parent | on HN ]
Re:Special Linear Systems and Cholesky Factorization for Pro... by reikonomusha on Tuesday September 19, 2017 @02:59PM (#15287764)
What about correctness? Numerical linear algebra is very difficult and bugs are hard to test for and weed out.
[ Parent | on HN ]
Re:Special Linear Systems and Cholesky Factorization for Pro... by dragandj on Tuesday September 19, 2017 @03:06PM (#15287826)
I guess when it comes to testing, Intel has more resources to do that than most other folks. That's one of the reasons why I use MKL despite it not being open source (but it is free as in free beer).
[ Parent | on HN ]