loop_eecoll_dke_yp

LUKE - Electron-electron Belaiev-Budker collision operator (loop function for distributed computing environment)

Electron-electron Belaiev-Budker collision operator (loop function
for distributed computing environment)

by Yves Peysson (CEA-IRFM,yves.peysson@cea.fr) and Joan Decker (CEA-IRFM,joan.decker@cea.fr)

0001 function [F1_1,F1_2,F2_1,B1_1,B1_2,B1_3,B1_4,B1_5,B2_1,B2_2,B2_3] = loop_eecoll_dke_yp(ip,pn,sfM,spn,spn2,sv,sgamma,sgamma2,ssigma,sz,sJ1,sJ2)
0002 %LUKE - Electron-electron Belaiev-Budker collision operator (loop function for distributed computing environment)
0003 %
0004 %Electron-electron Belaiev-Budker collision operator (loop function
0005 %for distributed computing environment)
0006 %
0007 %by Yves Peysson (CEA-IRFM,yves.peysson@cea.fr) and Joan Decker (CEA-IRFM,joan.decker@cea.fr)
0008 %
0009     smasku = spn <= pn(ip);
0010     smaskl = spn >= pn(ip);
0011     %
0012     F1_1 = trapz_dke_yp([spn,sfM.*spn.*sv.*smasku])';
0013     F1_2 = trapz_dke_yp([spn,sfM.*spn.*sv.*(1 - sgamma.*ssigma./sz).*smasku])';
0014     F2_1 = trapz_dke_yp([spn,sfM.*spn.*smaskl])';
0015     %
0016     B1_1 = trapz_dke_yp([spn,sfM.*smasku.*spn2])';
0017     B1_2 = trapz_dke_yp([spn,sfM.*smasku.*spn2.*spn2])';
0018     B1_3 = trapz_dke_yp([spn,sfM.*smasku.*spn2.*sJ1./sgamma])';
0019     B1_4 = trapz_dke_yp([spn,sfM.*smasku.*spn2.*sJ2./sgamma])';
0020     B1_5 = trapz_dke_yp([spn,sfM.*smasku.*spn2.*((sgamma - ssigma./sz)./sgamma)])';
0021     B2_1 = trapz_dke_yp([spn,sfM.*smaskl.*spn2./sv])';
0022     B2_2 = trapz_dke_yp([spn,sfM.*smaskl.*spn2./sv./sgamma2])';
0023     B2_3 = trapz_dke_yp([spn,sfM.*smaskl.*spn2.*sv])'; 
0024 end