Discussion :

[Aytonee]

Bonjour,

Je suis en train de bosser sur un petit algo et je me heurte un petit soucis (je suis limité en mémoire) et je cherche à regarder si en faisant les calculs en simple précision, la résolution de mes résultats sera toujours bonne.
J'avais développé un mex qui faisait ça pour moi, mais là je bosse sur une autre fonction qui est assez compliquée à passer en mex, du coup je cherchais à voir comment forcer le calcul en simple précision (notamment à l'intérieur des fonctions prédéfinies).

J'ai déjà regardé en passant les données d'entrée en single et en faisant tourner l'algo, mais à bien regarder, tous les calculs se font en double.

Une idée de comment je pourrai procéder ?

Merci d'avance

[Jerome Briot]

Sans plus d'information, je pense malheureusement que seuls les MEX peuvent répondre rapidement à ton besoin.

Ou alors il faudrait que tu nous montres ton code pour que nous puissions t'aider à l'optimiser.

Si tu possèdes la Parallel Computing Toolbox, et si l'algorithme répond aux contraintes requises, tu peux peut être aussi envisager de subdiviser ton problème en sous-problème moins gourmands ?

Bref, donne nous plus d'informations… et de surtout de code

[Aytonee]

Alors en fait j'ai quelques contraintes justement et l'une d'elles est de devoir diminuer substantiellement l'utilisation RAM et les temps de calcul pour pour, à terme, implémenter ça en C.

Pour donner un exemple de ce que je cherche à faire : J'ai un traitement qui consiste à interpoler mon signal sur une base de temps différente en utilisant simplement la fonction spline :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

samplingData = spline(timeData,Data,newTime);

J'ai joué sur le type en modifiant les entrées avant :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
 
Data = single(Data);
timeData = single(timeData);
newTime = single(newTime);
samplingData = spline(timeData,Data,newTime);

En procédant comme ça je ne fais rentrer que des variables en simple précision dans la fonction, et donc il me sort un output en simple précision, mais en regardant ce qui se passe à l'intérieur, j'ai constaté que dès le premier calcul intermédiaire, il repasse en double précision et ce jusqu'à la sortie qu'il converti en simple précision.

Mon idée serait de dupliquer la fonction et de la modifier pour que les calculs soient faits en simple précision (intégralement ou en en laissant certains en double pour conserver de la marge).

Est ce que c'est un peu plus clair comme ça ?

[Jerome Briot]

Que renvoi ceci ?

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

whos timeData Data newTime

Sinon, as-tu cherché sur le File Exchange pour des contributions d'interpolation faites avec des fichiers MEX ?
Il serait sans doute plus simple d'y faire les modifications…