📄 libearingfilm.for
字号:
P(JX)=P(JX)+BXY(M)*RTY(I,N,M)+BXZ(M)*RTZ(I,N,M)-AXX(M)*RUX(I,S,M)
P(JY)=P(JY)+BYY(M)*RTY(I,N,M)+BYZ(M)*RTZ(I,N,M)-AYX(M)*RUX(I,S,M)
P(JZ)=P(JZ)+BZY(M)*RTY(I,N,M)+BZZ(M)*RTZ(I,N,M)-AZX(M)*RUX(I,S,M)
GOTO 450
378 A(JX,KX)=A(JX,KX)+AXX(M)
A(JX,KY)=A(JX,KY)-BXY(M)
A(JX,KZ)=A(JX,KZ)-BXZ(M)
A(JY,KX)=A(JY,KX)+AYX(M)
A(JY,KY)=A(JY,KY)-BYY(M)
A(JY,KZ)=A(JY,KZ)-BYZ(M)
A(JZ,KX)=A(JZ,KX)+AZX(M)
A(JZ,KY)=A(JZ,KY)-BZY(M)
A(JZ,KZ)=A(JZ,KZ)-BZZ(M)
P(JX)=P(JX)+BXX(M)*RTX(I,N,M)-AXY(M)*RUY(I,S,M)-AXZ(M)*RUZ(I,S,M)
P(JY)=P(JY)+BYX(M)*RTX(I,N,M)-AYY(M)*RUY(I,S,M)-AYZ(M)*RUZ(I,S,M)
P(JZ)=P(JZ)+BZX(M)*RTX(I,N,M)-AZY(M)*RUY(I,S,M)-AZZ(M)*RUZ(I,S,M)
GOTO 450
380 A(JX,KX)=A(JX,KX)+AXX(M)
A(JX,KY)=A(JX,KY)-BXY(M)
A(JX,KZ)=A(JX,KZ)+AXZ(M)
A(JY,KX)=A(JY,KX)+AYX(M)
A(JY,KY)=A(JY,KY)-BYY(M)
A(JY,KZ)=A(JY,KZ)+AYZ(M)
A(JZ,KX)=A(JZ,KX)+AZX(M)
A(JZ,KY)=A(JZ,KY)-BZY(M)
A(JZ,KZ)=A(JZ,KZ)+AZZ(M)
P(JX)=P(JX)+BXX(M)*RTX(I,N,M)+BXZ(M)*RTZ(I,N,M)-AXY(M)*RUY(I,S,M)
P(JY)=P(JY)+BYX(M)*RTX(I,N,M)+BYZ(M)*RTZ(I,N,M)-AYY(M)*RUY(I,S,M)
P(JZ)=P(JZ)+BZX(M)*RTX(I,N,M)+BZZ(M)*RTZ(I,N,M)-AZY(M)*RUY(I,S,M)
GOTO 450
382 A(JX,KX)=A(JX,KX)-BXX(M)
A(JX,KY)=A(JX,KY)+AXY(M)
A(JX,KZ)=A(JX,KZ)-BXZ(M)
A(JY,KX)=A(JY,KX)-BYX(M)
A(JY,KY)=A(JY,KY)+AYY(M)
A(JY,KZ)=A(JY,KZ)-BYZ(M)
A(JZ,KX)=A(JZ,KX)-BZX(M)
A(JZ,KY)=A(JZ,KY)+AZY(M)
A(JZ,KZ)=A(JZ,KZ)-BZZ(M)
P(JX)=P(JX)+BXY(M)*RTY(I,N,M)-AXX(M)*RUX(I,S,M)-AXZ(M)*RUZ(I,S,M)
P(JY)=P(JY)+BYY(M)*RTY(I,N,M)-AYX(M)*RUX(I,S,M)-AYZ(M)*RUZ(I,S,M)
P(JZ)=P(JZ)+BZY(M)*RTY(I,N,M)-AZX(M)*RUX(I,S,M)-AZZ(M)*RUZ(I,S,M)
GOTO 450
385 A(JX,KY)=A(JX,KY)+AXY(M)
A(JX,KZ)=A(JX,KZ)+AXZ(M)
A(JY,KY)=A(JY,KY)+AYY(M)
A(JY,KZ)=A(JY,KZ)+AYZ(M)
A(JY,KY)=A(JY,KY)-BYY(M)
A(JY,KY)=A(JY,KY)-BYY(M)
P(JX)=P(JX)-AXX(M)*RUX(I,S,M)
P(JY)=P(JY)-AYX(M)*RUX(I,S,M)
P(JZ)=P(JZ)-AZX(M)*RUX(I,S,M)
GOTO 450
388 A(JX,KX)=A(JX,KX)+AXX(M)
A(JX,KZ)=A(JX,KZ)+AXZ(M)
A(JY,KX)=A(JY,KX)+AYX(M)
A(JY,KZ)=A(JY,KZ)+AYZ(M)
A(JZ,KX)=A(JZ,KX)+AZX(M)
A(JZ,KZ)=A(JZ,KZ)+AZZ(M)
P(JX)=P(JX)+BXX(M)*RTX(I,N,M)+BXY(M)*RTY(I,N,M)+BXZ(M)*RTZ(I,N,M)
P(JY)=P(JY)+BYX(M)*RTX(I,N,M)+BYY(M)*RTY(I,N,M)+BYZ(M)*RTZ(I,N,M)
P(JZ)=P(JZ)+BZX(M)*RTX(I,N,M)+BZY(M)*RTY(I,N,M)+BZZ(M)*RTZ(I,N,M)
P(JX)=P(JX)-AXY(M)*RUY(I,S,M)
P(JY)=P(JY)-AYY(M)*RUY(I,S,M)
P(JZ)=P(JZ)-AZY(M)*RUY(I,S,M)
GOTO 450
390 A(JX,KX)=A(JX,KX)+AXX(M)
A(JX,KY)=A(JX,KY)+AXY(M)
A(JY,KX)=A(JY,KX)+AYX(M)
A(JY,KY)=A(JY,KY)+AYY(M)
A(JZ,KX)=A(JZ,KX)+AZX(M)
A(JZ,KY)=A(JZ,KY)+AZY(M)
P(JX)=P(JX)+BXX(M)*RTX(I,N,M)+BXY(M)*RTY(I,N,M)+BXZ(M)*RTZ(I,N,M)
P(JY)=P(JY)+BYX(M)*RTX(I,N,M)+BYY(M)*RTY(I,N,M)+BYZ(M)*RTZ(I,N,M)
P(JZ)=P(JZ)+BZX(M)*RTX(I,N,M)+BZY(M)*RTY(I,N,M)+BZZ(M)*RTZ(I,N,M)
P(JX)=P(JX)-AXZ(M)*RUZ(I,S,M)
P(JY)=P(JY)-AYZ(M)*RUZ(I,S,M)
P(JZ)=P(JZ)-AZZ(M)*RUZ(I,S,M)
GOTO 450
392 A(JX,KY)=A(JX,KY)-BXY(M)
A(JX,KZ)=A(JX,KZ)-BXZ(M)
A(JY,KY)=A(JY,KY)-BYY(M)
A(JY,KZ)=A(JY,KZ)-BYZ(M)
A(JZ,KY)=A(JZ,KY)-BZY(M)
A(JZ,KZ)=A(JZ,KZ)-BZZ(M)
P(JX)=P(JX)+BXX(M)*RTX(I,N,M)
P(JY)=P(JY)+BYX(M)*RTX(I,N,M)
P(JZ)=P(JZ)+BZX(M)*RTX(I,N,M)
P(JX)=P(JX)-AXX(M)*RUX(I,S,M)-AXY(M)*RUY(I,S,M)-AXZ(M)*RUZ(I,S,M)
P(JY)=P(JY)-AYX(M)*RUX(I,S,M)-AYY(M)*RUY(I,S,M)-AYZ(M)*RUZ(I,S,M)
P(JZ)=P(JZ)-AZX(M)*RUX(I,S,M)-AZY(M)*RUY(I,S,M)-AZZ(M)*RUZ(I,S,M)
GOTO 450
395 A(JX,KX)=A(JX,KX)-BXX(M)
A(JX,KZ)=A(JX,KZ)-BXZ(M)
A(JY,KX)=A(JY,KX)-BYX(M)
A(JY,KZ)=A(JY,KZ)-BYZ(M)
A(JZ,KX)=A(JZ,KX)-BZX(M)
A(JZ,KZ)=A(JZ,KZ)-BZZ(M)
P(JX)=P(JX)+BXY(M)*RTY(I,N,M)
P(JY)=P(JY)+BYY(M)*RTY(I,N,M)
P(JZ)=P(JZ)+BZY(M)*RTY(I,N,M)
P(JX)=P(JX)-AXX(M)*RUX(I,S,M)-AXY(M)*RUY(I,S,M)-AXZ(M)*RUZ(I,S,M)
P(JY)=P(JY)-AYX(M)*RUX(I,S,M)-AYY(M)*RUY(I,S,M)-AYZ(M)*RUZ(I,S,M)
P(JZ)=P(JZ)-AZX(M)*RUX(I,S,M)-AZY(M)*RUY(I,S,M)-AZZ(M)*RUZ(I,S,M)
GOTO 450
398 A(JX,KX)=A(JX,KX)-BXX(M)
A(JX,KY)=A(JX,KY)-BXY(M)
A(JY,KX)=A(JY,KX)-BYX(M)
A(JY,KY)=A(JY,KY)-BYY(M)
A(JZ,KX)=A(JZ,KX)-BZX(M)
A(JZ,KY)=A(JZ,KY)-BZY(M)
P(JX)=P(JX)+BXZ(M)*RTZ(I,N,M)
P(JY)=P(JY)+BYZ(M)*RTZ(I,N,M)
P(JZ)=P(JZ)+BZZ(M)*RTZ(I,N,M)
P(JX)=P(JX)-AXX(M)*RUX(I,S,M)-AXY(M)*RUY(I,S,M)-AXZ(M)*RUZ(I,S,M)
P(JY)=P(JY)-AYX(M)*RUX(I,S,M)-AYY(M)*RUY(I,S,M)-AYZ(M)*RUZ(I,S,M)
P(JZ)=P(JZ)-AZX(M)*RUX(I,S,M)-AZY(M)*RUY(I,S,M)-AZZ(M)*RUZ(I,S,M)
GOTO 450
400 A(JX,KX)=A(JX,KX)+AXX(M)
A(JX,KZ)=A(JX,KZ)+AXZ(M)
A(JY,KX)=A(JY,KX)+AYX(M)
A(JY,KZ)=A(JY,KZ)+AYZ(M)
A(JY,KX)=A(JY,KX)-BYX(M)
A(JY,KZ)=A(JY,KZ)-BYZ(M)
A(JZ,KX)=A(JZ,KX)-BZX(M)
A(JZ,KZ)=A(JZ,KZ)-BZZ(M)
P(JX)=P(JX)+BXX(M)*RTX(I,N,M)+BXY(M)*RTY(I,N,M)
P(JY)=P(JY)+BYX(M)*RTX(I,N,M)+BYY(M)*RTY(I,N,M)
P(JZ)=P(JZ)+BZX(M)*RTX(I,N,M)+BZY(M)*RTY(I,N,M)
P(JX)=P(JX)-AXY(M)*RUY(I,S,M)-AXZ(M)*RUZ(I,S,M)
P(JY)=P(JY)-AYY(M)*RUY(I,S,M)-AYZ(M)*RUZ(I,S,M)
P(JZ)=P(JZ)-AZY(M)*RUY(I,S,M)-AZZ(M)*RUZ(I,S,M)
GOTO 450
402 A(JX,KX)=A(JX,KX)+AXX(M)
A(JX,KY)=A(JX,KY)+AXY(M)
A(JZ,KX)=A(JZ,KX)+AZX(M)
A(JZ,KY)=A(JZ,KY)+AZY(M)
A(JY,KX)=A(JY,KX)-BYX(M)
A(JY,KY)=A(JY,KY)-BYY(M)
A(JZ,KX)=A(JZ,KX)-BZX(M)
A(JZ,KY)=A(JZ,KY)-BZY(M)
P(JX)=P(JX)+BXX(M)*RTX(I,N,M)+BXZ(M)*RTZ(I,N,M)
P(JY)=P(JY)+BYX(M)*RTX(I,N,M)+BYZ(M)*RTZ(I,N,M)
P(JZ)=P(JZ)+BZX(M)*RTX(I,N,M)+BZZ(M)*RTZ(I,N,M)
P(JX)=P(JX)-AXY(M)*RUY(I,S,M)-AXZ(M)*RUZ(I,S,M)
P(JY)=P(JY)-AYY(M)*RUY(I,S,M)-AYZ(M)*RUZ(I,S,M)
P(JZ)=P(JZ)-AZY(M)*RUY(I,S,M)-AZZ(M)*RUZ(I,S,M)
GOTO 450
405 A(JY,KX)=A(JY,KX)+AYX(M)
A(JY,KY)=A(JY,KY)+AYY(M)
A(JZ,KX)=A(JZ,KX)+AZX(M)
A(JZ,KY)=A(JZ,KY)+AZY(M)
A(JX,KX)=A(JX,KX)-BXX(M)
A(JX,KY)=A(JX,KY)-BXY(M)
A(JZ,KX)=A(JZ,KX)-BZX(M)
A(JZ,KY)=A(JZ,KY)-BZY(M)
P(JX)=P(JX)+BXY(M)*RTY(I,N,M)+BXZ(M)*RTZ(I,N,M)
P(JY)=P(JY)+BYY(M)*RTY(I,N,M)+BYZ(M)*RTZ(I,N,M)
P(JZ)=P(JZ)+BZY(M)*RTY(I,N,M)+BZZ(M)*RTZ(I,N,M)
P(JX)=P(JX)-AXX(M)*RUX(I,S,M)-AXZ(M)*RUZ(I,S,M)
P(JY)=P(JY)-AYX(M)*RUX(I,S,M)-AYZ(M)*RUZ(I,S,M)
P(JZ)=P(JZ)-AZX(M)*RUX(I,S,M)-AZZ(M)*RUZ(I,S,M)
GOTO 450
408 A(JX,KY)=A(JX,KY)+AXY(M)
A(JX,KZ)=A(JX,KZ)+AXZ(M)
A(JY,KY)=A(JY,KY)+AYY(M)
A(JY,KZ)=A(JY,KZ)+AYZ(M)
A(JX,KY)=A(JX,KY)-BXY(M)
A(JX,KZ)=A(JX,KZ)-BXZ(M)
A(JZ,KY)=A(JZ,KY)-BZY(M)
A(JZ,KZ)=A(JZ,KZ)-BZZ(M)
P(JX)=P(JX)+BXX(M)*RTX(I,N,M)+BXY(M)*RTY(I,N,M)
P(JY)=P(JY)+BYX(M)*RTX(I,N,M)+BYY(M)*RTY(I,N,M)
P(JZ)=P(JZ)+BZX(M)*RTX(I,N,M)+BZY(M)*RTY(I,N,M)
P(JX)=P(JX)-AXX(M)*RUX(I,S,M)-AXZ(M)*RUZ(I,S,M)
P(JY)=P(JY)-AYX(M)*RUX(I,S,M)-AYZ(M)*RUZ(I,S,M)
P(JZ)=P(JZ)-AZX(M)*RUX(I,S,M)-AZZ(M)*RUZ(I,S,M)
GOTO 450
410 A(JX,KY)=A(JX,KY)+AXY(M)
A(JX,KZ)=A(JX,KZ)+AXZ(M)
A(JZ,KY)=A(JZ,KY)+AZY(M)
A(JZ,KZ)=A(JZ,KZ)+AZZ(M)
A(JX,KY)=A(JX,KY)-BXY(M)
A(JX,KZ)=A(JX,KZ)-BXZ(M)
A(JY,KY)=A(JY,KY)-BYY(M)
A(JY,KZ)=A(JY,KZ)-BYZ(M)
P(JX)=P(JX)+BXX(M)*RTX(I,N,M)+BXZ(M)*RTZ(I,N,M)
P(JY)=P(JY)+BYX(M)*RTX(I,N,M)+BYZ(M)*RTZ(I,N,M)
P(JZ)=P(JZ)+BZX(M)*RTX(I,N,M)+BZZ(M)*RTZ(I,N,M)
P(JX)=P(JX)-AXX(M)*RUX(I,S,M)-AXY(M)*RUY(I,S,M)
P(JY)=P(JY)-AYX(M)*RUX(I,S,M)-AYY(M)*RUY(I,S,M)
P(JZ)=P(JZ)-AZX(M)*RUX(I,S,M)-AZY(M)*RUY(I,S,M)
GOTO 450
412 A(JY,KX)=A(JY,KX)+AYX(M)
A(JY,KZ)=A(JY,KZ)+AYZ(M)
A(JZ,KX)=A(JZ,KX)+AZX(M)
A(JZ,KZ)=A(JZ,KZ)+AZZ(M)
P(JX)=P(JX)+BXY(M)*RTY(I,N,M)+BXZ(M)*RTZ(I,N,M)
P(JY)=P(JY)+BYY(M)*RTY(I,N,M)+BYZ(M)*RTZ(I,N,M)
P(JZ)=P(JZ)+BZY(M)*RTY(I,N,M)+BZZ(M)*RTZ(I,N,M)
A(JX,KX)=A(JX,KX)-BXX(M)
A(JX,KZ)=A(JX,KZ)-BXZ(M)
A(JY,KX)=A(JY,KX)-BYX(M)
A(JY,KZ)=A(JY,KZ)-BYZ(M)
P(JX)=P(JX)-AXX(M)*RUX(I,S,M)-AXY(M)*RUY(I,S,M)
P(JY)=P(JY)-AYX(M)*RUX(I,S,M)-AYY(M)*RUY(I,S,M)
P(JZ)=P(JZ)-AZX(M)*RUX(I,S,M)-AZY(M)*RUY(I,S,M)
GOTO 450
415 A(JX,KZ)=A(JX,KZ)+AXZ(M)
A(JY,KZ)=A(JY,KZ)+AYZ(M)
A(JZ,KZ)=A(JZ,KZ)+AZZ(M)
P(JX)=P(JX)+BXX(M)*RTX(I,N,M)+BXY(M)*RTY(I,N,M)+BXZ(M)*RTZ(I,N,M)
P(JY)=P(JY)+BYX(M)*RTX(I,N,M)+BYY(M)*RTY(I,N,M)+BYZ(M)*RTZ(I,N,M)
P(JZ)=P(JZ)+BZX(M)*RTX(I,N,M)+BZY(M)*RTY(I,N,M)+BZZ(M)*RTZ(I,N,M)
P(JX)=P(JX)-AXX(M)*RUX(I,S,M)-AXY(M)*RUY(I,S,M)
P(JY)=P(JY)-AYX(M)*RUX(I,S,M)-AYY(M)*RUY(I,S,M)
P(JZ)=P(JZ)-AZX(M)*RUX(I,S,M)-AZY(M)*RUY(I,S,M)
GOTO 450
418 A(JX,KX)=A(JX,KX)+AXX(M)
A(JY,KX)=A(JY,KX)+AYX(M)
A(JZ,KX)=A(JZ,KX)+AZX(M)
P(JX)=P(JX)+BXX(M)*RTX(I,N,M)+BXY(M)*RTY(I,N,M)+BXZ(M)*RTZ(I,N,M)
P(JY)=P(JY)+BYX(M)*RTX(I,N,M)+BYY(M)*RTY(I,N,M)+BYZ(M)*RTZ(I,N,M)
P(JZ)=P(JZ)+BZX(M)*RTX(I,N,M)+BZY(M)*RTY(I,N,M)+BZZ(M)*RTZ(I,N,M)
P(JX)=P(JX)-AXY(M)*RUY(I,S,M)-AXZ(M)*RUZ(I,S,M)
P(JY)=P(JY)-AYY(M)*RUY(I,S,M)-AYZ(M)*RUZ(I,S,M)
P(JZ)=P(JZ)-AZY(M)*RUY(I,S,M)-AZZ(M)*RUZ(I,S,M)
GOTO 450
420 A(JX,KY)=A(JX,KY)+AXY(M)
A(JY,KY)=A(JY,KY)+AYY(M)
A(JZ,KY)=A(JZ,KY)+AZY(M)
P(JX)=P(JX)+BXX(M)*RTX(I,N,M)+BXY(M)*RTY(I,N,M)+BXZ(M)*RTZ(I,N,M)
P(JY)=P(JY)+BYX(M)*RTX(I,N,M)+BYY(M)*RTY(I,N,M)+BYZ(M)*RTZ(I,N,M)
P(JZ)=P(JZ)+BZX(M)*RTX(I,N,M)+BZY(M)*RTY(I,N,M)+BZZ(M)*RTZ(I,N,M)
P(JX)=P(JX)-AXX(M)*RUX(I,S,M)-AXZ(M)*RUZ(I,S,M)
P(JY)=P(JY)-AYX(M)*RUX(I,S,M)-AYZ(M)*RUZ(I,S,M)
P(JZ)=P(JZ)-AZX(M)*RUX(I,S,M)-AZZ(M)*RUZ(I,S,M)
GOTO 450
422 A(JX,KZ)=A(JX,KZ)-BXZ(M)
A(JY,KZ)=A(JY,KZ)-BYZ(M)
A(JZ,KZ)=A(JZ,KZ)-BZZ(M)
P(JX)=P(JX)+BXX(M)*RTX(I,N,M)+BXY(M)*RTY(I,N,M)
P(JY)=P(JY)+BYX(M)*RTX(I,N,M)+BYY(M)*RTY(I,N,M)
P(JZ)=P(JZ)+BZX(M)*RTX(I,N,M)+BZY(M)*RTY(I,N,M)
P(JX)=P(JX)-AXX(M)*RUX(I,S,M)-AXY(M)*RUY(I,S,M)-AXZ(M)*RUZ(I,S,M)
P(JY)=P(JY)-AYX(M)*RUX(I,S,M)-AYY(M)*RUY(I,S,M)-AYZ(M)*RUZ(I,S,M)
P(JZ)=P(JZ)-AZX(M)*RUX(I,S,M)-AZY(M)*RUY(I,S,M)-AZZ(M)*RUZ(I,S,M)
GOTO 450
425 A(JX,KX)=A(JX,KX)-BXX(M)
A(JY,KX)=A(JY,KX)-BYX(M)
A(JZ,KX)=A(JZ,KX)-BZX(M)
P(JX)=P(JX)+BXY(M)*RTY(I,N,M)+BXZ(M)*RTZ(I,N,M)
P(JY)=P(JY)+BYY(M)*RTY(I,N,M)+BYZ(M)*RTZ(I,N,M)
P(JZ)=P(JZ)+BZY(M)*RTY(I,N,M)+BZZ(M)*RTZ(I,N,M)
P(JX)=P(JX)-AXX(M)*RUX(I,S,M)-AXY(M)*RUY(I,S,M)-AXZ(M)*RUZ(I,S,M)
P(JY)=P(JY)-AYX(M)*RUX(I,S,M)-AYY(M)*RUY(I,S,M)-AYZ(M)*RUZ(I,S,M)
P(JZ)=P(JZ)-AZX(M)*RUX(I,S,M)-AZY(M)*RUY(I,S,M)-AZZ(M)*RUZ(I,S,M)
GOTO 450
430 A(JX,KY)=A(JX,KY)-BXY(M)
A(JY,KY)=A(JY,KY)-BYY(M)
A(JZ,KY)=A(JZ,KY)-BZY(M)
P(JX)=P(JX)+BXX(M)*RTX(I,N,M)+BXZ(M)*RTZ(I,N,M)
P(JY)=P(JY)+BYX(M)*RTX(I,N,M)+BYZ(M)*RTZ(I,N,M)
P(JZ)=P(JZ)+BZX(M)*RTX(I,N,M)+BZZ(M)*RTZ(I,N,M)
P(JX)=P(JX)-AXX(M)*RUX(I,S,M)-AXY(M)*RUY(I,S,M)-AXZ(M)*RUZ(I,S,M)
P(JY)=P(JY)-AYX(M)*RUX(I,S,M)-AYY(M)*RUY(I,S,M)-AYZ(M)*RUZ(I,S,M)
P(JZ)=P(JZ)-AZX(M)*RUX(I,S,M)-AZY(M)*RUY(I,S,M)-AZZ(M)*RUZ(I,S,M)
GOTO 450
435 P(JX)=P(JX)+BXX(M)*RTX(I,N,M)+BXY(M)*RTY(I,N,M)+BXZ(M)*RTZ(I,N,M)
P(JY)=P(JY)+BYX(M)*RTX(I,N,M)+BYY(M)*RTY(I,N,M)+BYZ(M)*RTZ(I,N,M)
P(JZ)=P(JZ)+BZX(M)*RTX(I,N,M)+BZY(M)*RTY(I,N,M)+BZZ(M)*RTZ(I,N,M)
P(JX)=P(JX)-AXX(M)*RUX(I,S,M)-AXY(M)*RUY(I,S,M)-AXZ(M)*RUZ(I,S,M)
P(JY)=P(JY)-AYX(M)*RUX(I,S,M)-AYY(M)*RUY(I,S,M)-AYZ(M)*RUZ(I,S,M)
P(JZ)=P(JZ)-AZX(M)*RUX(I,S,M)-AZY(M)*RUY(I,S,M)-AZZ(M)*RUZ(I,S,M)
GOTO 450
440 KZ=KZ+ND(I)
KY=KZ-1
KX=KZ-2
DO 441 NN=1,MULTI
441 IF(NODE(I,L,NN).EQ.K) NB=NN
LZ=MS(I)+MULTI*(L-1)+NB
A(JX,KX)=A(JX,KX)+AXX(M)*COSBEX(I,L,1)+AXY(M)*COSBEY(I,L,1)
1 +AXZ(M)*COSBEZ(I,L,1)
A(JX,KY)=A(JX,KY)+AXX(M)*COSBEX(I,L,2)+AXY(M)*COSBEY(I,L,2)
1 +AXZ(M)*COSBEZ(I,L,2)
A(JX,KZ)=A(JX,KZ)+AXX(M)*COSBEX(I,L,3)+AXY(M)*COSBEY(I,L,3)
1 +AXZ(M)*COSBEZ(I,L,3)
A(JY,KX)=A(JY,KX)+AYX(M)*COSBEX(I,L,1)+AYY(M)*COSBEY(I,L,1)
1 +AYZ(M)*COSBEZ(I,L,1)
A(JY,KY)=A(JY,KY)+AYX(M)*COSBEX(I,L,2)+AYY(M)*COSBEY(I,L,2)
1 +AYZ(M)*COSBEZ(I,L,2)
A(JY,KZ)=A(JY,KZ)+AYX(M)*COSBEX(I,L,3)+AYY(M)*COSBEY(I,L,3)
1 +AYZ(M)*COSBEZ(I,L,3)
A(JZ,KX)=A(JZ,KX)+AZX(M)*COSBEX(I,L,1)+AZY(M)*COSBEY(I,L,1)
1 +AZZ(M)*COSBEZ(I,L,1)
A(JZ,KY)=A(JZ,KY)+AZX(M)*COSBEX(I,L,2)+AZY(M)*COSBEY(I,L,2)
1 +AZZ(M)*COSBEZ(I,L,2)
A(JZ,KZ)=A(JZ,KZ)+AZX(M)*COSBEX(I,L,3)+AZY(M)*COSBEY(I,L,3)
1 +AZZ(M)*COSBEZ(I,L,3)
A(JX,LZ)=-G(1)*(BXX(M)*COSBEX(I,L,3)+BXY(M)*COSBEY(I,L,3)
1 +BXZ(M)*COSBEZ(I,L,3))/G(I)
A(JY,LZ)=-G(1)*(BYX(M)*COSBEX(I,L,3)+BYY(M)*COSBEY(I,L,3)
1 +BYZ(M)*COSBEZ(I,L,3))/G(I)
A(JZ,LZ)=-G(1)*(BZX(M)*COSBEX(I,L,3)+BZY(M)*COSBEY(I,L,3)
1 +BZZ(M)*COSBEZ(I,L,3))/G(I)
450 SUMX(JX)=SUMX(JX)-AXX(M)
SUMY(JX)=SUMY(JX)-AXY(M)
SUMZ(JX)=SUMZ(JX)-AXZ(M)
SUMX(JY)=SUMX(JY)-AYX(M)
SUMY(JY)=SUMY(JY)-AYY(M)
SUMZ(JY)=SUMZ(JY)-AYZ(M)
SUMX(JZ)=SUMX(JZ)-AZX(M)
SUMY(JZ)=SUMY(JZ)-AZY(M)
460 SUMZ(JZ)=SUMZ(JZ)-AZZ(M)
DO 470 K=1,NUMBS(I)
DO 470 M=1,8
IF(NORD(I,K,M).NE.J) GOTO 470
L=K
N=M
IF(N.EQ.M) GOTO 475
470 CONTINUE
475 S=NUTYP(I,L)
L=MORD(I,J)
KZ=3*L
KY=KZ-1
KX=KZ-2
SUMX(JX)=SUMX(JX)+REAL(MTYP(I))
SUMY(JY)=SUMY(JY)+REAL(MTYP(I))
SUMZ(JZ)=SUMZ(JZ)+REAL(MTYP(I))
GOTO (480,485,490,495,500,505,510,515,520,520),KOD(I,J)
480 A(JX,KX)=SUMX(JX)
A(JX,KY)=SUMY(JX)
A(JX,KZ)=SUMZ(JX)
A(JY,KX)=SUMX(JY)
A(JY,KY)=SUMY(JY)
A(JY,KZ)=SUMZ(JY)
A(JZ,KX)=SUMX(JZ)
A(JZ,KY)=SUMY(JZ)
A(JZ,KZ)=SUMZ(JZ)
GOTO 600
485 P(JX)=P(JX)-SUMX(JX)*RUX(I,S,N)-SUMY(JX)*RUY(I,S,N)-SUMZ(JX)
1 *RUZ(I,S,N)
P(JY)=P(JY)-SUMX(JY)*RUX(I,S,N)-SUMY(JY)*RUY(I,S,N)-SUMZ(JY)
1 *RUZ(I,S,N)
P(JZ)=P(JZ)-SUMX(JZ)*RUX(I,S,N)-SUMY(JZ)*RUY(I,S,N)-SUMZ(JZ)
1 *RUZ(I,S,N)
GOTO 600
490 A(JX,KX)=SUMX(JX)
A(JX,KY)=SUMY(JX)
A(JY,KX)=SUMX(JY)
A(JY,KY)=SUMY(JY)
A(JZ,KX)=SUMX(JZ)
A(JZ,KY)=SUMY(JZ)
P(JX)=P(JX)-SUMZ(JX)*RUZ(I,S,N)
P(JY)=P(JY)-SUMZ(JY)*RUZ(I,S,N)
P(JZ)=P(JZ)-SUMZ(JZ)*RUZ(I,S,N)
GOTO 600
495 A(JX,KZ)=SUMZ(JX)
A(JY,KZ)=SUMZ(JY)
A(JZ,KZ)=SUMZ(JZ)
P(JX)=P(JX)-SUMX(JX)*RUX(I,S,N)-SUMY(JX)*RUY(I,S,N)
P(JY)=P(JY)-SUMX(JY)*RUX(I,S,N)-SUMY(JY)*RUY(I,S,N)
P(JZ)=P(JZ)-SUMX(JZ)*RUX(I,S,N)-SUMY(JZ)*RUY(I,S,N)
GOTO 600
500 A(JX,KY)=SUMY(JX)
A(JX,KZ)=SUMZ(JX)
A(JY,KY)=SUMY(JY)
A(JY,KZ)=SUMZ(JY)
A(JZ,KY)=SUMY(JZ)
A(JZ,KZ)=SUMZ(JZ)
P(JX)=P(JX)-SUMX(JX)*RUX(I,S,N)
P(JY)=P(JY)-SUMX(JY)*RUX(I,S,N)
P(JZ)=P(JZ)-SUMX(JZ)*RUX(I,S,N)
GOTO 600
505 A(JX,KX)=SUMX(JX)
A(JY,KX)=SUMX(JY)
A(JZ,KX)=SUMX(JZ)
P(JX)=P(JX)-SUMY(JX)*RUY(I,S,N)-SUMZ(JX)*RUZ(I,S,N)
P(JY)=P(JY)-SUMY(JY)*RUY(I,S,N)-SUMZ(JY)*RUZ(I,S,N)
P(JZ)=P(JZ)-SUMY(JZ)*RUY(I,S,N)-SUMZ(JZ)*RUZ(I,S,N)
GOTO 600
510 A(JX,KX)=SUMX(JX)
A(JX,KZ)=SUMZ(JX)
A(JY,KX)=SUMX(JY)
A(JY,KZ)=SUMZ(JY)
A(JZ,KX)=SUMX(JZ)
A(JZ,KZ)=SUMZ(JZ)
P(JX)=P(JX)-SUMY(JX)*RUY(I,S,N)
P(JY)=P(JY)-SUMY(JY)*RUY(I,S,N)
P(JZ)=P(JZ)-SUMY(JZ)*RUY(I,S,N)
GOTO 600
515 A(JX,KY)=SUMY(JX)
A(JY,KY)=SUMY(JY)
A(JZ,KY)=SUMY(JZ)
P(JX)=P(JX)-SUMX(JX)*RUX(I,S,N)-SUMZ(JX)*RUZ(I,S,N)
P(JY)=P(JY)-SUMX(JY)*RUX(I,S,N)-SUMZ(JY)*RUZ(I,S,N)
P(JZ)=P(JZ)-SUMX(JZ)*RUX(I,S,N)-SUMZ(JZ)*RUZ(I,S,N)
GOTO 600
520 KZ=KZ+ND(I)
KY=KZ-1
KX=KZ-2
A(JX,KX)=SUMX(JX)*COSBEX(I,L,1)+SUMY(JX)*COSBEY(I,L,1)
1 +SUMZ(JX)*COSBEZ(I,L,1)
A(JX,KY)=SUMX(JX)*COSBEX(I,L,2)+SUMY(JX)*COSBEY(I,L,2)
1 +SUMZ(JX)*COSBEZ(I,L,2)
A(JX,KZ)=SUMX(JX)*COSBEX(I,L,3)+SUMY(JX)*COSBEY(I,L,3)
1 +SUMZ(JX)*COSBEZ(I,L,3)
A(JY,KX)=SUMX(JY)*COSBEX(I,L,1)+SUMY(JY)*COSBEY(I,L,1)
1 +SUMZ(JY)*COSBEZ(I,L,1)
A(JY,KY)=SUMX(JY)*COSBEX(I,L,2)+SUMY(JY)*COSBEY(I,L,2)
1 +SUMZ(JY)*COSBEZ(I,L,2)
A(JY,KZ)=SUMX(JY)*COSBEX(I,L,3)+SUMY(JY)*COSBEY(I,L,3)
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -