
Listings 425
53 return 1+k;
54 }
55 }
56
57 for (Size k=0; k < n ; k++) {
58 iL(k,k) = 1 / L(k,k);
59
60 for (Size j=k+1; j<n; ++j){
61 iL(j , k) = − inner prod (
62 project(row(L, j ) , range(k, j )) ,
63 project(column(iL , k), range(k, j )) )
64 /L(j,j);
65 }
66 }
67
68 return 0;
69 }
70
71 }
B.8.10 Null Types
Listing B.105: The header file of null types.
1 // cl / u t i l i t i e s / n u l l . hpp
2 #ifndef CLUSLIB NULL HPP
3 #d e f i n e CLUSLIB NULL HPP
4
5 #include<c l / types . hpp>
6
7 namespace ClusLib {
8
9 template <class Type>
10 class Null ;
11
12 templa te <>
13 class Null<Integer> {
14 public :
15 Null () {}
16 operator Integer () con st {
17 return I n t e g e