cos (Core Functions)


Import

#include <archr/core/cos.h>
#include <archr/core/cos.hpp>
import archr.core
atomsLoads('archr')

Synopsis

 float archr_core_cos_1f32(const float i1);
(1)
 double archr_core_cos_1f64(const double i1);
(2)
 void archr_core_cos_f32(float *o0, const float *i1, size_t sz);
(3)
 void archr_core_cos_f64(double *o0, const double *i1, size_t sz);
(4)
 float cos(const float i1);
(1)
 double cos(const double i1);
(2)
 void cos(const float *i1, size_t sz, float *o0);
(3)
 void cos(const double *i1, size_t sz, double *o0);
(4)
template <typename Range> void cos(const Range& i1, Range& o0);
(5)
subroutine archr_core_cos_1f32(real(4) :: r, real(4), parameter :: i1)
(1)
subroutine archr_core_cos_1f64(real(8) :: r, real(8), parameter :: i1)
(2)
subroutine archr_core_cos_f32(real(4), dimension(*), parameter :: i1, integer(4) :: sz, real(4), dimension(*) :: o0)
(3)
subroutine archr_core_cos_f64(real(8), dimension(*), parameter :: i1, integer(4) :: sz, real(8), dimension(*) :: o0)
(4)
def cos(i1):
    return o0
(1)
function o0 = archr_core_cos(i1)
(1)

Description

Computes the Cosine in Radians:

Parameters

i1_first, i1_last

The range of input elements

o0

The beginning of the destination range, may be equal to i1

i1

The scalar/contiguous data input element

Example


#include <time.h>
#include <stdio.h>
#include <stdlib.h>
#include <stdint.h>
#include <stddef.h>
#include <archr/core/cos.h>

double rand_in(double min, double max) {
  return min + ((double)(max - min) * ((double)rand() / RAND_MAX));
}

int main() {
  time_t t;
  size_t i;
  size_t sz = 10;
  float* i1 = (float*)malloc(sz * sizeof(float));
  float* o1 = (float*)malloc(sz * sizeof(float));

  srand((unsigned int)time(&t));
  for (i = 0; i < sz; ++i) {
    i1[i] = (float)rand_in((float)-1, (float)1);
  }
  archr_core_cos_f32(o1, i1, sz);
  for (i = 0; i < sz; ++i) {
    printf("%4lu: cos(%f) = %f\n", i, (double)i1[i], (double)o1[i]);
  }
}

Possible Output

   0: cos(-0.782005) = 0.709502
   1: cos(-0.654290) = 0.793480
   2: cos(-0.937116) = 0.592115
   3: cos(0.569177) = 0.842345
   4: cos(0.616947) = 0.815649
   5: cos(0.746037) = 0.734385
   6: cos(-0.236809) = 0.972091
   7: cos(-0.672507) = 0.782262
   8: cos(0.890370) = 0.629124
   9: cos(-0.545268) = 0.854988

#include <ctime>
#include <cstdio>
#include <cstdlib>
#include <cstdint>
#include <cstddef>
#include <vector>
#include <algorithm>
#include <iomanip>
#include <iostream>
#include <archr/core/cos.hpp>

double rand_in(double min, double max) {
  return min + (double(max - min) * (double(std::rand()) / RAND_MAX));
}

int main() {
  std::size_t i;
  std::size_t sz = 10;
  std::vector<float> i1(sz);
  std::vector<float> o1(sz);

  std::srand(std::time(0));
  std::generate(i1.begin(), i1.end(), []() { return rand_in(float(-1), float(1)); });
  archr::core::cos(i1.data(), sz, o1.data());
  for (i = 0; i < sz; ++i) {
    std::cout << std::setw(4) << i << ": " << "cos" << "(" << i1[i] << ")" << " = " << o1[i] << std::endl;
  }
}

Possible Output

   0: cos(0.942286) = 0.58794
   1: cos(0.183788) = 0.983158
   2: cos(-0.326225) = 0.947259
   3: cos(-0.782217) = 0.709352
   4: cos(-0.500638) = 0.877276
   5: cos(-0.805106) = 0.693035
   6: cos(0.494578) = 0.880169
   7: cos(0.934714) = 0.594048
   8: cos(0.510651) = 0.872426
   9: cos(0.634508) = 0.805364

program main
  integer(4), parameter  :: sz = 10
  real(4), dimension(sz) :: o0
  real(4), dimension(sz) :: i0
  real(8)                :: r0, r1
  real(8)                :: min0, max0
  real(8)                :: min1, max1
  real(8)                :: min2, max2
  ! Init:
    min0 = -1
    max0 = 1
  do i=1,sz
    r0 = random_in(min0, max0)
    i0(i) = r0
  end do
  ! Example:
  call archr_core_cos_f32(o0, i0, size(i0))
  ! Output:
  do i=1,sz
    print *, i, ": ", o0(i)
  end do
contains
  ! Generate a random number within a range
  function random_in(mn, mx) result(r)
    real(8) :: r
    real(8), intent(in) :: mn, mx
    r = mn + (rand() * (mx - mn))
  end function random_in
end program

Possible Output

           1 :   0.540315151    
           2 :   0.740538716    
           3 :   0.872152925    
           4 :   0.996582329    
           5 :   0.997853398    
           6 :   0.846147418    
           7 :   0.616968930    
           8 :   0.936694443    
           9 :   0.936391711    
          10 :   0.645295918    
atomsLoad('archr');

// Increase the stacksize
v = getversion('scilab');
if v(1) < 6 stacksize('max'); end

// The number of internal iterations
iter = 2000;

//
sizes = [];
times = [];
speed = [];

// Main loop
for n = 10:17
    sz = 2^n;
    I  = rand(1, sz);
    timer();
    for j = 1:iter
        R1 = cos(I);
    end
    T1 = timer();
    
    timer();
    for j = 1:iter
        R2 = archr_core_cos(I);
    end
    T2 = timer();

    sizes = [sizes n];
    times = [times; T1 T2];
    speed = [speed; 1 (T1 / T2)];

    printf("\n");
    printf("============================\n");
    printf("Size is: %d (2^%d)\n", sz, n);
    printf("\n");
    printf("Raw timing:\n");
    printf("  SciLab: %f\n", T1);
    printf("  Arch-R: %f\n", T2);
    printf("----------------------------\n");
    printf("Speedup:\n");
    printf("  x%f\n", T1 / T2);
    printf("============================\n");
end

// Plot timings
scf(0);
bar(sizes, times);
xlabel("Size (power-of-2)");
ylabel("Time (s)");
legend(["SciLab"; "Archr-R"], 'in_upper_left');
title("Raw timings between Archr-R and SciLab");

// Plot speedup
scf(1);
bar(sizes, speed);
xlabel("Size (power-of-2)");
ylabel("Factor of speed (n times faster than)");
legend(["SciLab"; "Archr-R"], 'in_upper_left');
title("Speedup factor betweeb Archr-R and SciLab");

Possible Output

Start Toolbox archr
	Load macros
	Load gateways


============================
Size is: 1024 (2^10)

Raw timing:
  SciLab: 0.048000
  Arch-R: 0.020000
----------------------------
Speedup:
  x2.400000
============================

============================
Size is: 2048 (2^11)

Raw timing:
  SciLab: 0.100000
  Arch-R: 0.024000
----------------------------
Speedup:
  x4.166667
============================

============================
Size is: 4096 (2^12)

Raw timing:
  SciLab: 0.196000
  Arch-R: 0.040000
----------------------------
Speedup:
  x4.900000
============================

============================
Size is: 8192 (2^13)

Raw timing:
  SciLab: 0.404000
  Arch-R: 0.084000
----------------------------
Speedup:
  x4.809524
============================

============================
Size is: 16384 (2^14)

Raw timing:
  SciLab: 0.808000
  Arch-R: 0.164000
----------------------------
Speedup:
  x4.926829
============================

============================
Size is: 32768 (2^15)

Raw timing:
  SciLab: 1.632000
  Arch-R: 0.344000
----------------------------
Speedup:
  x4.744186
============================

============================
Size is: 65536 (2^16)

Raw timing:
  SciLab: 3.260000
  Arch-R: 0.648000
----------------------------
Speedup:
  x5.030864
============================

============================
Size is: 131072 (2^17)

Raw timing:
  SciLab: 6.532000
  Arch-R: 1.320000
----------------------------
Speedup:
  x4.948485
============================