initialize tensor module
subroutine initialize_tensors()
    call compute_tensor_coefficients()
    call compute_factorials()
    call compute_double_factorials()
    call compute_binomial_coefficients()
    call compute_trinomial_coefficients()
end subroutine initialize_tensors