edit makefile

doc/README.md

+## Praktikum 02
+# OpenMPI - Radix Sort
+
+# How To
+Run `make` then open `sort` using `mpirun`
+
+# Pembagian Tugas
+13516100 - Putu Gery Wahyu Nugraha
+* Counting Sort
+* Radix Sort
+* All realted with MPIs
+* Doc and README
+
+13516091 - Yasya Rusyda Aslina
+* Desain solusi paralel
+* Testing
+
+# Deskripsi Solusi Paralel
+OpenMPI merupakan sebuah message passing interface, oleh karena itu kemampuan terbesarnya adalah mengirimkan pesan antar proses.
+
+Perlu dibuat sebuah solusi radix sort yang menggunakan sistem pengiriman pesan antar proses. Salah satu cara termudahnya adalah
+
+1. Bagi data yang akan di sort kepada setiap proses secara adil (data lokal)
+2. Hitung nilai count dari setiap data yang akan di sort
+3. Sort data lokal pada tiap-tiap proses
+4. Gabungkan semua data lokal yang ada dari proses
+5. Gabungkan nilai count dari tiap proses
+6. Lakukan sort secara global
+7. Ulangi langkah 1 
+
+Dari langkah-langkah tersebut
+
+* Langkah 1 bisa kita selesaikan dengan menggunakan `MPI_Scatter`
+* Langkah 4 bisa kita selesaikan dengan menggunakan `MPI_Gather`
+* Langkah 5 bisa kita selesaikan dengan menggunakan `MPI_Reduce` dan `MPI_SUM`
+
+# Analisis Solusi
+Setelah penulis menulis kode diatas, penulis menemukan banyak hal yang dapat diefektifkan. Salah satunya adalah bagaimana program ini melakukan sorting. Jika kita melihat langkah yang telah dideskripsikan diatas, maka kita menemukan bahwa: 
+* Sort dilakukan dua kali, dalam scope lokal dan scope global. Tidak efektif, kemungkinan besar bisa dilakukan sekali saja
+
+# Jumlah Thread yang Digunakan
+## 2
+Ketika melakukan pengerjaan melalui komputer pribadi, kami menemukan bahwa jika nilai thread lebih dari 2, ada kemungkinan program MPI tersebut stuck. Walaupun hal ini tidak pernah kami temukan di komputer remote. Maka dari itu untuk mendapatkan data yang konsisten kami memutuskan untuk menggunakan 2 thread saja
+
+# Pengukuran Kinerja
+| N | Size | 1     | 2 | 3 | avg |
+|---|------|-------|---|---|-----|
+| 1 | 5000 | 15145 | 15131 | 15273 | 15183 | 
+| 2 | 5000 | 40787 | 11316 | 5893  | 19332 |
+| 1 | 50000| 187103| 148275| 105336 | 146904 |
+| 2 | 50000| 89551 | 130345| 119153 | 113016 |
+| 1 |100000| 224312|154683 |154609 |177868 | 
+| 2 |100000| 217054|197716|150975|188581 |
+| 1 |200000|621111|629161|706030|652100 |
+| 2 |200000|304229|358362|302638|321743 |
+| 1 |400000|1142259|728879|1230207|1033781 |
+| 2 |400000|1339009|1217159|1199709|1251959 |
+
+# Analisis Pengukuran Kinerja
+Dari data kinerja diatas, sebenarnya kami lebih cenderung merasa bahwa faktor terbesar berkurangnya waktu bukanlah jumlah thread yang digunakan, tapi bagaimana kondisi cache dari prosesor. Walaupun dari data diatas dapat kita temukan perbedaan, tapi perbedaan tersebut sifatnya random dan tidak memiliki arah yang jelas. Oleh karena itu setiap pengurangan waktu yang terjadi ketika jumlah thread ditambahkan (atau penambahan waktu) lebih cenderung disebabkan oleh kondisi cache pada saat itu.
+
+Tapi karena radix sort parallel jauh lebih unggul ketika jumlah data banyak maka mungkin saja penurunan drastis yang terjadi ketika menggunakan dua thread ketika size = 200000 disebabkan oleh penggunaan parallel programming dan bukan kondisi cache saja.

src/Makefile

-build:
-	gcc-8 -g -Wall -o radix_sort radix_sort.c -fopenmp -lm 
+all:
+	nvcc radix_sort_cuda.cu -o radix_sort_cuda
+
+clean:
+	rm -rf radix_sort_cuda
+	rm -rf radix_sort_cuda.exe
\ No newline at end of file

src/radix_par2.cpp

-// C++ implementation of Radix Sort 
-#include <iostream> 
-#include <cstdlib>
-#include <time.h>
-using namespace std; 
-  
-// A utility function to get maximum value in arr[] 
-int getMax(int arr[], int n) 
-{ 
-    int mx = arr[0]; 
-    for (int i = 1; i < n; i++) 
-        if (arr[i] > mx) 
-            mx = arr[i]; 
-    return mx; 
-} 
-  
-// A function to do counting sort of arr[] according to 
-// the digit represented by exp. 
-__global__ void countSort(int arr[], int n, int exp) 
-{ 
-    int output[n]; // output array 
-    int i, count[10] = {0}; 
-  
-    // Store count of occurrences in count[] 
-    for (i = 0; i < n; i++) 
-        count[ (arr[i]/exp)%10 ]++; 
-  
-    // Change count[i] so that count[i] now contains actual 
-    //  position of this digit in output[] 
-    for (i = 1; i < 10; i++) 
-        count[i] += count[i - 1]; 
-  
-    // Build the output array 
-    for (i = n - 1; i >= 0; i--) 
-    { 
-        output[count[ (arr[i]/exp)%10 ] - 1] = arr[i]; 
-        count[ (arr[i]/exp)%10 ]--; 
-    } 
-  
-    // Copy the output array to arr[], so that arr[] now 
-    // contains sorted numbers according to current digit 
-    for (i = 0; i < n; i++) 
-        arr[i] = output[i]; 
-} 
-  
-// The main function to that sorts arr[] of size n using  
-// Radix Sort 
-void radixsort(int arr[], int n) 
-{ 
-    int *d_arr;
-    // Find the maximum number to know number of digits 
-    int m = getMax(arr, n); 
-  
-    // Do counting sort for every digit. Note that instead 
-    // of passing digit number, exp is passed. exp is 10^i 
-    // where i is current digit number 
-
-    cudaMalloc((void**)&d_arr,sizeof(int)*n);
-    cudaMemcpy(d_arr, arr, sizeof(int)*n,cudaMemcpyHostToDevice);
- 
-    for (int exp = 1; m/exp > 0; exp *= 10) {
-        countSort<<<1,32>>>(arr, n, exp);
-        cudaMemcpy(arr, d_arr, sizeof(int)*n, cudaMemcpyDeviceToHost);
-    }
-
-    cudaFree(d_arr); 
-} 
-  
-// A utility function to print an array 
-void print(int arr[], int n) 
-{ 
-    for (int i = 0; i < n; i++) 
-        cout << arr[i] << "\n"; 
-} 
-  
-void rng(int* arr,int n){
-    int seed = 13516067;
-    srand(seed);
-    for (long i=0; i<n; i++){
-        arr[i] = (int) rand();
-    }
-}
-
- 
-timespec diff(timespec start, timespec end)
-{
-    timespec temp;
-    if ((end.tv_nsec - start.tv_nsec) < 0) {
-        temp.tv_sec = end.tv_sec - start.tv_sec - 1;
-        temp.tv_nsec = 1000000000 + end.tv_nsec - start.tv_nsec;
-    } else {
-        temp.tv_sec = end.tv_sec - start.tv_sec;
-        temp.tv_nsec = end.tv_nsec - start.tv_nsec;
-    }
-    return temp;
-}
-
-// Driver program to test above functions 
-int main() 
-{ 
-    timespec start, stop;
-    
-    int n;
-    cout<<"Masukkan nilai N\n";
-    cin>>n;
-    int arr[n];
-    rng(arr,n);
-    clock_gettime(CLOCK_REALTIME, &start);
-    radixsort(arr,n);
-    clock_gettime(CLOCK_REALTIME, &stop);
-    print(arr,n);
-    
-    timespec duration = diff(start, stop);
-    long time = duration.tv_sec * 1000000 + duration.tv_nsec/1000;
-    printf("\n%d.%09d s\n", duration.tv_sec, duration.tv_nsec);
-    
-    return 0; 
-} 
\ No newline at end of file