Changeset 789 in Sophya for trunk/SophyaProg/Tests/lpk.cc
- Timestamp:
- Mar 20, 2000, 7:24:31 PM (26 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SophyaProg/Tests/lpk.cc
r786 r789 29 29 for(j=0; j<n; j++) a(i,j,0) = GauRnd(0., 1.); 30 30 31 TArray<r_4> x( 1,n), b(n,n);31 TArray<r_4> x(n,1), b(n,1); 32 32 r_4 sum ; 33 for(i=0; i<n; i++) x( 0,i,0) = GauRnd(2., 1.5);33 for(i=0; i<n; i++) x(i,0,0) = GauRnd(2., 1.5); 34 34 for(i=0; i<n; i++) { 35 35 sum = 0.; 36 for(j=0; j<n; j++) sum += a( j,i,0)*x(0,j,0);37 b( 0,i,0) = sum;36 for(j=0; j<n; j++) sum += a(i,j,0)*x(j,0,0); 37 b(i,0,0) = sum; 38 38 } 39 39
Note:
See TracChangeset
for help on using the changeset viewer.