o 8Va>@sddlmZmZmZmZmZmZmZmZm Z m Z m Z m Z m Z mZmZmZmZmZmZmZmZmZmZmZmZmZmZmZmZmZmZddl m!Z!m"Z"m#Z#ddl$m%Z%m&Z&m'Z'm(Z(m)Z)m*Z*m+Z+m,Z,m-Z-m.Z.m/Z/m0Z0ddl1m2Z2ddl3m4Z4m5Z5m6Z6m7Z7m8Z8ddl9m:Z:m;Z;ddlZ>ddl?m@Z@mAZAdd lBmCZCdd lDmEZEmFZFd d ZGd dZHddZIddZJddZKddZLddZMddZNddZOddZPdd ZQd!d"ZRd#d$ZSd%d&ZTd'd(ZUd)d*ZVd+d,ZWd-d.ZXd/d0ZYd1d2ZZd3d4Z[d5d6Z\d7d8Z]d9d:Z^d;d<Z_d=d>Z`d?d@ZadAdBZbdCdDZcdEdFZddGdHZedIdJZfdKdLZgdMdNZhdOdPZidQdRZjdSdTZkdUdVZldWdXZmdYdZZnd[d\Zod]d^Zpd_S)`)sincosatan2logexpgamma conjugatesqrt factorialIntegral PiecewiseAdddiffsymbolsSFloatDummyEqRangeCatalan EulerGammaE GoldenRatioIpiFunctionRationalIntegerLambdasignMod)For Assignment aug_assign) DeclarationVariablefloat32float64 value_constrealbool_WhileFunctionPrototypeFunctionDefinitionintegerReturn) Relational)AndOrNot EquivalentXor)Matrix MatrixSymbol)fcode FCodePrinter) IndexedBaseIdx)implemented_function)raiseswarns_deprecated_sympycCs0td}Gdddt}t||dksJdS)Nxc@seZdZddZdS)ztest_printmethod..nintcSsd||jdS)Nznint(%s)r)_printargs)selfprinterrDC/usr/lib/python3/dist-packages/sympy/printing/tests/test_fortran.py_fcodesz%test_printmethod..nint._fcodeN)__name__ __module__ __qualname__rFrDrDrDrEnints rJz nint(x))rrr8)r?rJrDrDrEtest_printmethodsrKcstdtddd}tddd}ttddd d ksJtt|ddd d ks*Jtt|ddd d ks7Jttfd ddS)Nr?yTr.z)complex_freestandard source_formatz#merge(0d0, dsign(1d0, x), x == 0d0)zmerge(0, isign(1, y), y == 0)z/merge(cmplx(0d0, 0d0), z/abs(z), abs(z) == 0d0)cs ttSN)r8rrDr?rDrE$s z!test_fcode_sign..)rr8rr=NotImplementedError)rLrNrDrVrEtest_fcode_signs  rYcCs(td\}}tddd}t|ddksJt||ddks"Jtdt|d ||||d |d ksdS) NrZz sin(x)**cos(y)cttddSNBrSr8r rDr?rLrDrErWPz&test_fcode_functions..ctddSrnr8rDrrrDrErWQcrmNMrprqrDrrrDrErWRrscrtrwrurDrrrDrErWSrv)ZrPiirpz modulo(x, y))rr8rrr=rXr rprDrrrEtest_fcode_functionsMs rzc Cst}td\}}}}}}}tt|t||t||dks$Jtt|d|||tdddks;Jtt|t||t||dddksSJt|t|ddd ksaJ|j|t||d d d ksrJ|j|t|d d dksJ|j|dd dksJtddd\}} td} td}td}t d| } t d|} t| | | f|| || dddksJdS)Nzx,x_,x__,y,X,X_,Yz& exp(x_) + sin(x*y) + cos(X__*Y_)r^raz( 2*X_**(7.0d0/2.0d0)*Y*x + exp(x__)F)Z name_manglingz# exp(x_) + sin(x*y) + cos(X*Y)z x - cos(X)me assign_toz me = X*sin(x_) + x__Zmuz mu = X*sin(x_)Zadz ad = x__n,mTrMAr?rLirrQr}rTzndo i = 1, m y(i) = 0 end do do i = 1, m do I_ = 1, n y(i) = A(i, I_)*x(I_) + y(i) end do end do) r9rr8rrrrdoprintr:r;) Zobr?Zx_Zx__rLXZX_Yr[mrrrrDrDrE test_caseYs0"$&"    rcCstd}tdd}d}t|tdd|ksJt|tdd|ks)Jt|ttdd|ks9Jtttdd|ksGJttddtddksXJt|ttdd|kshJt|tttdd|kszJdS)Nr?r_z0.8340324452479558d0z x*%sd0z %sd0z x*%s)rrevalfr8rr)r?Zlog10_17Z loglog10_17rDrDrE"test_fcode_functions_with_integersxs " (rcCsd}t}ttdt|ksJttdt|ksJttdt|ks,Jttdt|ks9Jttdt|ksFJttdddtdksUJttd d t|t|fht d fkslJttd d t|t|fht d fksJttd d t|t|fht d fksJttd d t|t|fht dfksJttd d t|t|fht dfksJttdd dt|tdfht dfksJdS)Nrz. parameter (Catalan = %sd0) Catalanz4 parameter (EulerGamma = %sd0) EulerGammaz" parameter (E = %sd0) Ez6 parameter (GoldenRatio = %sd0) GoldenRatioz$ parameter (pi = %sd0) pi) precisionF)humanz Catalanz EulerGammaz Ez GoldenRatioz pi)rr) r9r8rrrrrrr@set)ZprecprDrDrEtest_fcode_NumberSymbolsH       rcCsttdksJtd}tdtdksJtddtdks"Jtddt|dks0Jtt|dks:Jtddt|d ksHJtdd d }td |d ksXJtt|dksbJtd|dkslJdS)Nz cmplx(0,1)r?z cmplx(0,4)r\z cmplx(3,4)z cmplx(3,4) + xz cmplx(0,1)*xz cmplx(3,4) - xT)Z imaginaryrz 5*xz x + 3)r8rrrVrDrDrEtest_fcode_complexs rcCsNtd\}}tt|dksJtt||dksJtt|dks%JdS)NrZz sin(x)z atan2(x, y)z conjg(x))rr8rrrrrrDrDrE test_implicits rcCsXtd}td}tt|}|dksJttt|dks Jt||dks*JdS)Nr?gz:C Not supported in Fortran: C gamma gamma(x)zHC Not supported in Fortran: C Integral Integral(sin(x), x)z2C Not supported in Fortran: C g g(x))rrr8rr r)r?rZgamma_frDrDrEtest_not_fortrans   rcCstd}tt|ddiddksJtd}tt|ddiddks$Jtd }t||d d idd ks6Jtd d d}tt|ddiddksJJdS)Nr?rZzsin)Zuser_functionsz zsin(x)rZmygammaz mygamma(x)rZgreatz great(x)r[TrMr Zfctz fct(n))rr8rrrr )r?rr[rDrDrEtest_user_functionss     rcCstd}tdt|d|}t||dksJtdt|dt|}t||dtdks3Jtd}tdtd d d }tdt||d |d|}t|||||d dks`JdS)Nr?rr^z 2*xz( parameter (pi = %sd0) 2*pi/xrrrr[TrMr]r|zI do i = 1, n A(i) = (A(i) + 1)*(A(i) + 2)*A(i) end do)rr<rr8rrr:r;)r?rrrrDrDrEtest_inline_functions    rcCs$td}tt|dddksJdS)Nr?sr|z s = sin(x))rr8rrVrDrDrEtest_assign_tosrcsZtd\}t|ddddksJfddtdD}tt|d ks+JdS) NrZr_varr|z var = x**10 + 10*x**9*y + 45*x**8*y**2 + 120*x**7*y**3 + 210*x**6* @ y**4 + 252*x**5*y**5 + 210*x**4*y**6 + 120*x**3*y**7 + 45*x**2*y @ **8 + 10*x*y**9 + y**10csg|]}|qSrDrD).0rrVrDrE sz&test_line_wrapping.. zR x**10 + x**9 + x**8 + x**7 + x**6 + x**5 + x**4 + x**3 + x**2 + x @ + 1)rr8expandranger )rLerDrVrEtest_line_wrappings   rcCstd\}}tt||k||dkdddksJtt||k||dkdddks,Jtt||k||dkddddd ksAJtt||k||dkddd ksTJdS) Nx yr]rQr`zx < y .and. y < x + 1zx < y .or. y < x + 1FZevaluatezx < y .neqv. y < x + 1zx < y .eqv. y < x + 1)rr8r1r2r5r4rrrDrDrEtest_fcode_precedences  rcCsttd\}}}tt|dddksJtt||dddks Jtt|t|dddks/Jttt||dddks>Jttt|t|dddksOJttt||d d ddd ks`Jtt||ddd ksmJtt|t|ddd ks|Jttt||dddksJttt|t|dddksJttt||d d dddksJttt|||dddksJttt|||dddksJttt|||dddksJttt|||dddksJttt|||dddksJttt|||dddksJtt|||dddksJtt||t|dddks.Jtt|t||dddks?Jttt|||dddksPJttt|||d d dddkscJtt|||dddksrJtt||t|dddksJtt|t||dddksJttt|||dddksJttt|||d d ddd ksJdS)!Nx y zrQr`z.not. xz x .and. yzx .and. .not. yzy .and. .not. xz.not. x .and. .not. yFrz.not. (x .and. y)zx .or. yzx .or. .not. yzy .or. .not. xz.not. x .or. .not. yz.not. (x .or. y)zx .and. (y .or. z)zy .and. (x .or. z)zz .and. (x .or. y)zx .or. y .and. zzy .or. x .and. zzz .or. x .and. yzx .and. y .and. zzx .and. y .and. .not. zzx .and. z .and. .not. yzy .and. z .and. .not. xz.not. (x .and. y .and. z)zx .or. y .or. zzx .or. y .or. .not. zzx .or. z .or. .not. yzy .or. z .or. .not. xz.not. (x .or. y .or. z))rr8r3r1r2r?rLrNrDrDrEtest_fcode_Logicalsh     "       rcCs(td\}}}tt||dddddksJtt|t|dddddks'Jttt||dddddks8Jttt|t|ddddd ksKJttt||ddddddd ks^Jtt||ddd kskJtt|t|ddd kszJttt||ddd ksJttt|t|dddksJttt||dddddksJttt|||dddksJttt|||dddksJttt|||dddksJttt|||dddksJttt|||dddksJttt|||dddks Jttt|||dddksJttt|||dddks.Jttt|||dddks?Jttt|||dddksPJttt|||dddksaJttt|||dddksrJttt||dd|dddksJttt||dd|dddksJttt||dd|dddksJttt|||dddddksJttt|||ddddd ksJttt|||ddddd!ksJttt|||ddddd"ksJttt|||ddddd#ks Jttt|||ddddd$ksJttt||dd|ddd%ks0Jttt||dd|ddd&ksCJttt||dd|ddd'ksVJttt|||ddddd(ksiJttt|||ddddd)ks|Jttt|||ddddd*ksJttt||dd|ddd+ksJttt||dd|ddd,ksJttt||dd|ddd-ksJtt|||ddddd.ksJtt||t|ddddd/ksJtt|t||ddddd0ksJttt|||ddddd1ksJdS)2NrFrrQr`z x .neqv. yzx .neqv. .not. yzy .neqv. .not. xz.not. x .neqv. .not. yz.not. (x .neqv. y)z x .eqv. yzx .eqv. .not. yzy .eqv. .not. xz.not. x .eqv. .not. yz.not. (x .eqv. y)zx .eqv. y .and. zzy .eqv. x .and. zzz .eqv. x .and. yzx .and. (y .eqv. z)zy .and. (x .eqv. z)zz .and. (x .eqv. y)zx .eqv. y .or. zzy .eqv. x .or. zzz .eqv. x .or. yzx .or. (y .eqv. z)zy .or. (x .eqv. z)zz .or. (x .eqv. y)zx .eqv. (y .neqv. z)zy .eqv. (x .neqv. z)zz .eqv. (x .neqv. y)zx .neqv. (y .eqv. z)zy .neqv. (x .eqv. z)zz .neqv. (x .eqv. y)zx .neqv. y .and. zzy .neqv. x .and. zzz .neqv. x .and. yzx .and. (y .neqv. z)zy .and. (x .neqv. z)zz .and. (x .neqv. y)zx .neqv. y .or. zzy .neqv. x .or. zzz .neqv. x .or. yzx .or. (y .neqv. z)zy .or. (x .neqv. z)zz .or. (x .neqv. y)zx .neqv. y .neqv. zzx .neqv. y .neqv. .not. zzx .neqv. z .neqv. .not. yzy .neqv. z .neqv. .not. x)rr8r5r3r4r1r2rrDrDrEtest_fcode_Xlogical+s*                            rcCstd\}}tt||ddddksJtt||ddddks"Jtt||dddd ks0Jtt||d ddd ks>Jtt||d ddd ksLJtt||ddddksZJdS)Nrz==rQr`zx == yz!=zx /= yz>=zx >= yz<=zx <= y>zx > y.rPrpz merge(x, x**2, x < 1)rr|zR if (x < 1) then var = x else var = x**2 end ifr_a/ if (x < 0) then weird_name = -cos(x)/x + 10*sin(x)/x**2 + 90*cos(x)/x**3 - 720* @ sin(x)/x**4 - 5040*cos(x)/x**5 + 30240*sin(x)/x**6 + 151200*cos(x @ )/x**7 - 604800*sin(x)/x**8 - 1814400*cos(x)/x**9 + 3628800*sin(x @ )/x**10 + 3628800*cos(x)/x**11 else weird_name = -sin(x)/x - 10*cos(x)/x**2 + 90*sin(x)/x**3 + 720* @ cos(x)/x**4 - 5040*sin(x)/x**5 - 30240*cos(x)/x**6 + 151200*sin(x @ )/x**7 + 604800*cos(x)/x**8 - 1814400*sin(x)/x**9 - 3628800*cos(x @ )/x**10 + 3628800*sin(x)/x**11 end ifrZ weird_namez1 merge(x, merge(x**2, sin(x), x > 1), x < 1)crrUrurDrrDrErWr) rr r=rXr8rrrr ValueError)r?codeexpectedabrrDrrEtest_fcode_Piecewises.         0 ,rcCsrt}gd}||}gd}|D] }t|dksJqt||D] \}}||ks,Jq"t|t|ks7JdS)N)zbC This is a long comment on a single line that must be wrapped properly to produce nice outputzf this = is + a + long + and + nasty + fortran + statement + that * must + be + wrapped + properlyzg this = is + a + long + and + nasty + fortran + statement + that * must + be + wrapped + properlyzh this = is + a + long + and + nasty + fortran + statement + that * must + be + wrapped + properlyzd this = is + a + long + and + nasty + fortran + statement + that*must + be + wrapped + properlyzf this = is + a + long + and + nasty + fortran + statement + that*must + be + wrapped + properlyzg this = is + a + long + and + nasty + fortran + statement + that*must + be + wrapped + properlyzh this = is + a + long + and + nasty + fortran + statement + that*must + be + wrapped + properlyze this = is + a + long + and + nasty + fortran + statement + that**must + be + wrapped + properlyzf this = is + a + long + and + nasty + fortran + statement + that**must + be + wrapped + properlyzg this = is + a + long + and + nasty + fortran + statement + that**must + be + wrapped + properlyzh this = is + a + long + and + nasty + fortran + statement + that**must + be + wrapped + properlyzi this = is + a + long + and + nasty + fortran + statement + that**must + be + wrapped + properlyzc this = is + a + long + and + nasty + fortran + statement(that)/must + be + wrapped + properlyzg this = is + a + long + and + nasty + fortran + statement(that)/must + be + wrapped + properly)zBC This is a long comment on a single line that must be wrappedz%C properly to produce nice outputzG this = is + a + long + and + nasty + fortran + statement + that *% @ must + be + wrapped + properlyzH this = is + a + long + and + nasty + fortran + statement + that *rG this = is + a + long + and + nasty + fortran + statement + thatz' @ * must + be + wrapped + properlyzF this = is + a + long + and + nasty + fortran + statement + that*rzH this = is + a + long + and + nasty + fortran + statement + that*rH this = is + a + long + and + nasty + fortran + statement + thatz& @ *must + be + wrapped + properly@ this = is + a + long + and + nasty + fortran + statement +z* @ that*must + be + wrapped + properlyzG this = is + a + long + and + nasty + fortran + statement + that**rzH this = is + a + long + and + nasty + fortran + statement + that**rr' @ **must + be + wrapped + properlyrrrz+ @ that**must + be + wrapped + properlyzE this = is + a + long + and + nasty + fortran + statement(that)/rzH this = is + a + long + and + nasty + fortran + statement(that)z& @ /must + be + wrapped + properlyH)r9 _wrap_fortranlenzip)rClinesZ wrapped_linesZexpected_lineslinewrrDrDrEtest_wrap_fortrans  rcCs,t}gd}gd}|||ksJdS)N)H this_variable_is_very_long_because_we_try_to_test_line_break=1.0d0zI this_variable_is_very_long_because_we_try_to_test_line_break =1.0d0zK this_variable_is_very_long_because_we_try_to_test_line_break = 1.0d0zL this_variable_is_very_long_because_we_try_to_test_line_break = 1.0d0zM this_variable_is_very_long_because_we_try_to_test_line_break = 1.0d0zK this_variable_is_very_long_because_we_try_to_test_line_break = 10.0d0) rD this_variable_is_very_long_because_we_try_to_test_line_break = @ 1.0d0zE this_variable_is_very_long_because_we_try_to_test_line_break =rzF this_variable_is_very_long_because_we_try_to_test_line_break =rzG this_variable_is_very_long_because_we_try_to_test_line_break =rrz @ 10.0d0r9rrCrrrDrDrEtest_wrap_fortran_keep_d0s rcCsttdddS)NcSsttdddS)NrZgarbage)method)r8rrDrDrDrErWrvztest_settings..)r= TypeErrorrDrDrDrE test_settingssrcCs0td\}}tt|t|dddksJdS)NrZrQr`zsin(y) + cos(x))rr8rrrrrDrDrEtest_free_form_code_lines $rcCs@td\}}tt|t|ddd}d}||ksJdS)NrZrarQr`zsin(y)**7 + 7*sin(y)**6*cos(x) + 21*sin(y)**5*cos(x)**2 + 35*sin(y)**4* & cos(x)**3 + 35*sin(y)**3*cos(x)**4 + 21*sin(y)**2*cos(x)**5 + 7* & sin(y)*cos(x)**6 + cos(x)**7)rr8rrr)r?rLresultrrDrDrE test_free_form_continuation_line"s  rcCs0tddi}dg}ddg}|||ksJdS)NrTrQz^! This is a long comment on a single line that must be wrapped properly to produce nice outputzG! This is a long comment on a single line that must be wrapped properlyz! to produce nice outputrrrDrDrEtest_free_form_comment_line-s rc Cstddd\}}td}td}td}td|}td|}d }t|||f||||d d }||d d iksS||d diksU||d diksW||d diksYJdSdSdSdS)Nr~TrMrr?rLrjz_do i = 1, m y(i) = 0 end do do i = 1, m do j = 1, n y(i) = %(rhs)s end do end dorQrrhszy(i) + A(i, j)*x(j)zy(i) + x(j)*A(i, j)zx(j)*A(i, j) + y(i)zA(i, j)*x(j) + y(i))rr:r;r8) r[rrr?rLrrrrrDrDrE test_loops6s   " rcCsftddtd\}}td}td}t||}d|jj|jd}t||||dd }||ks1JdS) Nzi mT)r.clsr?rLzMdo i_%(icount)i = 1, m_%(mcount)i y(i_%(icount)i) = x(i_%(icount)i) end do)ZicountZmcountrQr)rrr:r;ZlabelZ dummy_indexr8)rrr?rLrrrDrDrEtest_dummy_loopsPs rcCsd}td|fd}td|fd}td|dfd}td|d}t||||d||||d||}t|j|jdd }|d sIJdS) NrrL)shaper?Dyr]rF)r}Zcontractz+Dy(i) = (y(i + 1) - y(i))/(x(i + 1) - x(i)))r:r;rr8rZlhsendswith)Zlen_yrLr?rrrZcode0rDrDrE2test_fcode_Indexed_without_looking_for_contraction^s2rcCs:Gdddt}|}td}|t|ddksJdS)Nc@seZdZejZdS)z1test_derived_classes..MyFancyFCodePrinterN)rGrHrIr9Z_default_settingscopyrDrDrDrEMyFancyFCodePrinterjsrr?Zborkz bork = sin(x))r9rrr)rrCr?rDrDrEtest_derived_classesisrcCs.d}d}tddi}||}||ksJdS)Nausubroutine test(a) integer :: a, i, j do do do j = 1, 5 if (a>b) then if(b>0) then a = 3 donot_indent_me = 2 do_not_indent_me_either = 2 ifIam_indented_something_went_wrong = 2 if_I_am_indented_something_went_wrong = 2 end should not be unindented here end if endif end do end do enddo end subroutine subroutine test2(a) integer :: a do a = a + 1 end do end subroutine asubroutine test(a) integer :: a, i, j do do do j = 1, 5 if (a>b) then if(b>0) then a = 3 donot_indent_me = 2 do_not_indent_me_either = 2 ifIam_indented_something_went_wrong = 2 if_I_am_indented_something_went_wrong = 2 end should not be unindented here end if endif end do end do enddo end subroutine subroutine test2(a) integer :: a do a = a + 1 end do end subroutine rTrQ)r9Z indent_code)Z codelinesrrrrDrDrE test_indentrs  rc Cs,td\}}}t||td||dkf|dft|g}tddd}t||dks,Jtd|d |dkf|d dft|d |d }t|d d dksQJtddd}tddd}tt|d dt|d g|d |d |ddgd|d|d t|d ddgg}t||dksJdS)Nzx,y,zr^rTrr\r]z A(1, 1) = x*y if (y > 0) then A(2, 1) = x + 2 else A(2, 1) = y end if A(3, 1) = sin(z))r^r)r]rrrrPrpz? merge(2*A(3, 1), A(3, 1), x > 0) + sin(A(2, 1)) + A(1, 1)qrM)r\r)rrrz M(1, 1) = sin(q(2, 1)) M(2, 1) = q(2, 1) + q(3, 1) M(3, 1) = 2*q(5, 1)/q(2, 1) M(1, 2) = 0 M(2, 2) = q(4, 1) M(3, 2) = sqrt(q(1, 1)) + 4 M(1, 3) = cos(q(3, 1)) M(2, 3) = 5 M(3, 3) = 0)rr6r rr7r8rr ) r?rLrNZmatrrrrrrDrDrEtest_Matrix_printings&, 6   $ rcCsDtd\}}t|tdddt|||g}t|}|dks JdS)Nrrr_r^z3 do x = 0, 10, 2 y = x*y end do)rr!rr"r8)r?rLfZsolrDrDrEtest_fcode_Fors  rc Csdd}tddd}t|}t|}||dtd\}}t|tdthd }t|}||d t|td }t|} || d |td |td|td ttid|tdttiddS)Nc[s"t|fddd||ksJdS)NrPrQrRru)rrefkwargsrDrDrEchecks"z%test_fcode_Declaration..checkrTrMzinteger*4 :: ir*)valueZattrszreal*4, parameter :: x = 42typez logical :: yzreal*4zreal*8)Z type_aliases) rr%Zdeducedr$r&r(r*r'r)) rrZvar1Zdcl1r?rLZvar2Zdcl2bZvar3Zdcl3rDrDrEtest_fcode_Declarations          rcCsxtddd}tddd}tddd}t|ddksJtd|ddks(J|d|||}t|d ks:JdS) Nrr]r\BCrz A(1, 1)z 3*A(1, 1)z (A - B)(1, 1))r7r8Zsubs)rrrFrDrDrEtest_MatrixElement_printings   rcCs(td}tt|dddddksJdS)Nr?+r]rQr`z x = x + 1)rr8r#rVrDrDrEtest_aug_assigns rcCs8td}ttt|dkt|ddgdddksJdS)Nr?r]-rQr`z)do while (abs(x) > 1) x = x - 1 end do)rr8r+absr#rVrDrDrE test_While s" rcsRtd}tddd}t|td}t|td}ttd||gttfdddS) Nr?r[TrMrpowercrrUrurDfp1rDrErWrz.test_FunctionPrototype_print..)rr%r)r.r,r=rX)r?r[vxvnrDrrEtest_FunctionPrototype_prints    rcsjtd}tddd}t|td}t|td}t|||t|g}ttd||g|ttfdddS) Nr?r[TrMrrcrrUrurDZfd1rDrErW&rz/test_FunctionDefinition_print..) rr%r)r.r"r/r-r=rX)r?r[rrbodyrDrrEtest_FunctionDefinition_prints   rcCs4t ddl}WddS1swYdS)Nr)r>Zsympy.printing.fcode)sympyrDrDrEtest_fcode_submodule(s "rN)qrrrrrrrrr r r r r rrrrrrrrrrrrrrrrrrr Z sympy.codegenr!r"r#Zsympy.codegen.astr$r%r&r'r(r)r*r+r,r-r.r/Zsympy.core.relationalr0Zsympy.logic.boolalgr1r2r3r4r5Zsympy.matricesr6r7Zsympy.printing.fortranr8r9Z sympy.tensorr:r;Zsympy.utilities.lambdifyr<Zsympy.testing.pytestr=r>rKrYrbrfrjrlrzrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrDrDrDrEsh8          1c +<    ?#