Hab mal nen Moment nachgedacht und glaube es geht alles auch sehr viel schneller:
det ( A-1 * D-1 * A )
= 1/ det(a) * 1/ det(D) * det (A) = 1/det(D) = 1 / -2a
Spur geht zwar nicht faktorenweise aber es ist Spur(X*Y) = Spur(Y*X)
also auch Spur ( A-1 * D-1 * A) = Spur ( A-1 * A * D-1 )
= Spur ( E * D-1 )
= Spur ( D-1) = 1 -1/2 + 1/a
Damit bräuchte man A-1 noch nicht mal auszurechnen, und D-1 ist ja nun wirklich
schnell hingeschrieben.