18 template <class RT>
19 int rank_11(const RT& a0)  in rank_11()
24 template <class RT>
25 int rank_21(const RT& a0, const RT& a1)  in rank_21()
30 template <class RT>
31 int rank_12(const RT& a0, const RT& a1)  in rank_12()
36 template <class RT>
37 int rank_31(const RT& a0, const RT& a1, const RT& a2)  in rank_31()
42 template <class RT>
43 int rank_32(const RT& a0, const RT& b0,  in rank_32()
44             const RT& a1, const RT& b1,  in rank_32()
45             const RT& a2, const RT& b2)  in rank_32()
53         return rank_31<RT>(b0,b1,b2);  in rank_32()
57         return 1 + rank_21<RT>(b0, b1);  in rank_32()
62       return 1 + rank_21<RT>(b0, a1*b2-a2*b1);  in rank_32()
67     return 1 + rank_21<RT>(a0*b1-a1*b0, a0*b2-a2*b0);  in rank_32()
71 template <class RT>
72 int rank_22(const RT& a0, const RT& b0,  in rank_22()
73             const RT& a1, const RT& b1)  in rank_22()
79       return rank_21<RT>(b0,b1);  in rank_22()
82       return 1 + rank_11<RT>(b0);  in rank_22()
84   return 1 + rank_11<RT>(a0*b1-a1*b0);  in rank_22()
87 template <class RT>
88 int rank_33(const RT& a0, const RT& b0, const RT& c0,  in rank_33()
89             const RT& a1, const RT& b1, const RT& c1,  in rank_33()
90             const RT& a2, const RT& b2, const RT& c2)  in rank_33()
98         return rank_32<RT>(b0, c0, b1, c1, b2, c2);  in rank_33()
102         return 1 + rank_22<RT>(b0, c0, b1, c1);  in rank_33()
106       return 1 + rank_22<RT>(b0, c0, a1*b2-a2*b1, a1*c2-a2*c1);  in rank_33()
110     return 1 + rank_22<RT>(a0*b1-a1*b0, a0*c1-a1*c0, a0*b2-a2*b0, a0*c2-a2*c0);  in rank_33()
115 template <class RT>
116 int rank_23(const RT& a0, const RT& b0, const RT& c0,  in rank_23()
117             const RT& a1, const RT& b1, const RT& c1)  in rank_23()
123       return rank_22<RT>(b0, c0, b1, c1);  in rank_23()
126       return 1 + rank_12<RT>(b0,c0);  in rank_23()
130     return 1 + rank_12<RT>(a0*b1-a1*b0,a0*c1-a1*c0);  in rank_23()
135 template <class RT>
136 int rank_34(const RT& a0, const RT& b0, const RT& c0, const RT& d0,  in rank_34()
137             const RT& a1, const RT& b1, const RT& c1, const RT& d1,  in rank_34()
138             const RT& a2, const RT& b2, const RT& c2, const RT& d2)  in rank_34()
146         return rank_33<RT>(b0, c0, d0, b1, c1, d1, b2, c2, d2);  in rank_34()
150         return 1 + rank_23<RT>(b0, c0, d0, b1, c1, d1);  in rank_34()
154       return 1 + rank_23<RT>(b0, c0, d0,a1*b2-a2*b1, a1*c2-a2*c1, a1*d2 - a2*d1);  in rank_34()
158 …return 1 + rank_23<RT>(a0*b1-a1*b0, a0*c1-a1*c0, a0*d1-a1*d0, a0*b2-a2*b0, a0*c2-a2*c0, a0*d2-a2*d…  in rank_34()