Hello, If your matrices are not too big,
diag((M-C)*(N(:,1:2)-C)') may do the job: it gives the same result as Stéphane's solution HTH Denis [@@ THALES GROUP INTERNAL @@] Unité Mixte de Physique CNRS / THALES 1 Avenue Augustin Fresnel 91767 Palaiseau CEDEx - France Tel : +33 (0)1 69 41 58 52 Fax : +33 (0)1 69 41 58 78 e-mail : denis.cr...@thalesgroup.com<mailto:denis.cr...@thalesgroup.com> <mailto:%20denis.cr...@thalesgroup.com> http://www.trt.thalesgroup.com/ump-cnrs-thales http://www.research.thalesgroup.com 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 : 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
_______________________________________________ users mailing list users@lists.scilab.org http://lists.scilab.org/mailman/listinfo/users