c_cpp/etc/calc/cal/sumtimes.cal

   1 /*
   2  * sumtimes - runtimes evaluating sums & squares of large lists and mats
   3  *
   4  * Copyright (C) 2006  Ernest Bowen
   5  *
   6  * Calc is open software; you can redistribute it and/or modify it under
   7  * the terms of the version 2.1 of the GNU Lesser General Public License
   8  * as published by the Free Software Foundation.
   9  *
  10  * Calc is distributed in the hope that it will be useful, but WITHOUT
  11  * ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
  12  * or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General
  13  * Public License for more details.
  14  *
  15  * A copy of version 2.1 of the GNU Lesser General Public License is
  16  * distributed with calc under the filename COPYING-LGPL.  You should have
  17  * received a copy with calc; if not, write to Free Software Foundation, Inc.
  18  * 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA.
  19  *
  20  * @(#) $Revision: 30.1 $
  21  * @(#) $Id: sumtimes.cal,v 30.1 2007/03/16 11:09:54 chongo Exp $
  22  * @(#) $Source: /usr/local/src/bin/calc/cal/RCS/sumtimes.cal,v $
  23  *
  24  * Under source code control:   2006/06/22 17:29
  25  * File existed as early as:    2006
  26  *
  27  * Share and enjoy!  :-)        http://www.isthe.com/chongo/tech/comp/calc/
  28  */
  29
  30
  31 global sumtimes_t0, sumtimes_t1, sumtimes_t2, sumtimes_t3;
  32 global sumtimes_A, sumtimes_B;
  33 config("tilde", 0),;
  34
  35 define timematsum(N) {
  36     local n, s, p, ptop;
  37
  38     sumtimes_A = mat[N];
  39
  40     for (n = 0; n < N; n++) sumtimes_A[n] = rand(N);
  41
  42     ptop = &sumtimes_A[n-1];
  43     sumtimes_t0 = usertime();
  44     for (s = n = 0; n < N; n++) s += sumtimes_A[n];
  45     sumtimes_t1 = usertime();
  46     for (s = 0, p = &sumtimes_A[0]; p <= ptop; p++) s += *p;
  47     sumtimes_t2 = usertime();
  48     s = matsum(sumtimes_A);
  49     sumtimes_t3 = usertime();
  50
  51     print "Matrix sum runtimes";
  52     printf('\tStandard "for" loop:\t\t%.4f\n', sumtimes_t1 - sumtimes_t0);
  53     printf('\t"For" loop using pointers:\t\t%.4f\n', sumtimes_t2 - sumtimes_t1);
  54     printf('\tUsing builtin "matsum":\t\t%.4f\n', sumtimes_t3 - sumtimes_t2);
  55 }
  56
  57 define timelistsum(N) {
  58     local n, s;
  59
  60     sumtimes_A = makelist(N);
  61     for (n = 0; n < N; n++) sumtimes_A[n] = rand(N);
  62
  63     sumtimes_t0 = usertime();
  64     for (s = n = 0; n < N; n++) s += sumtimes_A[n];
  65     sumtimes_t1 = usertime();
  66     s = sum(sumtimes_A);
  67     sumtimes_t2 = usertime();
  68     print "List sum runtimes";
  69     printf('\tStandard "for" loop:\t\t%.4f\n', sumtimes_t1 - sumtimes_t0);
  70     printf('\tUsing builtin "sum":\t\t%.4f\n', sumtimes_t2 - sumtimes_t1);
  71 }
  72
  73
  74 define timematsort(N) {
  75     local n;
  76
  77     sumtimes_A = mat[N];
  78     for (n = 0; n < N; n++) sumtimes_A[n] = rand(N);
  79     sumtimes_t0 = usertime();
  80     sort(sumtimes_A);
  81     sumtimes_t1 = usertime();
  82     printf('\tMatrix sort runtime:\t\t%.4f\n', sumtimes_t1 - sumtimes_t0);
  83 }
  84
  85
  86 define timelistsort(N) {
  87     local n;
  88
  89     sumtimes_A = makelist(N);
  90     for (n = 0; n < N; n++) sumtimes_A[n] = rand(N);
  91     sumtimes_t0 = usertime();
  92     sort(sumtimes_A);
  93     sumtimes_t1 = usertime();
  94     printf('\tList sort runtime:\t\t%.4f\n', sumtimes_t1 - sumtimes_t0);
  95 }
  96
  97 define timematreverse(N) {
  98     local n;
  99
 100     sumtimes_A = mat[N];
 101     for (n = 0; n < N; n++) sumtimes_A[n] = rand(N);
 102     sumtimes_t0 = usertime();
 103     reverse(sumtimes_A);
 104     sumtimes_t1 = usertime();
 105     printf('\tMatrix reverse runtime        %.4f\n', sumtimes_t1 - sumtimes_t0);
 106 }
 107
 108 define timelistreverse(N) {
 109     local n;
 110
 111     sumtimes_A = makelist(N);
 112     for (n = 0; n < N; n++) sumtimes_A[n] = rand(N);
 113     sumtimes_t0 = usertime();
 114     reverse(sumtimes_A);
 115     sumtimes_t1 = usertime();
 116     printf('\tList reverse runtime:\t\t%.4f\n', sumtimes_t1 - sumtimes_t0);
 117 }
 118
 119 define timematssq(N) {
 120     local n, s, p, ptop;
 121
 122     sumtimes_A = mat[N];
 123
 124     for (n = 0; n < N; n++) sumtimes_A[n] = rand(N);
 125
 126     ptop = &sumtimes_A[n-1];
 127     sumtimes_t0 = usertime();
 128     for (s = n = 0; n < N; n++) s += sumtimes_A[n]^2;
 129     sumtimes_t1 = usertime();
 130     for (s = 0, p = &sumtimes_A[0]; p <= ptop; p++) s += (*p)^2;
 131     sumtimes_t2 = usertime();
 132
 133     print "Matrix sum of squares runtimes";
 134     printf('\tStandard "for" loop:\t\t%.4f\n', sumtimes_t1 - sumtimes_t0);
 135     printf('\t"For" loop using pointers:\t\t%.4f\n', sumtimes_t2 - sumtimes_t1);
 136 }
 137
 138 define timelistssq(N) {
 139     local n, s;
 140
 141     sumtimes_A = makelist(N);
 142     for (n = 0; n < N; n++) sumtimes_A[n] = rand(N);
 143
 144     sumtimes_t0 = usertime();
 145     for (s = n = 0; n < N; n++) s += sumtimes_A[n]^2;
 146     sumtimes_t1 = usertime();
 147     s = ssq(sumtimes_A);
 148     sumtimes_t2 = usertime();
 149     print "List sum of squares runtimes";
 150     printf('\tStandard "for" loop:\t\t%.4f\n', sumtimes_t1 - sumtimes_t0);
 151     printf('\tUsing builtin "ssq":\t\t%.4f\n', sumtimes_t2 - sumtimes_t1);
 152 }
 153
 154 define timehmean(N, M = 10) {
 155     local n, s, v1, v2;
 156
 157     sumtimes_A = makelist(N);
 158     for (n = 0; n < N; n++) sumtimes_A[n] = rand(1, M);
 159
 160     sumtimes_t0 = usertime();
 161     for (s = n = 0; n < N; n++) s += 1/sumtimes_A[n];
 162     v1 = N/s;
 163     sumtimes_t1 = usertime();
 164     v2 = hmean(sumtimes_A);
 165     sumtimes_t2 = usertime();
 166     print v1, v2;
 167     print "List harmonic meanruntimes";
 168     printf('\tStandard "for" loop:\t\t%.4f\n', sumtimes_t1 - sumtimes_t0);
 169     printf('\tUsing builtin "hmean":\t\t%.4f\n', sumtimes_t2 - sumtimes_t1);
 170 }
 171
 172 define doalltimes(N) {
 173     timematsum(N);
 174     print;
 175     timelistsum(N);
 176     print;
 177     timematssq(N);
 178     print;
 179     timelistssq(N);
 180     print;
 181     timematsort(N);
 182     timelistsort(N);
 183     timematreverse(N);
 184     timelistreverse(N);
 185     print;
 186 }