Global csl::factorial (int n): The factorial function returns a double for now and so is quickly out of bounds for exact numerical precision. The limit is between 20! and 25!.
Global CSL_INF: Will not work for many cases, in particular log(0)!=. Not yet implemented.
File indicial.h: For now setting an object fully symmetric or antisymmetric will cause bug if all the indices don't lie in the same space. When taking permutations, indices from different spaces will be exchanged and the tensor will not accept them.