|
Title: κλήση ρουτίνας σε C μέσα από Matlab Post by: snake5 on February 12, 2007, 22:45:07 pm γεια χαρά,
έχει προσπαθήσει κανενας να καλέσει ρουτίνα σε C μέσα απ το Matlab? Το πρόβλημά μου είναι ότι γράφω ένα πρόγραμμα σε C,η κλήση του γίνεται κανονικά από περιβάλλον Matlab,μου επιστρέφει κανονικά τις εξόδους που θέλω αλλά αργεί αρκετά. Ο κώδικάς μου αποτελείται από τρία nested for loops που τρέχουν το εξωτερικό για κάποιες χιλιάδες φορές,το μεσαίο γύρω στις 200 και το εσωτερικό γύρω στις 500.Ενδιάμεσα γίνονται κάποιοι υπολογισμοί και συγκρίσεις με double τιμές. Εξαρχής βάζω την C στο παιχνίδι γιατί πολύ απλά αυτό που θέλω να κάνω δεν μπορεί να τρέξει καν στο Matlab (κολλάει ανεπανόρθωτα..).Με προβληματίζει όμως ότι ο αντίστοιχος κώδικας σε Java πχ νομίζω ότι θα έτρεχε άνετα και χωρίς καμία καθυστέρηση.Έχει κανείς καμιά ιδέα τί μπορεί να γίνεται με τη χρήση της μνήμης σε αυτή την περίπτωση?Απ'ότι διάβασα,υπάρχει ένα θέμα χειρισμού της heap περιοχής της μνήμης και της stack όταν πας να συνδέσεις C κώδικα με το interface του Matlab.Μου φαίνεται ότι το τελευταίο είναι που κάνει όλη τη ζημιά γιατί ο απλός κώδικας σε C δε κολώνει με for loops... Thx Title: Re: κλήση ρουτίνας σε C μέσα από Matlab Post by: Netgull on February 13, 2007, 01:45:08 am Τα loops σου είναι στο Matlab ή στη C;
Αν είναι στο Matlab φρόντισε να αρχικοποιείς όλους τους πίνακες που χρησιμοποιείς πριν το loop. Title: Re: κλήση ρουτίνας σε C μέσα από Matlab Post by: snake5 on February 13, 2007, 02:30:00 am Όχι,όλα σε C εκτελούνται
Title: Re: κλήση ρουτίνας σε C μέσα από Matlab Post by: TheDoctor on February 17, 2007, 13:10:59 pm Ελπίζω να χρειάζεσαι ακόμα αυτή τη λύση γιατί είναι μακροσκελής και θα είναι κρίμα να την γράφω τζάμπα...
Για να τρέξεις από το matlab κώδικα C χρειάζεται να βάλεις τον κώδικά σου σε ένα DLL. Το Dll αυτό θα πρέπει να κάνει export μια μόνο συναρτηση που θα λέγεται mexFunction και θα συντάσσεται ως εξής: void mexFunction(int nlhs, mxArray *plhs[], int nrhs,const mxArray *prhs[]); η παράμετρος nlhs σημαίνει number of left hand side arguments, και είναι το πόσα αποτελέσματα περιμένει το matlab από εσένα. Τα αποτελέσματα θα πρέπει να τα βάλεις στον πίνακα *plhs[] όπως θα φανεί στο παρακάτω παράδειγμα. Η παράμετρος nrhs σημαίνει number of right hand side arguments και είναι το πόσες παραμέτρους έδωσε ο χρήστης. Οι παράμετροι αυτοί υπάρχουν στον πίνακα *prhs[] Εσύ τώρα πρέπει να δείς άν τα δεδομένα είναι όσα θέλεις και του τύπου που θέλεις, να δείς αν τα αποτελέσματα είναι αυτά που θέλεις, και μετά να κάνεις τις πράξεις σου. Για να αποκτήσεις πρόσβαση στα δεδομένα που σου έρχονται από το matlab χρησιμοποιείς τη συνάρτηση mxGetPr που σου δίνει έναν double pointer στα δεδομένα. Για να δείς τις διαστάσεις του πίνακα χρησιμοποιείς τις συναρτήσεις mxGetΜ και mxGetN. Ένα παράδειγμα συνάρτησης φαίνεται παρακάτω : void mexFunction(int nlhs,mxArray *plhs[],int nrhs,const mxArray *prhs[]) { int i; /* Βλέπει πόσες παράμετρους έδωσε ο χρήστης. */ mexPrintf("\n Yparxoun %d parametroi ", nrhs); /* Βλέπει πόσα αποτελέσματα ζήτησε ο χρήστης */ mexPrintf("\n\nYparxoun %d apotelesmata\n", nlhs); for (i=0; i<nlhs; i++) { plhs=mxCreateDoubleMatrix(1,1,mxREAL); *mxGetPr(plhs)=42; /* Βάζει σε όλες τις εξόδους το 42 (Douglas Adams, Hitchhiker's Guide to the galaxy, Η απάντηση στο πιο σύνθετο ερώτημα του σύμπαντος*/ } } στο C αρχείο κάνεις include to "mex.h" που σου δίνει πρόσβαση στις συναρτήσεις αυτές. Το compile πρέπει να γίνει μέσω της εντολής mex του matlab η οποία κανονίζει τα σύνθετα πράγματα του να κάνεις dll exports και λοιπά γράφοντας απλά mex myFile.c enjoy! Κώστας Πανούλας Title: Re: κλήση ρουτίνας σε C μέσα από Matlab Post by: fugiFOX on February 17, 2007, 13:48:15 pm Μην ανησυχείς γιατρέ.
Ακόμη και εάν το Φίδι5 έλυσε το πρόβλημά του, σίγουρα το ποστ σου δεν θα πάει χαμένο γιατί θα το χρησιμοποιήσει κάποιος άλλος που θα έχει το συγκεκριμένο πρόβλημα. Αυτή ειναι η δύναμη του γραπτου (και δημοσιευμένου) λόγου Title: Re: κλήση ρουτίνας σε C μέσα από Matlab Post by: snake5 on February 17, 2007, 16:58:23 pm Πανούλα,όντως μπράβο,καλή δουλειά αλλά τη δουλειά την έκανα...
είχα υπερεκτιμήσει τις δυνατότητες της C.Έδιωξα ένα for και από τριπλό το έκανα διπλό και πάει σφαίρα.Πάντως νομίζω είναι πολύ χρήσιμη φάση τα mex files γιατί άμα θες να κάνεις χοντρές δουλειές και θες να είσαι σε περιβάλλον Matlab,είναι ο μόνος τρόπος να ξεφύγεις από τον interpreted χαρακτήρα του τελευταίου και να μην περιμένεις με τις μέρες για να δεις έξοδο από τη δυνάρτησή σου... Μπράβο once more |