Suppose I have two representations (over $\mathbb{C}$) of a finite group $G$, $\rho : G \to GL(V)$ and $\tau : G \to GL(W)$. If I am given that these representations are isomorphic, i.e. there exists a linear map $A : V \to W$ such that for all $g \in G$, $\tau(g)A = A\rho(g)$, how can I find $A$?
The approach I have in mind is to compute $\tau(g)$ and $\rho(g)$ for all $g \in G$, then find the matrix that exhibits the simultaneous similarity of these. But this is just kind of rephrasing the problem.
One way of finding a solution would be to make each entry of $A$ a variable, then write out the linear equations given by $\tau(g)A = A\rho(g)$ for all $g \in G$. Then the solution space will have some invertible matrices in it, these are the ones I want. But I have no idea how to program this on a computer.
Another rephrasing of that idea is finding simultaneous fixed points of the linear maps given by:
$$A \mapsto \tau(g)A\rho(g^{-1}) $$
for each $g \in G$
This is going to be a very large matrix. And if I "unwrap" the matrix $A$ into vectors, the matrix of this linear map will be huge: if the representations have degree $n$, I'll end up with a matrix of size $n^2 \times n^2$ for which I have to find the eigenspace for value $1$.
Surely there is a smarter way of solving the problem than this?
I am using the GAP algebra system, so answers relating to this would be the most helpful, but any help is appreciated.
Some more information: usually I am dealing with the case where $G$ is a fairly large (tens of thousands of elements) permutation group, but I know the generators of $G$ and there are usually only a handful.
looking naively at the matter, the matrix $A$ has $n^2$ entries, which are initially unknown. and you have a set of $\#(G)$ linear equations in the $A_{ij}$:
$$ A_{ij} = \sum_{i} \sum_{j}\tau(g)_{ik}A_{kl}\rho(g^{-1})_{lj} $$