You need to either find a implementation in theano, or wrap the scipy function in theano using @as_op.