Thanks Stéphane Works fine and fast (I'm speaking about my application :-) )
Paul EXPORT CONTROL : Cet email ne contient pas de données techniques This email does not contain technical data De : users [mailto:users-boun...@lists.scilab.org] De la part de Stéphane Mottelet Envoyé : mercredi 7 novembre 2018 13:40 À : users@lists.scilab.org Objet : [EXTERNAL] Re: [Scilab-users] function and vectorization Hello, Try: function [Scar_P]=Scalar_product(C, N, M) // Scalar product Scar_P = (M(:,1) - C(:,1)).*(N(:,1) - C(:,1)) + (M(:,2) - C(:,2)).*(N(:,2) - C(:,2)); //printf("C = (%g,%g)\n",C(1),C(2)); printf("N = (%g,%g)\n",N(1),N(2)); printf("M = (%g,%g)\n",M(1),M(2)); endfunction with (do not transpose N(i,1:2)) i = 1 : n; Scal(i) = Scalar_product(C(i,:),N(i,1:2),M(i,:)) S. Le 07/11/2018 à 12:35, Carrico, Paul a écrit : ##################################################################### mode(0) clear function [Scar_P]=Scalar_product(C, N, M) // Scalar product Scar_P = (M(1) - C(1))*(N(1) - C(1)) + (M(2) - C(2))*(N(2) - C(2)); printf("C = (%g,%g)\n",C(1),C(2)); printf("N = (%g,%g)\n",N(1),N(2)); printf("M = (%g,%g)\n",M(1),M(2)); endfunction n = 10; C = rand(n,2); M = rand(n,2); N = rand(n,5); Scal = zeros(n); Scal2 = Scal printf(" ****************\n"); i = 1 : n; Scal(i) = Scalar_product(C(i,:),N(i,1:2)',M(i,:)) printf("\n ****************\n"); for i = 1 : n Scal2(i) = Scalar_product(C(i,:),N(i,1:2)',M(i,:)); end -- Stéphane Mottelet Ingénieur de recherche EA 4297 Transformations Intégrées de la Matière Renouvelable Département Génie des Procédés Industriels Sorbonne Universités - Université de Technologie de Compiègne CS 60319, 60203 Compiègne cedex Tel : +33(0)344234688 http://www.utc.fr/~mottelet<https://urldefense.proofpoint.com/v2/url?u=http-3A__www.utc.fr_-7Emottelet&d=DwMG-g&c=0hKVUfnuoBozYN8UvxPA-w&r=4TCz--8bXfJhZZvIxJAemAJyz7Vfx78XvgYu3LN7eLo&m=4nq72IuY_O4BD4bc-E25iQkCiD8-ZEBjCl6xzxjHbWw&s=9-ZrmcA1S1rslVYW5SkH4zF_TANcPxDeWItXi95Gqls&e=>
_______________________________________________ users mailing list users@lists.scilab.org http://lists.scilab.org/mailman/listinfo/users