**Question 1** What is the number of pairs of commuting elements in GL_n(F_q) ?

I am aware of many results concerning commuting elements in Mat_n(F_q), but I am interested in GL i.e. non-degenerate matrices. In particular Feit, Fine 1960, results for many kinds of Lie algebras obtained recently Jason Fulman, Robert Guralnick, and even Motivic Donaldson... is related to counting nilpotent matrices - see beautiful answer on math.se by Olivier Schiffmann.

**[EDIT after answers]**
Answers say: the number of commuting pairs in ANY group equals to number of conjugacy classes multiply order of group. Number of conjugacy classes in $GL(n,F_q)$ has nice generating function
MO8415, MO104457 (or links in answers below):
$$ \sum_{n\geq 0} |number~of~conjugacy~classes~in~GL(n,F_q)| x^n = \prod_{j\geq 1}
\frac{1-x^j}{1-qx^j} $$

So the question is done.
**[end EDIT]**

**Question 1b** If we take limit q->1 can we get number of commuting elements in S_n ? (Standard analogy S_n = GL_n(F_1) ).

**[EDIT after answers]**
There is beautiful formula for $S_n$ which can count m-tuples of commuting elements, number of involutions, number of elements of order k and so on - *just one formula can do it all* - see MO272045.
From the answer there we know that there is certain q-analog of that formula for GL(n,F_q). However it does not immediately cover the case of commuting pairs (m-tuples)
and also the limit q->1 is not clear. In particular it is not clear to me
how generating function above can in the limit q->1 give generating functions for conjugacy classes in $S_n$ - partition function. So there is place to think more (imho).
**[end EDIT]**

**Question 2** what about triples, m-tuples ?

**[EDIT after answers]**
Is it open problem ? (Since it is not covered by powerful analysis in papers quoted in MO answer. It suggest that should be nice generating function):

$$ \sum_{n\geq 0} \frac{|commuting~~(m-tuples)~in~GL(n,F_q)|}{|GL(n,F_q)|} x^n = ??? $$

For commuting pairs we will get generating function above, since order of $GL(n,F_q)$ in numerator cancel one in denominator.
**[end EDIT]**

**Question 3** what about other finite algebraic groups ? (In view of results
Fulman, Guralnick for Lie algebras it seems natural).

**Question 4** for compact simple groups over R one can put for example Killing Riemnanian metric and consider volumes for such manifolds - are there some results ?

PS By the way nice fact on pairs of commuting elements is here: 5/8 bound in group theory

PSPS

**Question 1c**
By the way can computer algebra systems like GAP or MAGMA get the answer on such question for small "n", but in the form of polynomial in "q",
or they can deal only with fixed q ?

(See answer in comment by Alexander Konovalov ).