fork download 
  1. import numpy
  2.  
  3. a = numpy.random.rand(100, 30)*10
  4. b = numpy.random.rand(30, 50)
  5.  
  6. out1 = numpy.exp(numpy.matmul(numpy.log(a), b))
  7. out = numpy.ones((100, 50))
  8. for i in range(100):
  9.     for j in range(50):
  10.         for k in range(30):
  11.             out[i, j] *= a[i, k]**b[k, j]
  12.  
  13. print(numpy.abs(out1-out).max())
Success #stdin #stdout 0.19s 93248KB
comments (?)
stdin 
Standard input is empty
stdout 
2.25
https://ideone.com/nQKvjj
language:
Python 3 (python  3.12)
created:
8 years ago
visibility:
public

Share or Embed source code

Discover > Sphere Engine API

The brand new service which powers Ideone!

Discover > IDE Widget

Widget for compiling and running the source code in a web browser!