[Παράλληλα & Διανεμημένα] 3η εργασία

[Thmmy.gr portal]

  Thmmy.gr portal
   Forum
   Downloads
   Ενεργ. Λογαριασμού
   Επικοινωνία
  
  Χρήσιμα links
   Σελίδα τμήματος
   Βιβλιοθήκη Τμήματος
   Elearning
   Φοιτητικά fora
   Πρόγραμμα Λέσχης
   Πρακτική Άσκηση
   Ηλεκτρονική Εξυπηρέτηση Φοιτητών
   Διανομή Συγγραμμάτων
   Ψηφιακό Καταθετήριο Διπλωματικών
   Πληροφορίες Καθηγητών
   Instagram @thmmy.gr
   mTHMMY
  
  Φοιτητικές Ομάδες
   ACM
   Aristurtle
   ART
   ASAT
   BEAM
   BEST Thessaloniki
   EESTEC LC Thessaloniki
   EΜΒ Auth
   IAESTE Thessaloniki
   IEEE φοιτητικό παράρτημα ΑΠΘ
   SpaceDot
   VROOM
   Panther
  

THMMY.gr Web

Helios_MultiScribbles2SMF Default ThemeSMFone_Blue

[Ανεβάζετε τα θέματα των εξετάσεων στον τομέα Downloads]

Είσαι πρωτοετής;... Καλώς ήρθες! Μπορείς να βρεις πληροφορίες εδώ. Βοήθεια για τους καινούργιους μέσω χάρτη.

Κατεβάστε εδώ το Android Application για εύκολη πρόσβαση στο forum.

Ανεβάζετε τα θέματα των εξετάσεων στον τομέα Downloads με προσοχή στα ονόματα των αρχείων!

[12:00-13:30]

Η γραμματεία είναι ανοιχτή καθημερινά 12:00-13:30

[DoomGuard]

Παιδια τι παιζει με τα ορια της shared? Πως ξερουμε ποσα bytes χωραει?
 Π.χ. ολον τον Τ τον χωραει ή πολλα ζηταω...

ο διάδης έχει 16384 bytes  shared memory ανά block. Οπότε ναι.

στη shared εχουμε: 24*24 για τον Τ  + 16*16 για τον Υ  +  [2*(24-1)+16] * [2*(24-1)+16] για τον pF.. (16*16 το μεγεθος του block)
συνολο= 4676 * 4 bytes για τους float = 18704..

που κανω λαθος?? κατι απο ολα αυτα μηπως δεν χρειαζεται??

δεν ξέρω πως το υλοποίησες αλλά αν κάθε thread υπολογίζει ενα στοιχείο του πίνακα Y τότε αντι για ένα shared πίνακα Υ μπορείς να ορισεις απλά εναν

Code:

   int result;

που αυτόματα αποθηκεύεται σε register και όταν τελειώσει ο υπολογισμός των γράφεις στη θέση που πρέπει στον Y που έιναι στην global memory.

Επίσης ο buffer για πίνακα pF δεν πρέπει να είναι [2*(24-1)+16] * [2*(24-1)+16]

αλλα [24-1 + 16][24-1 + 16] αφού στον υπολογισμό το κάθε thread χρειάζεται ένα υπο-πίνακα 24χ24 του pF.

Αλλά αυτός είναι διαφορετικό για κάθε thread και έχει ως αρχικό σημειό 0,0 τα i,j του συγκεκριμένου thread.

Άρα για το μέγιστο μέγεθος λογικά είναι το παραπάνω.

Αυτά ισχύουν μόνον αν αντιμετωπίσαμε το πρόβλημα με τον ίδιο τρόπο.

Ερώτηση:

Στην εργασία λέει

συγκρίνει τη μέγιστη διαφορά των 2 αποτελεσμάτων

Τι εννοεί?? να αφαιρέσουμε το πίνακα Υ που υπολογίσαμε με την CPU από αυτόν που

βρήκαμε με την GPU?? Ή αναφέρεται σε χρόνους εκτέλεσης ??

[AgentCain]

To 1o
γιατί όπως μας είπε και ο Πιτς, στην CUDA δεν ασχολήθηκαν τόσο με την ακρίβεια όσο με την ταχύτητα, οπότε υπάρχει μια απόκλιση της ταξης των E-5

[DoomGuard]

To 1o
γιατί όπως μας είπε και ο Πιτς, στην CUDA δεν ασχολήθηκαν τόσο με την ακρίβεια όσο με την ταχύτητα, οπότε υπάρχει μια απόκλιση της ταξης των E-5

Ok thx. 

Ακόμη μία ερώτηση.

Σας είπε ο Πιτσ. ότι οι πίνακες πρέπει να έχουν μέσα float και όχι πχ int??

PS: εγώ δεν έχω δει κάποια απόκλιση πάντως.

[tomshare3]

Επίσης ο buffer για πίνακα pF δεν πρέπει να είναι [2*(24-1)+16] * [2*(24-1)+16]

αλλα [24-1 + 16][24-1 + 16] αφού στον υπολογισμό το κάθε thread χρειάζεται ένα υπο-πίνακα 24χ24 του pF.

Αλλά αυτός είναι διαφορετικό για κάθε thread και έχει ως αρχικό σημειό 0,0 τα i,j του συγκεκριμένου thread.

Άρα για το μέγιστο μέγεθος λογικά είναι το παραπάνω.

Αυτά ισχύουν μόνον αν αντιμετωπίσαμε το πρόβλημα με τον ίδιο τρόπο.

ετσι σκεφτομουν οτι πρεπει να ειναι ο pF (attach) ωστε ολα τα threads να εξυπηρετηθουν σωστα.. εσυ παιρνεις μονο το κομματι στα αριστερα?? ή κατι αλλο?

[AgentCain]

Καινούργιοι χρονοι:

μέγιστο προβλημα 2000χ2000 24χ24
CPU 27.9
GPU 22.4 (χρήση Global)
GPU 9.6 (χρήση Shared)


ερώτηση προς κάθε γνώστη:
Εχει νόημα να αναθέσουμε σε κάθε thread του block να ανεβάσει στη shared ενα κομματάκι του υποπινακα του pF που χρειάζεται?
Προς το παρόν μονο το thread 00 ανεβαζει αλλα σκέφτομαι μήπως δε γίνεται ολα τα thread να γράψουν παραλληλα στη shared, οπότε και να τα βάλω να το κανουν όλα, να μη κερδίσω τπτ γιατί σε κατασταση εγγραφής να τελειώνουν σειριακά

[AgentCain]

Σας είπε ο Πιτσ. ότι οι πίνακες πρέπει να έχουν μέσα float και όχι πχ int??

PS: εγώ δεν έχω δει κάποια απόκλιση πάντως.

Νομίζω μόνο στους float εμφανίζει προβλημα. Εκεί με την μαντισα κτλ

[AgentCain]

ερώτηση προς κάθε γνώστη:
Εχει νόημα να αναθέσουμε σε κάθε thread του block να ανεβάσει στη shared ενα κομματάκι του υποπινακα του pF που χρειάζεται?
Προς το παρόν μονο το thread 00 ανεβαζει αλλα σκέφτομαι μήπως δε γίνεται ολα τα thread να γράψουν παραλληλα στη shared, οπότε και να τα βάλω να το κανουν όλα, να μη κερδίσω τπτ γιατί σε κατασταση εγγραφής να τελειώνουν σειριακά

Όσον αφορα το παραπάνω ειχα μια 2η σκέψη
Ακόμα και αν η shared δεν επιτρέπει παραλληλη εγγραφή, σίγουρα θα γίνει κάποιο παραλληλο fetch από τη global
Οπότε με το να μοιράσω τη δουλεια στα threads, επεσα από το 9.6 στο 1.9

[DoomGuard]

ερώτηση προς κάθε γνώστη:
Εχει νόημα να αναθέσουμε σε κάθε thread του block να ανεβάσει στη shared ενα κομματάκι του υποπινακα του pF που χρειάζεται?
Προς το παρόν μονο το thread 00 ανεβαζει αλλα σκέφτομαι μήπως δε γίνεται ολα τα thread να γράψουν παραλληλα στη shared, οπότε και να τα βάλω να το κανουν όλα, να μη κερδίσω τπτ γιατί σε κατασταση εγγραφής να τελειώνουν σειριακά

Όσον αφορα το παραπάνω ειχα μια 2η σκέψη
Ακόμα και αν η shared δεν επιτρέπει παραλληλη εγγραφή, σίγουρα θα γίνει κάποιο παραλληλο fetch από τη global
Οπότε με το να μοιράσω τη δουλεια στα threads, επεσα από το 9.6 στο 1.9

Ναι τ κάθε thread αν βοηθάει στην αντιγραφή έχει μεγάλη βελτίωση.

Επίσης αν το κάθε thread δεν αντιγράφει συνεχόμενα στοιχειά αλλά
αν πχ αντιγράφει το i στοιχείο βοηθάει το επόμενο να είναι το i + BLOCK_SIZE
για να γίνονται σε bursts οι αντιγραφές όταν τα threads εκτελούνται ταυτόχρονα

Άλλαξα σε float και έχω μια διαφορά 0.000076

Επίσης για κάποιο λόγο με float έχω καλύτερη απόδοση

ετσι σκεφτομουν οτι πρεπει να ειναι ο pF (attach) ωστε ολα τα threads να εξυπηρετηθουν σωστα.. εσυ παιρνεις μονο το κομματι στα αριστερα?? ή κατι αλλο?

Στο attach έχω πως καταλαβαίνω τον buffer για τον pF.

Είναι έτσι επειδη αν δεις την πράξη

Code:

Y[i*(nF+nT-1)+j]+=pF[(i+ii)*(nF+2*nT-2)+j+jj]*T[ii*nT+ jj];

για συγκεκριμένα i και j τα ii και jj αυξάνουν μόνο άρα ο πίνακας θα είναι προς μόνο μία κατεύθυνση.

[Time_Tripper]

για να τρεξουμε το cuda πρεπει να βαλουμε καμια βιβλιοθηκη? Επειδη δεν εχω καρτα με cuda και τρεχω στο Διαδη, δεν μου αναγνωριζει τις εντολες __global__ κλπ. 

[tomshare3]

δεν ξερω αν κανω κατι λαθος αλλα αυτη τη στιγμη ο διαδης δεν τρεχει τιποτα σε cuda... (αν και κανει compile με nvcc!!)

[vasso]

για να τρεξουμε το cuda πρεπει να βαλουμε καμια βιβλιοθηκη? Επειδη δεν εχω καρτα με cuda και τρεχω στο Διαδη, δεν μου αναγνωριζει τις εντολες __global__ κλπ. 

<cuda.h> 
έβαλες;

[DoomGuard]

δεν ξερω αν κανω κατι λαθος αλλα αυτη τη στιγμη ο διαδης δεν τρεχει τιποτα σε cuda... (αν και κανει compile με nvcc!!)

πάτα top και δες αν κάποιος άλλος τρέχει κάτι γιατί η κάρτα γραφικών μόνο μία εργασία τη φορά μπορεί να τρέχει.

για να τρεξουμε το cuda πρεπει να βαλουμε καμια βιβλιοθηκη? Επειδη δεν εχω καρτα με cuda και τρεχω στο Διαδη, δεν μου αναγνωριζει τις εντολες __global__ κλπ. 

Φρόντισε να τα αρχεία σου να τελειώνουν σε .cu και όχι .c επίσης να κάνεις compile με nvcc oxi gcc

[AgentCain]

απο σήμερα το πρωί η κάρτα δε τη παλεύει
Χθες το βράδυ παντως κατα τις 1-2 που τον άφησα, δούλευε κανονικά 

[Time_Tripper]

Αχα! Ωστε εσυ την χαλασες Αλεξ!!!

[AgentCain]

you got me!
της έβαλα να κάνει τόσο intensive δουλειά με το υπεργαμάτο μου πρόγραμμα που κάηκαν οι καταχωρητές της 

μουαχα!