| 800 | } |
| 801 | |
| 802 | void |
| 803 | PFEMElement2DBubble::getL(Matrix& l) const |
| 804 | { |
| 805 | Matrix Gbub(2,3); |
| 806 | getGbub(Gbub); |
| 807 | double Mbub = getMbub(); |
| 808 | |
| 809 | Matrix Kbub(2,2); |
| 810 | getKbub(Kbub); |
| 811 | |
| 812 | if (ops_Dt > 0) { |
| 813 | Kbub(0,0) += Mbub/ops_Dt; |
| 814 | Kbub(1,1) += Mbub/ops_Dt; |
| 815 | } |
| 816 | |
| 817 | Matrix invKbub(2,2); |
| 818 | Kbub.Invert(invKbub); |
| 819 | |
| 820 | l.resize(3,3); |
| 821 | l.addMatrixTripleProduct(0.0, Gbub, invKbub, 1.0); |
| 822 | } |
| 823 | |
| 824 | void |
| 825 | PFEMElement2DBubble::getFp(Vector& fp) const |
nothing calls this directly
no test coverage detected