Ignore:
Timestamp:
Dec 21, 2025, 1:35:46 PM (2 days ago)
Author:
Peter A. Buhr <pabuhr@…>
Branches:
master
Children:
4904b05
Parents:
eb0d9b7
Message:

last proofread array chapter

Location:
doc/theses/mike_brooks_MMath/programs/ar-bchk
Files:
6 edited

Legend:

Unmodified
Added
Removed
  • doc/theses/mike_brooks_MMath/programs/ar-bchk/control.c

    reb0d9b7 r80e83b6c  
    2222        size_t n, float x[] ) {
    2323        double sum = 0;
    24         for( size_t i = 0;
     24        for ( size_t i = 0;
    2525            i < BND( n ); i++ )
    26         sum += x[i];
     26        sum += @x[i]@;
    2727        return sum;
    2828}
     
    4343int main() {
    4444        float x[EXPSZ];
    45         for( size_t i = 0; i < EXPSZ; i++ ) x[i] = 0.1 * (i + 1);
    46         for( size_t i = 0; i < EXPSZ; i++ ) printf("elm %zd %g\n", i, x[i]);
     45        for ( size_t i = 0; i < EXPSZ; i++ ) x[i] = 0.1 * (i + 1);
     46        for ( size_t i = 0; i < EXPSZ; i++ ) printf("elm %zd %g\n", i, x[i]);
    4747        double sum_ret = sum( EXPSZ, x );
    4848        printf( "sum   %2g\n", sum_ret );
  • doc/theses/mike_brooks_MMath/programs/ar-bchk/control.cc

    reb0d9b7 r80e83b6c  
    2222        vector<float> & x ) {
    2323        double sum = 0;
    24         for( size_t i = 0;
     24        for ( size_t i = 0;
    2525            i < BND( x.size() ); i++ )
    26         sum += x.at(i);
     26        sum += @x.at(i)@;
    2727        return sum;
    2828}
     
    4343int main() {
    4444        vector<float> x( EXPSZ );
    45         for( size_t i = 0; i < EXPSZ; i++ ) x.at(i) = 0.1 * (i + 1);
    46         for( size_t i = 0; i < EXPSZ; i++ ) cout << "elm " << i << " " <<  x[i] << endl;
     45        for ( size_t i = 0; i < EXPSZ; i++ ) x.at(i) = 0.1 * (i + 1);
     46        for ( size_t i = 0; i < EXPSZ; i++ ) cout << "elm " << i << " " <<  x[i] << endl;
    4747        double sum_ret = sum( x );
    4848        cout << "sum   " << sum_ret << endl;
  • doc/theses/mike_brooks_MMath/programs/ar-bchk/control.cfa

    reb0d9b7 r80e83b6c  
    2222                array(float, N) & x ) {
    2323        double sum = 0;
    24         for( i;
     24        for ( i;
    2525                        BND( N ) )
    26                 sum += x[i];
     26                sum += @x[i]@;
    2727        return sum;
    2828}
  • doc/theses/mike_brooks_MMath/programs/ar-bchk/treatment.c

    reb0d9b7 r80e83b6c  
    2424        float lhs[m][p],
    2525        float rhs[p][n] ) {
    26     for( size_t i = 0;
     26    for ( size_t i = 0;
    2727            i < m; i++ )
    28         for( size_t j = 0;
     28        for ( size_t j = 0;
    2929                j < n; j++ ) {
    3030            res[i][j] = 0.0;
    31             for( size_t k = 0;
     31            for ( size_t k = 0;
    3232                    k < p; k++ )
    33                 res[i][j] +=
    34                     lhs[i][k] *
    35                     rhs[k][j];
     33                @res[i][j] +=@
     34                    @lhs[i][k] *@
     35                    @rhs[k][j];@
    3636        }
    3737}
     
    4949
    5050static void zero( size_t r, size_t c, float mat[r][c] ) {
    51     for( size_t i = 0; i < r; i++ )
    52         for( size_t j = 0; j < c; j++ )
     51    for ( size_t i = 0; i < r; i++ )
     52        for ( size_t j = 0; j < c; j++ )
    5353            mat[i][j] = 0.0;
    5454}
     
    5656
    5757static void fill( size_t r, size_t c, float mat[r][c] ) {
    58     for( size_t i = 0; i < r; i++ )
    59         for( size_t j = 0; j < c; j++ )
     58    for ( size_t i = 0; i < r; i++ )
     59        for ( size_t j = 0; j < c; j++ )
    6060            mat[i][j] = 1.0 * (i + 1) + 0.1 * (j+1);
    6161}
     
    6363
    6464static void print( size_t r, size_t c, float mat[r][c] ) {
    65     for( size_t i = 0; i < r; i++ ) {
    66         for( size_t j = 0; j < c; j++ )
     65    for ( size_t i = 0; i < r; i++ ) {
     66        for ( size_t j = 0; j < c; j++ )
    6767            printf("%2g ", mat[i][j]);
    6868        printf("\n");
  • doc/theses/mike_brooks_MMath/programs/ar-bchk/treatment.cc

    reb0d9b7 r80e83b6c  
    2424        vector<vector<float>> & lhs,
    2525        vector<vector<float>> & rhs ) {
    26     for( size_t i = 0;
     26    for ( size_t i = 0;
    2727            i < res.size(); i++ )
    28         for( size_t j = 0;
     28        for ( size_t j = 0;
    2929                j < res.at(i).size(); j++ ) {
    3030            res.at(i).at(j) = 0.0;
    31             for( size_t k = 0;
     31            for ( size_t k = 0;
    3232                    k < rhs.size(); k++ )
    33                 res.at(i).at(j) +=
    34                     lhs.at(i).at(k) *
    35                     rhs.at(k).at(j);
     33                @res.at(i).at(j) +=@
     34                    @lhs.at(i).at(k) *@
     35                    @rhs.at(k).at(j);@
    3636        }
    3737}
     
    4949
    5050static void zero( vector<vector<float>> & mat ) {
    51     for( size_t i = 0; i < mat.size(); i++ )
    52         for( size_t j = 0; j < mat.at(i).size(); j++ )
     51    for ( size_t i = 0; i < mat.size(); i++ )
     52        for ( size_t j = 0; j < mat.at(i).size(); j++ )
    5353            mat.at(i).at(j) = 0.0;
    5454}
     
    5656
    5757static void fill( vector<vector<float>> & mat ) {
    58     for( size_t i = 0; i < mat.size(); i++ )
    59         for( size_t j = 0; j < mat.at(i).size(); j++ )
     58    for ( size_t i = 0; i < mat.size(); i++ )
     59        for ( size_t j = 0; j < mat.at(i).size(); j++ )
    6060            mat.at(i).at(j) = 1.0 * (i + 1) + 0.1 * (j+1);
    6161}
     
    6363
    6464static void print( vector<vector<float>> & mat ) {
    65     for( size_t i = 0; i < mat.size(); i++ ) {
    66         for( size_t j = 0; j < mat.at(i).size(); j++ )
     65    for ( size_t i = 0; i < mat.size(); i++ ) {
     66        for ( size_t j = 0; j < mat.at(i).size(); j++ )
    6767            cout << mat.at(i).at(j) << " ";
    6868        cout << endl;
  • doc/theses/mike_brooks_MMath/programs/ar-bchk/treatment.cfa

    reb0d9b7 r80e83b6c  
    2424        array(float, M, P) & lhs,
    2525        array(float, P, N) & rhs ) {
    26     for( i; M )
     26    for ( i; M )
    2727
    28         for( j; N ) {
     28        for ( j; N ) {
    2929
    3030            res[i, j] = 0.0;
    31             for( k; P )
     31            for ( k; P )
    3232
    33                 res[i, j] +=
    34                     lhs[i, k] *
    35                     rhs[k, j];
     33                @res[i, j] +=@
     34                    @lhs[i, k] *@
     35                    @rhs[k, j];@
    3636        }
    3737}
Note: See TracChangeset for help on using the changeset viewer.