65           PetscScalar     Vmf, Vmt, thetaf, thetat, thetaft, thetatf;  in FormFunction_Subnet()  local
91           thetatf = thetat - thetaf;  in FormFunction_Subnet()
97 …] += Gtt * Vmt * Vmt + Vmt * Vmf * (Gtf * PetscCosScalar(thetatf) + Btf * PetscSinScalar(thetatf));  in FormFunction_Subnet()
98 …+= -Btt * Vmt * Vmt + Vmt * Vmf * (-Btf * PetscCosScalar(thetatf) + Gtf * PetscSinScalar(thetatf));  in FormFunction_Subnet()
230           PetscScalar     Vmf, Vmt, thetaf, thetat, thetaft, thetatf;  in FormJacobian_Subnet()  local
261           thetatf = thetat - thetaf;  in FormJacobian_Subnet()
303 …values[0] = Vmt * Vmf * (Gtf * -PetscSinScalar(thetatf) + Btf * PetscCosScalar(thetaft));         …  in FormJacobian_Subnet()
304 …alues[1] = 2.0 * Gtt * Vmt + Vmf * (Gtf * PetscCosScalar(thetatf) + Btf * PetscSinScalar(thetatf))…  in FormJacobian_Subnet()
305 …values[2] = Vmt * Vmf * (Gtf * PetscSinScalar(thetatf) + Btf * -PetscCosScalar(thetatf));         …  in FormJacobian_Subnet()
306 …values[3] = Vmt * (Gtf * PetscCosScalar(thetatf) + Btf * PetscSinScalar(thetatf));                …  in FormJacobian_Subnet()
317 …           values[0] = Vmt * Vmf * (Btf * PetscSinScalar(thetatf) + Gtf * PetscCosScalar(thetatf));  in FormJacobian_Subnet()
318 …ues[1] = -2.0 * Btt * Vmt + Vmf * (-Btf * PetscCosScalar(thetatf) + Gtf * PetscSinScalar(thetatf));  in FormJacobian_Subnet()
319 …         values[2] = Vmt * Vmf * (-Btf * PetscSinScalar(thetatf) + Gtf * -PetscCosScalar(thetatf));  in FormJacobian_Subnet()
320               values[3] = Vmt * (-Btf * PetscCosScalar(thetatf) + Gtf * PetscSinScalar(thetatf));  in FormJacobian_Subnet()