\hB"SSKJr SSKJrJrJrJrJrJr SSK J r SSK J r SSK Jr SSKJr SSKJrJr SSKJrJrJrJrJrJrJrJrJrJrJr SS K J!r! SS K"J#r# SGS jr$S r%S r&Sr'Sr(Sr)Sr*Sr+Sr,Sr-Sr.Sr/Sr0Sr1Sr2Sr3Sr4Sr5Sr6Sr7Sr8S r9S!r:S"r;S#rS&r?S'r@S(rAS)rBS*rCS+rDS,rES-rFS.rGS/rHS0rIS1rJS2rKS3rLS4rMS5rNS6rOS7rPS8rQS9rRS:rSS;rTS<rUS=rVS>rWS?rXS@rYSArZSBr[SCr\SDr]SEr^gF)H)StringIO)symbolsEqpiCatalanLambdaDummyEquality)Symbol)erf)Integral)Matrix MatrixSymbol) codegen make_routineCCodeGen C89CodeGen C99CodeGen InputArgument CodeGenErrorFCodeGenCodeGenArgumentListErrorOutputArgument InOutArgument)raises)implemented_functioncn[5nU"XX#U5 UR5nUR5 U$)aWrapper for dump_fn. dump_fn writes its results to a stream object and this wrapper returns the contents of that stream as a string. This auxiliary function is used by many tests below. The header and the empty lines are not generated to facilitate the testing of the output. )rgetvalueclose)dump_fnroutinesprefixheaderemptyoutputsources Z/var/www/auris/envauris/lib/python3.13/site-packages/sympy/utilities/tests/test_codegen.py get_stringr)s1ZF Hfe4 __ F LLN McZ^^ ^ ^ ^ [S5umm m m T T -T -m [[U U U 4Sj5 [[UU U U U 4Sj5 [S[ TT 5T T TT /S9nUR Vs/sHoR PM snT T TT /:XdeUR Vs/sHn[U5PM sn[[[[/:Xde[S[ T T 5T T T /S9nUR Vs/sHn[U5PM sn[[[/:XdeSSK J nJ n [USS /5upE[S S S 9nU"S U5n[S[ XGXW5XTU/S9nUR Vs/sHoR PM snURURU/:Xde[!T T -T -T SS4T SS45m [S[ TT 5T T TT /S9nUR Vs/sHoR PM snT T TT /:Xdegs snfs snfs snfs snfs snf)Nza x y zc>[STTT/S9$Ntestargument_sequencer)exprxzsr(-test_Routine_argument_order..&s\&$V.%r*c4>[S[TT5TTT/S9$r-)rr)ar2r3yr4sr(r5r6(s"\&"Q C%&1I./r*r.r/r IndexedBaseIdxABmTintegeri)rrrrr argumentsnametyperrr sympy.tensorr;r<maplabelr) rargr;r<r=r>r?rBr8r2r3r9r4s @@@@@r(test_Routine_argument_orderrN#s#JAq!Q E19D #&%& #&/&/0VR4[Q1aLIA!" .#XX .1aA, >> >"#++ /+3T#Y+ /}nm4G GG GVR4[Q1IFA"#++ /+3T#Y+ /}m46 66 6. {S#J 'DAT"A C AVRad^ayIA!" .#XX .177AGGQ2G GG G AaCEAq!9q!Qi 0DVR4[Q1aLIA!" .#XX .1aA, >> >! / / 0 / /s3H#H;H<H#1H(cT[5n[UR/5nUS:Xdeg)Nz$#include "file.h" #include rr)dump_ccode_genr's r(test_empty_c_coderT>s(|H  ,F ? ?? ?r*cn[5n[UR/SS9nUSSS:XdeUSSS:Xdeg)NTr$RzR/****************************************************************************** *a* * * * See http://www.sympy.org/ for more information. * * * * This file is part of 'project' * ******************************************************************************/ #include "file.h" #include rPrRs r(test_empty_c_code_with_commentrYDsV|H D 9F #2;]   #$<["  r*cT[5n[UR/5nUS:Xdeg)Nz9#ifndef PROJECT__FILE__H #define PROJECT__FILE__H #endif )rr)dump_hrRs r(test_empty_c_headerr\Vs(|H  ,F S SS Sr*c[S5upnX-U-n[SU5n[5n[URU/5nSnXg:Xdeg)Nx,y,zr.#include "file.h" #include double test(double x, double y, double z) { double test_result; test_result = z*(x + y); return test_result; } rrrr)rQr3r9r4r2routinerSr'expecteds r(test_simple_c_coderd\sXgGA! E19D64(G|H ' 3F     r*c[S5upnX-U-n[SU5n[5n[URU/5nSnXg:Xdeg)Nzif, typedef, whiler.z#include "file.h" #include double test(double if_, double typedef_, double while_) { double test_result; test_result = while_*(if_ + typedef_); return test_result; } )rrrr)rQras r(test_c_code_reserved_wordsrfnsY*+GA! EQ;D64(G|H ' 3F     r*c[S[[-5n[5n[ UR U/5nS[R "S5-nX#:Xdeg)Nr.z#include "file.h" #include double test() { double test_result; double const Catalan = %s; test_result = pow(M_PI, Catalan); return test_result; } )rrrrr)rQevalfrbrSr'rcs r(test_numbersymbol_c_coderksS62w;/G|H ' 3F   b H   r*c[S5upnX-n[SX2X/S9n[5n[URU/5nSnXg:Xdeg)Nr^r.r/z#include "file.h" #include double test(double z, double x, double y) { double test_result; test_result = x + y; return test_result; } r`ras r(test_c_code_argument_orderrmsXgGA! 5D64q9EG|H ' 3F     r*c[S5upnX-U-n[SU5n[5n[URU/5nSnXg:Xdeg)Nr^r.d#ifndef PROJECT__FILE__H #define PROJECT__FILE__H double test(double x, double y, double z); #endif )rrrr)r[ras r(test_simple_c_headerrpsXgGA! E19D64(G|H ' 3F    r*cb[S5upnX-U-nSS/n[SU4SSSSS9nXT:Xdeg) Nr^)file.cr_)file.hror.CfileFr$r%rr)r3r9r4r2rcresults r(test_simple_c_codegenrysPgGA! E19D  Hfd^S&e LF   r*c^^[S5upnX-U-nX- U-n[SX4/5m[5m[[UU4Sj5 g)Nr^r.c2>[TRT/5$Nr)r[rSrbsr(r5)test_multiple_results_c..HOOgY!Gr*)rrrrrr3r9r4expr1expr2rSrbs @@r(test_multiple_results_crsMgGA! UAIE UAIE G|H <GHr*c&[[S5 g)Nc[S/5$Nr.r1r*r(r5#test_no_results_c.. |FB7r*r ValueErrorrr*r(test_no_results_cr  :78r*cSSKJn SSKJn SSKJnJnJn SSKJ nJ n SSK J n SSK JnJn Jn Jn Jn Jn [)S5nS U"U54S U"U54S U "U54S U "U54S U"U54SU "U54SU"U54SU"U54SU"U54SU"U54SU "U54SU"U54SU"U54SU "U54SU"U54/n[+USSSSS9nUSSS:XdeUSSS:XdeUSSS:XdeUSSS :Xdeg)!NrAbslogcoshsinhtanh)ceilingfloorsqrtacosasinatancossintanr3 test_fabs test_acos test_asin test_atan test_ceiltest_cos test_cosh test_floortest_logtest_lntest_sin test_sinh test_sqrttest_tan test_tanhC89ruFrvrrrCa#include "file.h" #include double test_fabs(double x) { double test_fabs_result; test_fabs_result = fabs(x); return test_fabs_result; } double test_acos(double x) { double test_acos_result; test_acos_result = acos(x); return test_acos_result; } double test_asin(double x) { double test_asin_result; test_asin_result = asin(x); return test_asin_result; } double test_atan(double x) { double test_atan_result; test_atan_result = atan(x); return test_atan_result; } double test_ceil(double x) { double test_ceil_result; test_ceil_result = ceil(x); return test_ceil_result; } double test_cos(double x) { double test_cos_result; test_cos_result = cos(x); return test_cos_result; } double test_cosh(double x) { double test_cosh_result; test_cosh_result = cosh(x); return test_cosh_result; } double test_floor(double x) { double test_floor_result; test_floor_result = floor(x); return test_floor_result; } double test_log(double x) { double test_log_result; test_log_result = log(x); return test_log_result; } double test_ln(double x) { double test_ln_result; test_ln_result = log(x); return test_ln_result; } double test_sin(double x) { double test_sin_result; test_sin_result = sin(x); return test_sin_result; } double test_sinh(double x) { double test_sinh_result; test_sinh_result = sinh(x); return test_sinh_result; } double test_sqrt(double x) { double test_sqrt_result; test_sqrt_result = sqrt(x); return test_sqrt_result; } double test_tan(double x) { double test_tan_result; test_tan_result = tan(x); return test_tan_result; } double test_tanh(double x) { double test_tanh_result; test_tanh_result = tanh(x); return test_tanh_result; } rsa#ifndef PROJECT__FILE__H #define PROJECT__FILE__H double test_fabs(double x); double test_acos(double x); double test_asin(double x); double test_atan(double x); double test_ceil(double x); double test_cos(double x); double test_cosh(double x); double test_floor(double x); double test_log(double x); double test_ln(double x); double test_sin(double x); double test_sinh(double x); double test_sqrt(double x); double test_tan(double x); double test_tanh(double x); #endif )$sympy.functions.elementary.complexesr&sympy.functions.elementary.exponentialr%sympy.functions.elementary.hyperbolicrrr#sympy.functions.elementary.integersrr(sympy.functions.elementary.miscellaneousr(sympy.functions.elementary.trigonometricrrrrrrrr)rrrrrrrrrrrrrrr3 name_exprrxs r(test_ansi_math1_codegenrs_8:HHD=ZZ A c!f d1g d1g d1g gaj! SV d1g uQx SV CF SV d1g d1g SV d1gI"Yve5 IF !9Q<8 ## # !9Q< F  $ !9Q<8 ## # !9Q< 0   r*cSSKJn [S5upSU"X54SX-4/n[USSSSS 9nUSSS :XdeUSS S :XdeUS SS :XdeUS S S:Xdeg)Nratan2x,y test_atan2test_powrruFrvrrrCa*#include "file.h" #include double test_atan2(double x, double y) { double test_atan2_result; test_atan2_result = atan2(x, y); return test_atan2_result; } double test_pow(double x, double y) { double test_pow_result; test_pow_result = pow(x, y); return test_pow_result; } rsz#ifndef PROJECT__FILE__H #define PROJECT__FILE__H double test_atan2(double x, double y); double test_pow(double x, double y); #endif rrrr)rr3r9rrxs r(test_ansi_math2_codegenrs> 5>DA uQ{# QTIYve5 IF !9Q<8 ## # !9Q< N   !9Q<8 ## # !9Q<   r*cSSKJnJnJn [ S5up4nSU"U5U"U5-U"U5-S-R 54SU"U"U"U"U"U"U"U"X4-U-555555554/n[ USSS S S 9nUSSS :XdeUSS S :XdeUS SS:XdeUS S S:Xdeg)Nrrrrr^test1test2rruFrvrrrCa-#include "file.h" #include double test1(double x, double y, double z) { double test1_result; test1_result = pow(sin(x), 7) + 7*pow(sin(x), 6)*cos(y) + 7*pow(sin(x), 6)*tan(z) + 21*pow(sin(x), 5)*pow(cos(y), 2) + 42*pow(sin(x), 5)*cos(y)*tan(z) + 21*pow(sin(x), 5)*pow(tan(z), 2) + 35*pow(sin(x), 4)*pow(cos(y), 3) + 105*pow(sin(x), 4)*pow(cos(y), 2)*tan(z) + 105*pow(sin(x), 4)*cos(y)*pow(tan(z), 2) + 35*pow(sin(x), 4)*pow(tan(z), 3) + 35*pow(sin(x), 3)*pow(cos(y), 4) + 140*pow(sin(x), 3)*pow(cos(y), 3)*tan(z) + 210*pow(sin(x), 3)*pow(cos(y), 2)*pow(tan(z), 2) + 140*pow(sin(x), 3)*cos(y)*pow(tan(z), 3) + 35*pow(sin(x), 3)*pow(tan(z), 4) + 21*pow(sin(x), 2)*pow(cos(y), 5) + 105*pow(sin(x), 2)*pow(cos(y), 4)*tan(z) + 210*pow(sin(x), 2)*pow(cos(y), 3)*pow(tan(z), 2) + 210*pow(sin(x), 2)*pow(cos(y), 2)*pow(tan(z), 3) + 105*pow(sin(x), 2)*cos(y)*pow(tan(z), 4) + 21*pow(sin(x), 2)*pow(tan(z), 5) + 7*sin(x)*pow(cos(y), 6) + 42*sin(x)*pow(cos(y), 5)*tan(z) + 105*sin(x)*pow(cos(y), 4)*pow(tan(z), 2) + 140*sin(x)*pow(cos(y), 3)*pow(tan(z), 3) + 105*sin(x)*pow(cos(y), 2)*pow(tan(z), 4) + 42*sin(x)*cos(y)*pow(tan(z), 5) + 7*sin(x)*pow(tan(z), 6) + pow(cos(y), 7) + 7*pow(cos(y), 6)*tan(z) + 21*pow(cos(y), 5)*pow(tan(z), 2) + 35*pow(cos(y), 4)*pow(tan(z), 3) + 35*pow(cos(y), 3)*pow(tan(z), 4) + 21*pow(cos(y), 2)*pow(tan(z), 5) + 7*cos(y)*pow(tan(z), 6) + pow(tan(z), 7); return test1_result; } double test2(double x, double y, double z) { double test2_result; test2_result = cos(cos(cos(cos(cos(cos(cos(cos(x + y + z)))))))); return test2_result; } rsz#ifndef PROJECT__FILE__H #define PROJECT__FILE__H double test1(double x, double y, double z); double test2(double x, double y, double z); #endif rrrrrexpandr)rrrr3r9r4rrxs r(test_complicated_codegenr,sHHgGA! CFSVOc!f,q088:; #c#c#c#c!%!)n*=&>"?@ABCDIYve5 IF !9Q<8 ## # !9Q<. 0 00 b !9Q<8 ## # !9Q<   r*c SSKJnJn SSKJn U"SSS9up4U"S5nU"S5nU"S 5nU"S U5nU"S U5n [ S [ XxXXU 4Xi-54S SSSS9uupupU S:XdeSnXSSX-U --0-:Xd8XSSXU---0-:Xd&XSSX-U --0-:XdXSSXU---0-:XdeU S:XdeU S:Xdeg)Nrr:rn mTr@r=r3r9rBj matrix_vectorC99ruFrvrra#include "file.h" #include void matrix_vector(double *A, int m, int n, double *x, double *y) { for (int i=0; i>   |qsQw'?@@ @ |qQ3w'?@@ @ |qsQw'?@@ @ |qQ3w'?@@ @B A >>    r*cSSKJnJn [SS[S9up#U"S5nU"S5nU"X#5nSUR R UR S .-n[S [XRXB55n[5n[5n [U RU/5n X:Xde[[5 [URU/5 SSS5 g!,(df  g=f) Nrr:i mTrAclsr3r9z#include "file.h" #include void test_dummies(int m_%(mno)i, double *x, double *y) { for (int i_%(ino)i=0; i_%(ino)i void matrix_vector(double *A, int m, int n, int o, int p, double *x, double *y) { for (int i=o; i<%(upperi)s; i++){ y[i] = 0; } for (int i=o; i<%(upperi)s; i++){ for (int j=0; j> q5+ ,H H|qsQw'?@@ @ H|qQ3w'?@@ @ H|qsQw'?@@ @ H|qQ3w'?@@ @B A >>    r*c SSKJn SSKJnJn [ S5up4n[ SU"XB"U55U"U5/5n[5nURU/SSSS9nUSSS :XdeS n USS U :Xdeg) Nrr rrr^foor.Frvtest.cz#include "test.h" #include double foo(double x, double *y) { (*y) = sin(x); double foo_result; foo_result = cos(x); return foo_result; } rC sympy.core.relationalr rrrrrrwrite r rrr3r9r4rLcrxrcs r(test_output_arg_crs.CgGA!UXaQ0#a&9:A A WWaS&eW double foo(double if_, double *while_) { (*while_) = sin(if_); double foo_result; foo_result = cos(if_); return foo_result; } rCrrs r( test_output_arg_c_reserved_wordsrs.Cn%GA!UXaQ0#a&9:A A WWaS&eW void c(double *A, double *b, double *out) { out[0] = A[0]*b[0] + A[1]*b[1] + A[2]*b[2]; out[1] = A[3]*b[0] + A[4]*b[1] + A[5]*b[2]; out[2] = A[6]*b[0] + A[7]*b[1] + A[8]*b[2]; } ) rrrrbrF result_var_namer)rQ) A_symb_symr=rrcgenrLrrcs r(test_multidim_c_argument_cser s a #E a #Eu Au A A F -D S!A!&AKKO!AKKO dkkA3v 6D     r*c [S5upn[S5up4[SSS5n[U[SSU//55n[X@U-U-5n[USU-5nSXgU/4n Sn [ U S SS S XX!XS4S 9n U S Sn X:Xdeg) Nr^B,Cr=rCrErDr.z#include "test.h" #include void test(double x, double *C, double z, double y, double *A, double *B) { (*C) = z*(x + y); A[0] = 1; A[1] = 2; A[2] = x; (*B) = 2*x; } rFr$r%r0rrrr rr) r3r9r4r>rtr=rrexpr3rrcrxr's r( test_ccode_results_named_orderedrsgGA! 5>DAS!QA QAq {+ ,E QQ "E Q! E%./I  YVE()aA'9;F AYq\F   r*c [SSS5n[SSS5n[SSS5n[SSS5nS[XS SS245[X SSS245[X0SS2S 45/4n[US SS S S 9nUS SnSnXg:Xdeg)Nr=rrEr>rCrtDr.rrDrFrva#include "test.h" #include void test(double *A, double *B, double *C, double *D) { B[0] = A[0]; B[1] = A[1]; B[2] = A[2]; C[0] = A[3]; C[1] = A[4]; C[2] = A[5]; D[0] = A[2]; D[1] = A[5]; D[2] = A[8]; D[3] = A[11]; D[4] = A[14]; } rr rr=r>rtrrrxr'rcs r(test_ccode_matrixsymbol_slicer8sS!QAS!QAS!QAS!QA(11g."11g."11g.01IYve5 IF AYq\F  "   r*c [S5upp#[SSS5nS[U[X-/X-X#--/X-U-U-//55/4n[ SS9n[ XVSSS 9nUS SnS n X:Xdeg) Nza b c derErCr.T)rF)rSr$r%rz#include "test.h" #include void test(double a, double b, double c, double d, double *e) { const double x0 = a*b; const double x1 = c*d; e[0] = x0; e[1] = x0 + x1; e[2] = x0*x1; } )rrr rrr) r8rrdrr generatorrxr'rcs r(test_ccode_cserUs#JA!S!QA(1fqseacACi[13q57)-L&MNOPIT"I Y5 NF AYq\F     r*cr[SSS5nSn[5n[XSSU4S9nUSSnSnXE:Xdeg) Nr3rDrC)r.g?F)rSr$r%r0rz#include "test.h" #include double test(double *x) { double test_result; test_result = 1.0; return test_result; } )rrr)r3rrrxr'rcs r(test_ccode_unused_array_argrisUS!QAI I Y5bcae fF AYq\F     r*c [SSS5n[SSS5n[SSS9nSUS US-4n[US S S XU4S 9nUS SnS nXV:Xdeg)NXrErCYr4Tr@testBugrrtF)languager$r%r0z#include "testBug.h" #include double testBug(double *X, double *Y, int z) { double testBug_result; testBug_result = X[0] + X[1]; return testBug_result; } )rrr)r r!r4rrxr'rcs r( test_ccode_unused_array_arg_funcr$|s~S1AS1Ad#AAaD1Q4K(I YU%\]bc[d eF AYq\F     r*cT[5n[UR/5nUS:XdegNrr)dump_f95rRs r(test_empty_f_coder*s'zH ))2 .F R<[TRT/5$r|r}r~sr(r5)test_multiple_results_f.."rr*)rrrrrrs @@r(test_multiple_results_frDsMgGA! UAIE UAIE GzH <GHr*c&[[S5 g)Nc[S/5$rr1rr*r(r5#test_no_results_f..&rr*rrr*r(test_no_results_frH%rr*cSSKJn SSKJn SSKJnJnJn SSKJ n SSK J nJ nJ nJn Jn Jn [#S5n SU"U 54S U"U 54S U"U 54S U"U 54S U "U 54S U"U 54SU"U 54SU"U 54SU "U 54SU"U 54SU"U 54SU "U 54SU"U 54/ n [%U SSSSS9nUSSS:XdeSnUSSU:XdeUSSS:XdeSnUSSU:Xdeg)Nrrrrrrr3test_absrrrrrrrrrrrrr>ruFrvr?aREAL*8 function test_abs(x) implicit none REAL*8, intent(in) :: x test_abs = abs(x) end function REAL*8 function test_acos(x) implicit none REAL*8, intent(in) :: x test_acos = acos(x) end function REAL*8 function test_asin(x) implicit none REAL*8, intent(in) :: x test_asin = asin(x) end function REAL*8 function test_atan(x) implicit none REAL*8, intent(in) :: x test_atan = atan(x) end function REAL*8 function test_cos(x) implicit none REAL*8, intent(in) :: x test_cos = cos(x) end function REAL*8 function test_cosh(x) implicit none REAL*8, intent(in) :: x test_cosh = cosh(x) end function REAL*8 function test_log(x) implicit none REAL*8, intent(in) :: x test_log = log(x) end function REAL*8 function test_ln(x) implicit none REAL*8, intent(in) :: x test_ln = log(x) end function REAL*8 function test_sin(x) implicit none REAL*8, intent(in) :: x test_sin = sin(x) end function REAL*8 function test_sinh(x) implicit none REAL*8, intent(in) :: x test_sinh = sinh(x) end function REAL*8 function test_sqrt(x) implicit none REAL*8, intent(in) :: x test_sqrt = sqrt(x) end function REAL*8 function test_tan(x) implicit none REAL*8, intent(in) :: x test_tan = tan(x) end function REAL*8 function test_tanh(x) implicit none REAL*8, intent(in) :: x test_tanh = tanh(x) end function rCrsaAinterface REAL*8 function test_abs(x) implicit none REAL*8, intent(in) :: x end function end interface interface REAL*8 function test_acos(x) implicit none REAL*8, intent(in) :: x end function end interface interface REAL*8 function test_asin(x) implicit none REAL*8, intent(in) :: x end function end interface interface REAL*8 function test_atan(x) implicit none REAL*8, intent(in) :: x end function end interface interface REAL*8 function test_cos(x) implicit none REAL*8, intent(in) :: x end function end interface interface REAL*8 function test_cosh(x) implicit none REAL*8, intent(in) :: x end function end interface interface REAL*8 function test_log(x) implicit none REAL*8, intent(in) :: x end function end interface interface REAL*8 function test_ln(x) implicit none REAL*8, intent(in) :: x end function end interface interface REAL*8 function test_sin(x) implicit none REAL*8, intent(in) :: x end function end interface interface REAL*8 function test_sinh(x) implicit none REAL*8, intent(in) :: x end function end interface interface REAL*8 function test_sqrt(x) implicit none REAL*8, intent(in) :: x end function end interface interface REAL*8 function test_tan(x) implicit none REAL*8, intent(in) :: x end function end interface interface REAL*8 function test_tanh(x) implicit none REAL*8, intent(in) :: x end function end interface )rrrrrrrrrrrrrrrrrrr)rrrrrrrrrrrrr3rrxrcs r(test_intrinsic_math_codegenrK)sK8:HH=ZZ A SV d1g d1g d1g SV d1g SV CF SV d1g d1g SV d1gIYve5 IF !9Q<: %% %@  F !9Q<8 ## # !9Q<8 ## #M  ` !9Q<8 ## #r*cSSKJn [S5upSU"X54SX-4/n[USSSSS 9nUSSS :XdeS nUSS U:XdeUS SS :XdeSnUS S U:Xdeg)Nrrrrrr>ruFrvr?zREAL*8 function test_atan2(x, y) implicit none REAL*8, intent(in) :: x REAL*8, intent(in) :: y test_atan2 = atan2(x, y) end function REAL*8 function test_pow(x, y) implicit none REAL*8, intent(in) :: x REAL*8, intent(in) :: y test_pow = x**y end function rCrsainterface REAL*8 function test_atan2(x, y) implicit none REAL*8, intent(in) :: x REAL*8, intent(in) :: y end function end interface interface REAL*8 function test_pow(x, y) implicit none REAL*8, intent(in) :: x REAL*8, intent(in) :: y end function end interface r)rr3r9rrxrcs r(test_intrinsic_math2_codegenrMs> 5>DA uQ{# QTIYve5 IF !9Q<: %% %   !9Q<8 ## # !9Q<8 ## #  !9Q<8 ## #r*cSSKJnJnJn [ S5up4nSU"U5U"U5-U"U5-S-R 54SU"U"U"U"U"U"U"U"X4-U-555555554/n[ USSS S S 9nUSSS :XdeS nUSS U:XdeUS SS:XdeSnUS S U:Xdeg)Nrrr^rrrr>ruFrvr?a{REAL*8 function test1(x, y, z) implicit none REAL*8, intent(in) :: x REAL*8, intent(in) :: y REAL*8, intent(in) :: z test1 = sin(x)**7 + 7*sin(x)**6*cos(y) + 7*sin(x)**6*tan(z) + 21*sin(x) & **5*cos(y)**2 + 42*sin(x)**5*cos(y)*tan(z) + 21*sin(x)**5*tan(z) & **2 + 35*sin(x)**4*cos(y)**3 + 105*sin(x)**4*cos(y)**2*tan(z) + & 105*sin(x)**4*cos(y)*tan(z)**2 + 35*sin(x)**4*tan(z)**3 + 35*sin( & x)**3*cos(y)**4 + 140*sin(x)**3*cos(y)**3*tan(z) + 210*sin(x)**3* & cos(y)**2*tan(z)**2 + 140*sin(x)**3*cos(y)*tan(z)**3 + 35*sin(x) & **3*tan(z)**4 + 21*sin(x)**2*cos(y)**5 + 105*sin(x)**2*cos(y)**4* & tan(z) + 210*sin(x)**2*cos(y)**3*tan(z)**2 + 210*sin(x)**2*cos(y) & **2*tan(z)**3 + 105*sin(x)**2*cos(y)*tan(z)**4 + 21*sin(x)**2*tan & (z)**5 + 7*sin(x)*cos(y)**6 + 42*sin(x)*cos(y)**5*tan(z) + 105* & sin(x)*cos(y)**4*tan(z)**2 + 140*sin(x)*cos(y)**3*tan(z)**3 + 105 & *sin(x)*cos(y)**2*tan(z)**4 + 42*sin(x)*cos(y)*tan(z)**5 + 7*sin( & x)*tan(z)**6 + cos(y)**7 + 7*cos(y)**6*tan(z) + 21*cos(y)**5*tan( & z)**2 + 35*cos(y)**4*tan(z)**3 + 35*cos(y)**3*tan(z)**4 + 21*cos( & y)**2*tan(z)**5 + 7*cos(y)*tan(z)**6 + tan(z)**7 end function REAL*8 function test2(x, y, z) implicit none REAL*8, intent(in) :: x REAL*8, intent(in) :: y REAL*8, intent(in) :: z test2 = cos(cos(cos(cos(cos(cos(cos(cos(x + y + z)))))))) end function rCrsa4interface REAL*8 function test1(x, y, z) implicit none REAL*8, intent(in) :: x REAL*8, intent(in) :: y REAL*8, intent(in) :: z end function end interface interface REAL*8 function test2(x, y, z) implicit none REAL*8, intent(in) :: x REAL*8, intent(in) :: y REAL*8, intent(in) :: z end function end interface r) rrrr3r9r4rrxrcs r(test_complicated_codegen_f95rO sHHgGA! CFSVOc!f,q088:; #c#c#c#c!%!)n*=&>"?@ABCDIYve5 IF !9Q<: %% %  < !9Q<8 ## # !9Q<8 ## #  $ !9Q<8 ## #r*c *SSKJnJn SSKJn U"SSS9up4[ US5upVnU"SU5nU"S U5n [ S [XxXXU 4Xi-54S S S S S9uupupU S:XdeSnXSS0-:Xd XSS0-:XdeU S:XdeU S:Xdeg)Nrr:rzn,mTr@AxyrBrrr>ruFrvr?a~subroutine matrix_vector(A, m, n, x, y) implicit none INTEGER*4, intent(in) :: m INTEGER*4, intent(in) :: n REAL*8, intent(in), dimension(1:m, 1:n) :: A REAL*8, intent(in), dimension(1:n) :: x REAL*8, intent(out), dimension(1:m) :: y INTEGER*4 :: i INTEGER*4 :: j do i = 1, m y(i) = 0 end do do i = 1, m do j = 1, n y(i) = %(rhs)s + y(i) end do end do end subroutine r A(i, j)*x(j) x(j)*A(i, j)rsainterface subroutine matrix_vector(A, m, n, x, y) implicit none INTEGER*4, intent(in) :: m INTEGER*4, intent(in) :: n REAL*8, intent(in), dimension(1:m, 1:n) :: A REAL*8, intent(in), dimension(1:n) :: x REAL*8, intent(out), dimension(1:m) :: y end subroutine end interface )rIr;r<rrrJrrrs r( test_loopsrTGs-) 5$ 'DA+u%GA! C A C A") "QT1T714<015&V[#]JR     * un55 5 E>2223 3 >>     r*c6SSKJnJn [SS[S9up#U"S5nU"S5nU"X#5nSUR R UR S .-n[S [XRXB55n[5n[URU/5n X:Xdeg) Nrr:rTrr3r9a>subroutine test_dummies(m_%(mcount)i, x, y) implicit none INTEGER*4, intent(in) :: m_%(mcount)i REAL*8, intent(in), dimension(1:m_%(mcount)i) :: x REAL*8, intent(out), dimension(1:m_%(mcount)i) :: y INTEGER*4 :: i_%(icount)i do i_%(icount)i = 1, m_%(mcount)i y(i_%(icount)i) = x(i_%(icount)i) end do end subroutine )icountmcountr) rIr;r<rr rKrrrrr)r)) r;r<rBr?r3r9rcrLrrs r(test_dummy_loops_f95rXzs- 5$E 2DACACA A A 77&&!--@ AH ^Rad^4A A ajj1# &D   r*c hSSKJnJn SSKJn U"SSS9up4pVU"S5upxn U"U5U"X65U"XE54nU"U5U"XE5nU"U 5U"X65n [ S[ XXx--54S S S S S 9uupupU S :XdeSnXSS0-:Xd XSS0-:XdeU S:XdeU S:Xdeg)Nrr:rzi,j,n,mTr@zA,x,yrr>ruFrvr?aasubroutine matrix_vector(A, m, n, x, y) implicit none INTEGER*4, intent(in) :: m INTEGER*4, intent(in) :: n REAL*8, intent(in), dimension(1:m, 1:n) :: A REAL*8, intent(in), dimension(1:n) :: x REAL*8, intent(inout), dimension(1:m) :: y INTEGER*4 :: i INTEGER*4 :: j do i = 1, m do j = 1, n y(i) = %(rhs)s + y(i) end do end do end subroutine rrRrSrsainterface subroutine matrix_vector(A, m, n, x, y) implicit none INTEGER*4, intent(in) :: m INTEGER*4, intent(in) :: n REAL*8, intent(in), dimension(1:m, 1:n) :: A REAL*8, intent(in), dimension(1:n) :: x REAL*8, intent(inout), dimension(1:m) :: y end subroutine end interface r)r;r<rrBrrr?r=r3r9rrrrrcs r(test_loops_InOutrZs-)D1JA!gGA!As1y#a)+,AAs1y!AAs1y!A") "QAC.)5&e#UJR     $ ~66 6 ~66 68 7 >>     r*c zSSKJnJn SSKJn U"SSS9up4pVU"SXF4S9nU"S 5nU"S 5n U"S XTS - 45n U"S U5n [ S[ XXzU 4X-54SSSSS9uupupSS[US- 5[SU-5[US- U- 5S.-nU USS0-:XdU USS0-:Xdegg)Nrr:rrTr@r=rr3r9rBrrrr>ruFrvasubroutine matrix_vector(A, m, n, o, p, x, y) implicit none INTEGER*4, intent(in) :: m INTEGER*4, intent(in) :: n INTEGER*4, intent(in) :: o INTEGER*4, intent(in) :: p REAL*8, intent(in), dimension(1:m, 1:p) :: A REAL*8, intent(in), dimension(1:n) :: x REAL*8, intent(out), dimension(1:%(iup-ilow)s) :: y INTEGER*4 :: i INTEGER*4 :: j do i = %(ilow)s, %(iup)s y(i) = 0 end do do i = %(ilow)s, %(iup)s do j = 1, n y(i) = %(rhs)s + y(i) end do end do end subroutine rrrC)riupilowziup-ilowrrRrS)rIr;r<rrrrstrrs r(test_partial_loops_fr_s.)D1JA!Cv&ACACA C!UA C A") "QT1T714<015&V[#]JR *1q5zAE A N  +H8 8un55 5 E>2223 32 6r*c SSKJn SSKJnJn [ S5up4n[ SU"XB"U55U"U5/5n[5nURU/SSSS9nUSSS :XdeUSS S :Xdeg) Nrr rr^rr.Frvztest.f90rCz~REAL*8 function foo(x, y) implicit none REAL*8, intent(in) :: x REAL*8, intent(out) :: y y = sin(x) foo = cos(x) end function ) rr rrrrrrr) r rrr3r9r4rLrrxs r(test_output_arg_fras.CgGA!UXaQ0#a&9:A A WWaS&eW ?G ajj7) ,D   "D d? " 4":%%& &% #r*c[S5nSnUHnX- nM [SU5n[5n[URU/5nSnXV:Xdeg)Nzx:20rr.a REAL*8 function test(x0, x1, x10, x11, x12, x13, x14, x15, x16, x17, x18, & x19, x2, x3, x4, x5, x6, x7, x8, x9) implicit none REAL*8, intent(in) :: x0 REAL*8, intent(in) :: x1 REAL*8, intent(in) :: x10 REAL*8, intent(in) :: x11 REAL*8, intent(in) :: x12 REAL*8, intent(in) :: x13 REAL*8, intent(in) :: x14 REAL*8, intent(in) :: x15 REAL*8, intent(in) :: x16 REAL*8, intent(in) :: x17 REAL*8, intent(in) :: x18 REAL*8, intent(in) :: x19 REAL*8, intent(in) :: x2 REAL*8, intent(in) :: x3 REAL*8, intent(in) :: x4 REAL*8, intent(in) :: x5 REAL*8, intent(in) :: x6 REAL*8, intent(in) :: x7 REAL*8, intent(in) :: x8 REAL*8, intent(in) :: x9 test = x0 + x1 + x10 + x11 + x12 + x13 + x14 + x15 + x16 + x17 + x18 + & x19 + x2 + x3 + x4 + x5 + x6 + x7 + x8 + x9 end function r1)r3r2symrbrSr'rcs r(test_f_code_call_signature_wrapri s]A D 64(GzH ))G9 5FH8   r*cN^^[S5umm[[UU4Sj5 g)Nzx,Xc(>[STT-4SS5$)Nr.f95r#)rr r3sr(r5!test_check_case..Js&!A#x!Hr*)rrrrms@@r(test_check_caseroHs 5>DAq <HIr*c[S5n[SSS9n[SX-4SS5 g![a8nURSR S5(a [ S 5eSnAgSnAff=f) Nr3T) my_assumptionr.rlr#rzFortran ignores case.z$This exception should not be raised!)rrrre startswithAssertionError)x1x2rs r(test_check_case_false_positivervMsn B D )BI1 I 66!9   7 8 8 !GH H 9Is) A+.A&&A+ct[S5upSSU-4/n[USSSS9n[USSSSS9nUSSUSS:XdeSX-4n[USSSS9n[USSSSS9nUSSUSS:XdeS[X/X-X- //54n[US SSS9n[US SSSS9nUSSUSS:Xdeg) NrrrDr>FrvrrCr)rrr)r3r9rrx expresults r( test_c_fortran_omit_routine_namery[s 5>DA1I Ye5 AF 5%UKI !9Q<9Q<? ** * I Ye5 AF 5%UKI !9Q<9Q<? ** *ac 345I Ye5 AF 5%UKI !9Q<9Q<? ** *r*c[S5upnX-n[X/US//5nSX444n[USSSSS9nUSSnS nUR5S n U R S 5n U Sn US U 0-nXx:Xdeg) Nr^r.rlFrvrrCa2REAL*8 function test(x, y, z, out_%(hash)s) implicit none REAL*8, intent(in) :: x REAL*8, intent(in) :: y REAL*8, intent(in) :: z REAL*8, intent(out), dimension(1:2, 1:2) :: out_%(hash)s out_%(hash)s(1, 1) = x out_%(hash)s(2, 1) = z out_%(hash)s(1, 2) = y out_%(hash)s(2, 2) = 16 test = x + y end function r_hash)rrr splitlinessplit) r3r9r4e1e2rrxr'rcr8rrs r(test_fcode_matrix_outputrmsgGA! B !!R! "B""I Yve5 IF AYq\F   AA  A A$C63-'H   r*c [S5upn[S5up4[SSS5n[U[SSU//55n[X@U-U-5n[USU-5nSXgU/4n [ U SSS S XXXS4S 9n U S Sn S n X:Xdeg) Nr^rr=rCrErDr.rlFrra!subroutine test(x, z, y, C, A, B) implicit none REAL*8, intent(in) :: x REAL*8, intent(in) :: z REAL*8, intent(in) :: y REAL*8, intent(out) :: C REAL*8, intent(out) :: B REAL*8, intent(out), dimension(1:1, 1:3) :: A C = z*(x + y) A(1, 1) = 1 A(1, 2) = 2 A(1, 3) = x B = 2*x end subroutine r) r3r9r4r>rtr=rrrrrxr'rcs r( test_fcode_results_named_orderedrsgGA! 5>DAS!QA QAq {+ ,E QQ "E Q! E%./I Yve5()aA'9;F AYq\F    r*c [SSS5n[SSS5n[SSS5n[SSS5nS[XS SS245[X SSS245[X0SS2S45/4n[US SS S S 9nUS SnS nXg:Xdeg)Nr=rDrEr>rCrtrr.rrlFrvasubroutine test(A, B, C, D) implicit none REAL*8, intent(in), dimension(1:2, 1:3) :: A REAL*8, intent(out), dimension(1:1, 1:3) :: B REAL*8, intent(out), dimension(1:1, 1:3) :: C REAL*8, intent(out), dimension(1:2, 1:1) :: D B(1, 1) = A(1, 1) B(1, 2) = A(1, 2) B(1, 3) = A(1, 3) C(1, 1) = A(2, 1) C(1, 2) = A(2, 2) C(1, 3) = A(2, 3) D(1, 1) = A(1, 3) D(2, 1) = A(2, 3) end subroutine rrs r(test_fcode_matrixsymbol_slicersS!QAS!QAS!QAS!QA(11g."11g."11g.01IYve5 IF AYq\F  "   r*c[SSS5nSUSS2S44n[USSSSS9nUS SnS nUR5SnURS 5nUSnUS U0-nX4:Xdeg) Nr=rDrEr.rCrlFrvrzsubroutine test(A, out_%(hash)s) implicit none REAL*8, intent(in), dimension(1:2, 1:3) :: A REAL*8, intent(out), dimension(1:2, 1:1) :: out_%(hash)s out_%(hash)s(1, 1) = A(1, 2) out_%(hash)s(2, 1) = A(2, 2) end subroutine r|r})rrr~r)r=rrxr'rcr8rrs r(&test_fcode_matrixsymbol_slice_autonamersS!QA1a4!I Yve5 IF AYq\F   AA  A A$C63-'H   r*c[S5upp#[SX-4SSSU4S9nUSSnSnXV:XdeS n[SX-U-4S SSX#4S9nUSSnXV:Xdeg) Nzx y z tfr>F)r$r% global_varsrrCzPREAL*8 function f(x) implicit none REAL*8, intent(in) :: x f = x*y end function z#include "f.h" #include double f(double x, double y) { double f_result; f_result = x*y + z; return f_result; } rtrw)r3r9r4trxr'rcs r(test_global_varsrs#JA! c13ZuE"#'F AYq\F       c13q5\3uE"#)F AYq\F   r*c8SSKJn SSKJn U"SSS00S9n[ S5up4U"X4-5n[ US /S 9nS n[ S U4S S US9nUSSn X:Xde[ US9nURRS 5 Sn[ S U4S S US9nUSSn X:Xdeg)NrC99CodePrinter)expuser_functionsrfastexp)settingszx yz#include "fastexp.h")printerpreprocessor_statementsz#include "expr.h" #include "fastexp.h" double expr(double x, double y) { double expr_result; expr_result = fastexp(x + y); return expr_result; } r2F)r$r%rSrC)rz#include "expr.h" #include #include "fastexp.h" double expr(double x, double y) { double expr_result; expr_result = fastexp(x + y); return expr_result; } ) sympy.printing.crrrrrrrappend) rrrr3r9r2genrcrxr's r(test_custom_codegenrs/:'7%9K&LMG 5>DA qu:D W.D-E GC  fd^E MF AYq\F    W %C&&'=>  fd^E MF AYq\F   r*c SSKJn "SSU5n[S5nUS-nSS/n[S U4S S S S U"5S 9nXT:Xdeg)Nrrc\rSrSrSrSrg)*test_c_with_printer..CustomPrinteri*cSRURUR5URUR55$)Nzfastpow({}, {}))format_printbaser)selfr2s r( _print_Pow5test_c_with_printer..CustomPrinter._print_Pow+s6$++DKK ,B,0KK,AC Cr*rN)__name__ __module__ __qualname____firstlineno__r__static_attributes__rr*r( CustomPrinterr*s Cr*rr3rE)rrz#include "file.h" #include double test(double x) { double test_result; test_result = fastpow(x, 3); return test_result; } )rszP#ifndef PROJECT__FILE__H #define PROJECT__FILE__H double test(double x); #endif r.rtruF)r$r%r)rrrr)rrr3r2rcrxs r(test_c_with_printerr'sb/CC  A a4D  Hfd^SUVcVe fF   r*cRSSKnSURRl[ SSS9n[ SSS9n[SX-4SSSSS 9nUSS nS nXE:Xde[ S5n[ SSS9n[SX-4SSSSS 9nUSS nS nXE:XdeSURRlg) NrTr3)realr9r.rlFrvrCzsREAL*8 function test(x, y) implicit none REAL*8, intent(in) :: x REAL*8, intent(in) :: y test = x + y end function z{COMPLEX*16 function test(x, y) implicit none COMPLEX*16, intent(in) :: x REAL*8, intent(in) :: y test = x + y end function )sympy.utilities.codegen utilitiesrCOMPLEX_ALLOWEDr )sympyr3r9rxr'rcs r(test_fcode_complexrDs".2EOO+sAsA fQS\5&e LFQilF     s AsA fQS\5&e LFQilF     .3EOO+r*N)ruFF)_ior sympy.corerrrrrr rr rr 'sympy.functions.special.error_functionsr sympy.integrals.integralsrsympy.matricesrrrrrrrrrrrrrrsympy.testing.pytestrsympy.utilities.lambdifyrr)rNrTrYr\rdrfrkrmrpryrrrrrrrrrrr rrrrr$r*r,r.r2r4r7r9r<r@rDrHrKrMrOrTrXrZr_rarfrirorvryrrrrrrrrrr*r(rst>>*$7./(9 ?6@ $T $$"$ , I99x0AH'T$.)X$*$*.4:(&&   $  &$&6 I9o$d+$\;$|0f0.b,3^&&6%PJ I+$:<<.8,\:4r*