version 1.13 | version 1.14 |
---|
| |
////////////////////////////////////////////////// | ////////////////////////////////////////////////// |
// Force dAi_drk | // Force dAi_drk |
////////////////////////////////////////////////// | ////////////////////////////////////////////////// |
BigReal dAjkdrjk = PI*rj*rjk_1*(rjk_1*rjk_1*(rj*rj-ri*ri) - 1.f);//ef' | BigReal dAjkdrjk = PI*rj*rjk_1*(rjk_1*rjk_1*(rj*rj-rk*rk) - 1.f);//ef' |
BigReal dAjkdrjkdxj = -dAjkdrjk*dxjk; // e f h' | BigReal dAjkdrjkdxj = -dAjkdrjk*dxjk; // e f h' |
BigReal dAjkdrjkdyj = -dAjkdrjk*dyjk; | BigReal dAjkdrjkdyj = -dAjkdrjk*dyjk; |
BigReal dAjkdrjkdzj = -dAjkdrjk*dzjk; | BigReal dAjkdrjkdzj = -dAjkdrjk*dzjk; |