Как можно добавить в цикл вычисление определителя матрицы?

Подскажите пожалуйста , как можно запихнуть в этот код вычисление определителя, т.е. до того как происходит деление на temp и в тот момент когда матрица приведена к "почти единичному" виду (когда на диагонали стоят не единицы, а случайные числа)

Вычисление определителя:

double det=1;
for(int k =0;k<size;k++) det*=A[k][k];

Код:

for (int k = 0; k < size; k++)
    {
        temp = A[k][k];

        for (int j = 0; j < size; j++)
        {
            A[k][j] /= temp;
            T[k][j] /= temp;
        }

        for (int i = k + 1; i < size; i++)
        {
            temp = A[i][k];

            for (int j = 0; j < size; j++)
            {
                A[i][j] -= A[k][j] * temp;
                T[i][j] -= T[k][j] * temp;
            }
        }
    }

    for (int k = size - 1; k > 0; k--)
    {
        for (int i = k - 1; i >= 0; i--)
        {
            temp = A[i][k];

            for (int j = 0; j < size; j++)
            {
                A[i][j] -= A[k][j] * temp;
                T[i][j] -= T[k][j] * temp;
            }
        }
    }
```