 		}
 	}
 
-	for ( i = 0; i < M->cols; i++ ) {
-		elem(a, 0, i) /= M->rows;
-	}
+	m_elem_mult(a, 1.0f / M->rows);
 
 	return a;
 }


	// print debug information
	if ( LOGGER(LL_DEBUG) ) {
		printf("%s [%d,%d], %s [%d,%d] <- eig(%s [%d,%d])\n",
		       "V", V->rows, V->cols,
		       "D", D->rows, D->cols,
		       "M", M->rows, M->cols);

		}
	}

	m_elem_mult(a, 1.0f / M->cols);

	return a;
}
