I have an algorithm for calibrating a vector magnetometer. The input is N readings of the x, y, z axes: , , and .
The algorithm fits an ellipsoid to the data by estimating a symmetric matrix A. In order to calibrate the system, it needs to calculate . I am adapting the algorithm for a microcontroller with very little memory, so cannot load standard matrix manipulation libraries.
Is there an explicit formula for calculating the square root of positive definite matrix?