Question: The dchange command returns zero

Hello Maple primers

I am trying to do a coordinate transformation which involves a number of partial derivatives. I turned to the PDETools[dchange] to accomplish this but it will just return 0 when run through.

The problem is as such. First two functions are defined which contain a lot of "stuff". Then the forward and reverse transformations are defined between the coordinates; and finally the transformation is done.

   del_1:=Diff(Y,r1,r1)+2/r1*(Diff(Y,r1)) + Diff(Y,r2,r2)+2/r2*(Diff(Y,r2))+2*(Diff(Y,r3,r3))+4/r3*(Diff(Y,r3))+((r1^2-r2^2+r3^2)/(r1*r3))*(Diff(Y,r1,r3))+(r2^2-r1^2+r3^2)/(r2*r3)*(Diff(Y,r2,r3)):

   del_2:=Diff(Y,r1,r1)+2/r1*(Diff(Y,r1))+(r1^2-r2^2+r3^2)/(r1*r3)*Diff(Y,r1,r3)+Diff(Y,r2,r2)+2/r2*(Diff(Y,r2))+(-r1^2+r2^2+r3^2)/(r2*r3)*(Diff(Y,r2,r3))+4/r3*(Diff(Y,r3))+2*(Diff(Y,r3,r3)):

#Forward transformation

R1:=0.5*(v)+0.25*(w):
R2:=0.5*(u)+0.25*(w):
R3:=0.5*(u)+0.5*(v):

tr:={r1=R1,r2=R2,r3=R3}:

#Reverse transformation
rR1:=(-r1+r2+r3):
rR2:=(r1-r2+r3):
rR3:=2*(r1+r2-r3):

rtr:={u=rR1,v=rR2,w=rR3}:
nv:={u,v,w}:
AA:=simplify(del_1+del_2);
PDEtools[dchange](tr,AA,nv,rtr);

This will return zero. Is there something obvious I am missing here? I have used the dchange tool before in a similar manner and it has worked without issue.

Please Wait...