Question: Long waiting time of plotting the chen–gackstatter surface

Hello do u know how to make this code more optimal? this is very long time to plot.

g := evalf((Beta(1/4, 1/4)/2)^4);

A := Pi*r*cos(theta) - Re(WeierstrassZeta(r*exp(theta*I), g, 0) + Pi*WeierstrassPPrime(r*exp(theta*I), g, 0)/g);
B := Pi*r*sin(theta) + Im(WeierstrassZeta(r*exp(theta*I), g, 0) - Pi*WeierstrassPPrime(r*exp(theta*I), g, 0)/g);
C := sqrt(6*Pi/g)*Re(WeierstrassP(r*exp(theta*I), g, 0));

plot3d([A, B, C], r = 1/5 .. 4/5, theta = -Pi .. Pi, view = [-8 .. 8, -8 .. 8, -8 .. 8], shading = zhue, grid = [200, 200]);

Please Wait...