¿La forma más rápida de hacer esta doble suma?

¿Cuál es la forma más rápida de hacer este resumen?

Esto es lo que tengo hasta ahora

ans = 0
for (i in 1:dimx[1]){
  for (j in 1:dimx[2]){
    ans = ans + ((x[i,j] - parameters$mu)^2)/(parameters$omega_2[i]*parameters$sigma_2[j])
  }
}

dóndeomega_2ysigma_2 son omega ^ 2 y sigma ^ 2 respectivamente.