INTEGER(kind = INT8), DIMENSION(0:5), PARAMETER                                 :: const_factorial_INT8 = (/ &
                                                                                    1_INT8, &
                                                                                    1_INT8, &
                                                                                    2_INT8, &
                                                                                    6_INT8, &
                                                                                    24_INT8, &
                                                                                    120_INT8  &
                                                                                /)
