Revision 1107 trunk/src/df1b2separable/f1b2vc6.cpp
f1b2vc6.cpp (revision 1107)  

2  2 
* $Id$ 
3  3 
* 
4  4 
* Author: David Fournier 
5 
* Copyright (c) 20082012 Regents of the University of California


5 
* Copyright (c) 20082012 Regents of the University of California 

6  6 
*/ 
7  7 
/** 
8  8 
* \file 
...  ...  
20  20 
int rmin=M.indexmin(); 
21  21 
int rmax=M.indexmax(); 
22  22 
df1b2matrix tmp(rmin,rmax); 
23 
for (int i=rmin;i<=rmax;i++)


23 
for (int i=rmin;i<=rmax;i++) 

24  24 
{ 
25  25 
int cmin=M(i).indexmin(); 
26  26 
int cmax=M(i).indexmax(); 
27  27 
tmp(i).noallocate(cmin,cmax); 
28 
for (int j=cmin;j<=cmax;j++)


28 
for (int j=cmin;j<=cmax;j++) 

29  29 
tmp(i,j)=exp(M(i,j)); 
30  30 
} 
31  31 
return tmp; 
...  ...  
44  44 
int cmax=M(rmin).indexmax(); 
45  45 
df1b2vector tmp(cmin,cmax); 
46  46 
tmp.initialize(); 
47 
for (int i=rmin;i<=rmax;i++)


47 
for (int i=rmin;i<=rmax;i++) 

48  48 
{ 
49 
for (int j=cmin;j<=cmax;j++)


49 
for (int j=cmin;j<=cmax;j++) 

50  50 
tmp(j)+=M(i,j); 
51  51 
} 
52  52 
return tmp; 
...  ...  
63  63 
int rmax=M.indexmax(); 
64  64 
df1b2vector tmp(rmin,rmax); 
65  65 
tmp.initialize(); 
66 
for (int i=rmin;i<=rmax;i++)


66 
for (int i=rmin;i<=rmax;i++) 

67  67 
{ 
68  68 
int cmin=M(rmin).indexmin(); 
69  69 
int cmax=M(rmin).indexmax(); 
70 
for (int j=cmin;j<=cmax;j++)


70 
for (int j=cmin;j<=cmax;j++) 

71  71 
tmp(i)+=M(i,j); 
72  72 
} 
73  73 
return tmp; 
Also available in: Unified diff