X-Git-Url: https://ruderich.org/simon/gitweb/?a=blobdiff_plain;f=maxima%2Fmaxima-init.mac;h=f270c10f3f5bef5153c78a901547447d4fc16c57;hb=8165a8ace54a29a17fa975278d83609c179d78de;hp=438989f1e567c4262931af244ee9e017065cd7ac;hpb=29daa783d21800160e6b262b932f984e67556e9e;p=config%2Fdotfiles.git

diff --git a/maxima/maxima-init.mac b/maxima/maxima-init.mac
index 438989f..f270c10 100644
--- a/maxima/maxima-init.mac
+++ b/maxima/maxima-init.mac
@@ -5,16 +5,17 @@
  */
 
 
-/* PACKAGES */
-
-/* Provides maxi() and mini() (and more). */
-load (descriptive);
+/* CUSTOM FUNCTIONS */
 
+/* Didn't find any euclid norm in Maxima, so here it is. */
+norm(x) := sqrt(transpose(x) . x);
 
-/* CUSTOM FUNCTIONS */
+/* Maxima's mat_norm() can't calculate the 2-norm. Thanks to Wolfgang Lindner
+ * (http://www.ma.utexas.edu/pipermail/maxima/2007/006300.html) for an
+ * improved version. */
+mat_norm2(x) := sqrt(lmax(eigenvalues(transpose(conjugate(x)).x)[1]));
 
-/* Maxima's norm() can't calculate the 2-norm. Only works for non-complex
- * matrices. */
-norm_2(x) := sqrt(maxi(eigenvalues(transpose(x) . x)[1]));
+/* Shortcut to perform lu factorization. */
+lu(x) := block(x : lu_factor(x), get_lu_factors(x));
 
 /* vim: ft=maxima */