ó ¡¼™\c@sQddlmZmZmZmZmZmZmZmZm 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/m0Z0ddl1m2Z2ddl3m4Z4ddl5m6Z6m7Z7m8Z8m9Z9m:Z:ddl;m<Z<m=Z=ddl>m?Z?m@Z@dd lAmBZBmCZCdd lDmEZEdd lFmGZGd „ZHd „ZId„ZJd„ZKd„ZLd„ZMd„ZNd„ZOd„ZPd„ZQd„ZRd„ZSd„ZTd„ZUd„ZVd„ZWd„ZXd„ZYd„ZZd„Z[d „Z\d!„Z]d"„Z^d#„Z_d$„Z`d%„Zad&„Zbd'„Zcd(„Zdd)„Zed*„Zfd+„Zgd,„Zhd-„Zid.„Zjd/„Zkd0„Zld1„Zmd2„Znd3„Zod4„Zpd5S(6iÿÿÿÿ(tsintcostatan2tlogtexptgammat conjugatetsqrtt factorialtIntegralt PiecewisetAddtdifftsymbolstStFloattDummytEqtRangetCatalant EulerGammatEt GoldenRatiotItpitFunctiontRationaltIntegertLambdatsigntMod(tFort Assignmentt aug_assign( t DeclarationtVariabletfloat32tfloat64t value_consttrealtbool_tWhiletFunctionPrototypetFunctionDefinitiontintegertReturn(trange(t Relational(tAndtOrtNott EquivalenttXor(tMatrixt MatrixSymbol(tfcodet FCodePrinter(t IndexedBasetIdx(timplemented_function(traisescCsDtdƒ}dtfd„ƒY}t||ƒƒdks@t‚dS(NtxtnintcBseZd„ZRS(cSsd|j|jdƒS(Nsnint(%s)i(t_printtargs(tselftprinter((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyt_fcodes(t__name__t __module__RC(((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyR>ss nint(x)(R RR7tAssertionError(R=R>((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyttest_printmethods csÈtdƒ‰tddtƒ}tddtƒ}ttˆƒdddd ƒd ksZt‚tt|ƒdddd ƒd ks„t‚tt|ƒdddd ƒd ks®t‚tt‡fd †ƒdS(NR=tyR,tztcomplextstandardi_t source_formattfrees#merge(0d0, dsign(1d0, x), x == 0d0)smerge(0, isign(1, y), y == 0)s/merge(cmplx(0d0, 0d0), z/abs(z), abs(z) == 0d0)csttˆƒƒS(N(R7R((R=(s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyt%t(R tTrueR7RRFR<tNotImplementedError(RHRI((R=s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyttest_fcode_signs ***cCsœtdƒ\}}tddtƒ}t|dƒdks@t‚t||dƒdks`t‚tdt|ƒd||||d |ƒd ksžt‚tt|ƒƒd ks¼t‚tt|ƒƒd ksÚt‚t|d ƒd ksöt‚tt|ƒƒd kst‚ttdƒƒdks2t‚t|dƒdksNt‚t|ddddƒdksst‚t|tddƒƒdks˜t‚dS(Nsx,ytnR,is x**3s x**(y**3)ig @is, (3.5d0*sin(x))**(-x + y**x)/(x**2 + y)s sqrt(x)s sqrt(dble(n))gà?i s sqrt(10.0d0)gð¿s 1d0/xgÀRHRLRMsy = x**(-2.0d0)is x**(3.0d0/7.0d0)(R RPR7RFRRR(R=RHRS((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyttest_fcode_Pow(s /%cCsÞtdƒ}ttddƒƒdks-t‚ttddƒƒdksNt‚ttddƒƒd ksot‚ttd dƒƒdkst‚t|tddƒƒd ksµt‚ttddƒ|ƒd ksÚt‚dS( NR=iis 3.0d0/7.0d0ii s 2iùÿÿÿs -3.0d0/7.0d0iýÿÿÿs x + 3.0d0/7.0d0s (3.0d0/7.0d0)*x(R R7RRF(R=((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyttest_fcode_Rational:s !!!!%cCs@ttdƒƒdkst‚ttdƒƒdks<t‚dS(NiCs 67iÿÿÿÿs -1(R7RRF(((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyttest_fcode_IntegerDscCs@ttdƒƒdkst‚ttdƒƒdks<t‚dS(NgE@s 42.0000000000000d0g@Œµx¯Äs -1.00000000000000d+20(R7RRF(((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyttest_fcode_FloatIscstdƒ\‰‰ttˆƒtˆƒƒdks:t‚tt‡‡fd†ƒtt‡‡fd†ƒtt‡‡fd†ƒtt‡‡fd†ƒxcddd d gD]O}ttˆˆƒd |ƒd ksÞt‚tˆˆd |ƒd ks±t‚q±WdS( Nsx,ys sin(x)**cos(y)csttˆˆƒddƒS(NRKiB(R7R((R=RH(s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyRNQROcstˆˆddƒS(NRKiB(R7((R=RH(s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyRNRROcsttˆˆƒddƒS(NRKiM(R7R((R=RH(s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyRNSROcstˆˆddƒS(NRKiM(R7((R=RH(s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyRNTROiZi_iÓiØRKs modulo(x, y)(R R7RRRFR<RQR(RK((R=RHs>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyttest_fcode_functionsNs('c Cstƒ}tdƒ\}}}}}}}tt|ƒt||ƒt||ƒƒdksdt‚tt|ƒd|||tddƒƒdksŸt‚tt|ƒt||ƒt||ƒdtƒdksßt‚t|t|ƒdtƒdkst‚|j |t|ƒ|d d ƒd ks6t‚|j |t|ƒd d ƒd ksat‚|j |d dƒdks‚t‚tddt ƒ\}} t dƒ} t dƒ}t dƒ}t d| ƒ} t d|ƒ} t| | | f|| d || ddƒdkst‚dS(Nsx,x_,x__,y,X,X_,Ys& exp(x_) + sin(x*y) + cos(X__*Y_)iis( 2*X_**(7.0d0/2.0d0)*Y*x + exp(x__)t name_manglings# exp(x_) + sin(x*y) + cos(X*Y)s x - cos(X)t assign_totmes me = X*sin(x_) + x__tmus mu = X*sin(x_)tads ad = x__sn,mR,tAR=RHtiRRLRMsndo 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( R8R R7RRRRFRtFalsetdoprintRPR9R:( tobR=tx_tx__RHtXtX_tYRStmR^R_R((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyt test_caseZs( !+,1(/+!   +cCsZtdƒ}tdƒjdƒ}d}t|tdƒƒd|ksMt‚t|tdƒƒd|ksst‚t|ttdƒƒƒd|ksŸt‚tttdƒƒƒd|ksÇt‚ttdƒƒdtdƒjdƒksøt‚t|ttdƒƒƒd|ks$t‚t|tttdƒƒƒƒd|ksVt‚dS(NR=i is0.8340324452479558d0s x*%sd0s %sd0s x*%s(R RtevalfR7RFRR(R=tlog10_17t loglog10_17((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyt"test_fcode_functions_with_integersys &&,(1,cCsÑd}tƒ}ttƒdtj|ƒks4t‚ttƒdtj|ƒksYt‚ttƒdtj|ƒks~t‚ttƒdtj|ƒks£t‚ttƒdtj|ƒksÈt‚ttddƒdtjdƒksót‚ttd t ƒt t|j tj|ƒƒfgƒt gƒd fksAt‚ttd t ƒt t|j tj|ƒƒfgƒt gƒd fkst‚ttd t ƒt t|j tj|ƒƒfgƒt gƒd fksÝt‚ttd t ƒt t|j tj|ƒƒfgƒt gƒd fks+t‚ttd t ƒt t|j tj|ƒƒfgƒt gƒdfksyt‚ttddd t ƒt t|j tjdƒƒfgƒt gƒdfksÍt‚dS(Nis. parameter (Catalan = %sd0) Catalans4 parameter (EulerGamma = %sd0) EulerGammas" parameter (E = %sd0) Es6 parameter (GoldenRatio = %sd0) GoldenRatios$ parameter (pi = %sd0) pit precisionithumans Catalans EulerGammas Es GoldenRatios pi( R8R7RRjRFRRRRR`tsetR?(tprectp((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyttest_fcode_NumberSymbol†s* %%%%%(<3?3?cCs.ttƒdkst‚tdƒ}tdtƒdks@t‚tddtƒdks`t‚tddt|ƒdks„t‚tt|ƒdks t‚tddt|ƒd ksÄt‚tdd tƒ}td |ƒd ksòt‚tt|ƒdkst‚td|ƒd ks*t‚dS(Ns cmplx(0,1)R=is cmplx(0,4)is cmplx(3,4)s cmplx(3,4) + xs cmplx(0,1)*xs cmplx(3,4) - xt imaginaryis 5*xs x + 3(R7RRFR RP(R=((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyttest_fcode_complexžs  $$cCsstdƒ\}}tt|ƒƒdks0t‚tt||ƒƒdksQt‚tt|ƒƒdksot‚dS(Nsx,ys sin(x)s atan2(x, y)s conjg(x)(R R7RRFRR(R=RH((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyt test_implicit¬s!cCs‚tdƒ}tdƒ}tt|ƒƒ}|dks<t‚ttt|ƒƒƒdks`t‚t||ƒƒdks~t‚dS(NR=tgs:C Not supported in Fortran: C gamma gamma(x)sHC Not supported in Fortran: C Integral Integral(sin(x), x)s2C Not supported in Fortran: C g g(x)(R RR7RRFR R(R=Rwtgamma_f((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyttest_not_fortran³s   $cCsætdƒ}tt|ƒdidd6ƒdks7t‚tdƒ}tt|ƒdidd6ƒdksnt‚td ƒ}t||ƒdid d 6ƒd ks¥t‚td d tƒ}tt|ƒdidd6ƒdksât‚dS(NR=tuser_functionstzsinRs zsin(x)tmygammaRs mygamma(x)Rwtgreats great(x)RSR,tfctRs fct(n)(R R7RRFRRRPR(R=RwRS((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyttest_user_functions¼s + ( +cCstdƒ}tdt|d|ƒƒ}t||ƒƒdksFt‚tdt|dt|ƒƒ}t||ƒƒdtjdƒks‘t‚tdƒ}tdtd d t ƒƒ}tdt||d |d|ƒƒ}t|||ƒd ||ƒd ks t‚dS(NR=Rwis 2*xs( parameter (pi = %sd0) 2*pi/xiR^R_RSR,iRZsI do i = 1, n A(i) = (A(i) + 1)*(A(i) + 2)*A(i) end do( R R;RR7RFRRjR9R:RP(R=RwR^R_((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyttest_inline_functionÉs   (cCs4tdƒ}tt|ƒddƒdks0t‚dS(NR=RZtss s = sin(x)(R R7RRF(R=((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyttest_assign_toÜs cCsƒtdƒ\}}t||djƒddƒdks>t‚gtdƒD]}||^qK}tt|Œƒdkst‚dS(Nsx,yi RZtvars¯ 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**10i sR x**10 + x**9 + x**8 + x**7 + x**6 + x**5 + x**4 + x**3 + x**2 + x @ + 1(R R7texpandRFR.R (R=RHR_te((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyttest_line_wrappingás #cCsøtdƒ\}}tt||k||dkƒddƒdksIt‚tt||k||dkƒddƒdks€t‚tt||k||dkdtƒddƒdks½t‚tt||k||dkƒddƒd ksôt‚dS( Nsx yiRLRMsx < y .and. y < x + 1sx < y .or. y < x + 1tevaluatesx < y .neqv. y < x + 1sx < y .eqv. y < x + 1(R R7R0RFR1R4R`R3(R=RH((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyttest_fcode_precedenceïs((((cCs tdƒ\}}}tt|ƒddƒdks9t‚tt||ƒddƒdks`t‚tt|t|ƒƒddƒdkst‚ttt|ƒ|ƒddƒdksºt‚ttt|ƒt|ƒƒddƒdksít‚ttt||ƒd tƒddƒd ks t‚tt||ƒddƒd ksGt‚tt|t|ƒƒddƒd kstt‚ttt|ƒ|ƒddƒd ks¡t‚ttt|ƒt|ƒƒddƒdksÔt‚ttt||ƒd tƒddƒdkst‚ttt||ƒ|ƒddƒdks7t‚ttt||ƒ|ƒddƒdksgt‚ttt||ƒ|ƒddƒdks—t‚ttt||ƒ|ƒddƒdksÇt‚ttt||ƒ|ƒddƒdks÷t‚ttt||ƒ|ƒddƒdks't‚tt|||ƒddƒdksQt‚tt||t|ƒƒddƒdkst‚tt|t|ƒ|ƒddƒdks±t‚ttt|ƒ||ƒddƒdksát‚ttt|||ƒd tƒddƒdkst‚tt|||ƒddƒdksAt‚tt||t|ƒƒddƒdksqt‚tt|t|ƒ|ƒddƒdks¡t‚ttt|ƒ||ƒddƒdksÑt‚ttt|||ƒd tƒddƒdkst‚dS( Nsx y zRLRMs.not. xs x .and. ysx .and. .not. ysy .and. .not. xs.not. x .and. .not. yR‡s.not. (x .and. y)sx .or. ysx .or. .not. ysy .or. .not. xs.not. x .or. .not. ys.not. (x .or. y)sx .and. (y .or. z)sy .and. (x .or. z)sz .and. (x .or. y)sx .or. y .and. zsy .or. x .and. zsz .or. x .and. ysx .and. y .and. zsx .and. y .and. .not. zsx .and. z .and. .not. ysy .and. z .and. .not. xs.not. (x .and. y .and. z)sx .or. y .or. zsx .or. y .or. .not. zsx .or. z .or. .not. ysy .or. z .or. .not. xs.not. (x .or. y .or. z)(R R7R2RFR0R`R1(R=RHRI((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyttest_fcode_LogicalûsP$'--$$'--$$000000*!!!'*!!!'cCsãtdƒ\}}}tt||dtƒddƒdksBt‚tt|t|ƒdtƒddƒdksut‚ttt|ƒ|dtƒddƒdks¨t‚ttt|ƒt|ƒdtƒddƒdksát‚ttt||dtƒdtƒddƒd kst‚tt||ƒddƒd ksAt‚tt|t|ƒƒddƒd ksnt‚ttt|ƒ|ƒddƒd ks›t‚ttt|ƒt|ƒƒddƒd ksÎt‚ttt||ƒdtƒddƒdkst‚ttt||ƒ|ƒddƒdks1t‚ttt||ƒ|ƒddƒdksat‚ttt||ƒ|ƒddƒdks‘t‚ttt||ƒ|ƒddƒdksÁt‚ttt||ƒ|ƒddƒdksñt‚ttt||ƒ|ƒddƒdks!t‚ttt||ƒ|ƒddƒdksQt‚ttt||ƒ|ƒddƒdkst‚ttt||ƒ|ƒddƒdks±t‚ttt||ƒ|ƒddƒdksát‚ttt||ƒ|ƒddƒdkst‚ttt||ƒ|ƒddƒdksAt‚ttt||dtƒ|ƒddƒdkswt‚ttt||dtƒ|ƒddƒdks­t‚ttt||dtƒ|ƒddƒdksãt‚ttt||ƒ|dtƒddƒdkst‚ttt||ƒ|dtƒddƒdksOt‚ttt||ƒ|dtƒddƒd ks…t‚ttt||ƒ|dtƒddƒd!ks»t‚ttt||ƒ|dtƒddƒd"ksñt‚ttt||ƒ|dtƒddƒd#ks't‚ttt||dtƒ|ƒddƒd$ks]t‚ttt||dtƒ|ƒddƒd%ks“t‚ttt||dtƒ|ƒddƒd&ksÉt‚ttt||ƒ|dtƒddƒd'ksÿt‚ttt||ƒ|dtƒddƒd(ks5t‚ttt||ƒ|dtƒddƒd)kskt‚ttt||dtƒ|ƒddƒd*ks¡t‚ttt||dtƒ|ƒddƒd+ks×t‚ttt||dtƒ|ƒddƒd,ks t‚tt|||dtƒddƒd-ks=t‚tt||t|ƒdtƒddƒd.ksst‚tt|t|ƒ|dtƒddƒd/ks©t‚ttt|ƒ||dtƒddƒd0ksßt‚dS(1Nsx y zR‡RLRMs x .neqv. ysx .neqv. .not. ysy .neqv. .not. xs.not. x .neqv. .not. ys.not. (x .neqv. y)s x .eqv. ysx .eqv. .not. ysy .eqv. .not. xs.not. x .eqv. .not. ys.not. (x .eqv. y)sx .eqv. y .and. zsy .eqv. x .and. zsz .eqv. x .and. ysx .and. (y .eqv. z)sy .and. (x .eqv. z)sz .and. (x .eqv. y)sx .eqv. y .or. zsy .eqv. x .or. zsz .eqv. x .or. ysx .or. (y .eqv. z)sy .or. (x .eqv. z)sz .or. (x .eqv. y)sx .eqv. (y .neqv. z)sy .eqv. (x .neqv. z)sz .eqv. (x .neqv. y)sx .neqv. (y .eqv. z)sy .neqv. (x .eqv. z)sz .neqv. (x .eqv. y)sx .neqv. y .and. zsy .neqv. x .and. zsz .neqv. x .and. ysx .and. (y .neqv. z)sy .and. (x .neqv. z)sz .and. (x .neqv. y)sx .neqv. y .or. zsy .neqv. x .or. zsz .neqv. x .or. ysx .or. (y .neqv. z)sy .or. (x .neqv. z)sz .or. (x .neqv. y)sx .neqv. y .neqv. zsx .neqv. y .neqv. .not. zsx .neqv. z .neqv. .not. ysy .neqv. z .neqv. .not. x( R R7R4R`RFR2R3R0R1(R=RHRI((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyttest_fcode_Xlogical,s°$$$$'$!!!!!!!!!!!!!!!!!!''''''''''''!'''cCstdƒ\}}tt||dƒddƒdks<t‚tt||dƒddƒdksft‚tt||dƒddƒd kst‚tt||d ƒddƒd ksºt‚tt||d ƒddƒd ksät‚tt||dƒddƒdkst‚dS(Nsx ys==RLRMsx == ys!=sx /= ys>=sx >= ys<=sx <= yt>sx > ytlib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyttest_fcode_Relationals*****cstdƒ}t||dkf|dtfƒ‰tt‡fd†ƒtˆddƒ}d}||ksqt‚tt||dkf|dtfƒdd ƒd ks®t‚t|ƒ|}t|ƒ|}x2t d ƒD]$}t ||ƒ}t ||ƒ}qÛWd }tt||d kf|tfƒddƒ}||ksHt‚tt||dkf|d|dkft|ƒtfƒddƒ}d}||ks¦t‚t||dkf|d|dkft|ƒ|d kfƒ‰tt ‡fd†ƒdS(NR=iics tˆƒS(N(R7((texpr(s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyRNRORKi_s merge(x, x**2, x < 1)RZRƒsR if (x < 1) then var = x else var = x**2 end ifi s/ 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 ifit weird_names1 merge(x, merge(x**2, sin(x), x > 1), x < 1)cs tˆƒS(N(R7((RŽ(s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyRNÁRO( R R RPR<RQR7RFRRR.R t ValueError(R=tcodetexpectedtatbR_((RŽs>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyttest_fcode_Piecewise™s* %. -F@cCs(tƒ}ddddddddd d d d d ddg}|j|ƒ}dddddddddddddddddddddddddd d!dd"d#g}x&|D]}t|ƒd$ks²t‚q²Wx/t||ƒD]\}}||ksät‚qäWt|ƒt|ƒks$t‚dS(%NsbC This is a long comment on a single line that must be wrapped properly to produce nice outputsf this = is + a + long + and + nasty + fortran + statement + that * must + be + wrapped + properlysg this = is + a + long + and + nasty + fortran + statement + that * must + be + wrapped + properlysh this = is + a + long + and + nasty + fortran + statement + that * must + be + wrapped + properlysd this = is + a + long + and + nasty + fortran + statement + that*must + be + wrapped + properlysf this = is + a + long + and + nasty + fortran + statement + that*must + be + wrapped + properlysg this = is + a + long + and + nasty + fortran + statement + that*must + be + wrapped + properlysh this = is + a + long + and + nasty + fortran + statement + that*must + be + wrapped + properlyse this = is + a + long + and + nasty + fortran + statement + that**must + be + wrapped + properlysf this = is + a + long + and + nasty + fortran + statement + that**must + be + wrapped + properlysg this = is + a + long + and + nasty + fortran + statement + that**must + be + wrapped + properlysh this = is + a + long + and + nasty + fortran + statement + that**must + be + wrapped + properlysi this = is + a + long + and + nasty + fortran + statement + that**must + be + wrapped + properlysc this = is + a + long + and + nasty + fortran + statement(that)/must + be + wrapped + properlysg this = is + a + long + and + nasty + fortran + statement(that)/must + be + wrapped + properlysBC This is a long comment on a single line that must be wrappeds%C properly to produce nice outputsG this = is + a + long + and + nasty + fortran + statement + that *s% @ must + be + wrapped + properlysH this = is + a + long + and + nasty + fortran + statement + that *sG this = is + a + long + and + nasty + fortran + statement + thats' @ * must + be + wrapped + properlysF this = is + a + long + and + nasty + fortran + statement + that*sH this = is + a + long + and + nasty + fortran + statement + that*sH this = is + a + long + and + nasty + fortran + statement + thats& @ *must + be + wrapped + properlys@ this = is + a + long + and + nasty + fortran + statement +s* @ that*must + be + wrapped + properlysG this = is + a + long + and + nasty + fortran + statement + that**sH this = is + a + long + and + nasty + fortran + statement + that**s' @ **must + be + wrapped + properlys+ @ that**must + be + wrapped + properlysE this = is + a + long + and + nasty + fortran + statement(that)/sH this = is + a + long + and + nasty + fortran + statement(that)s& @ /must + be + wrapped + properlyiH(R8t _wrap_fortrantlenRFtzip(RBtlinest wrapped_linestexpected_linestlinetwR…((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyttest_wrap_fortranÄsh    c Csgtƒ}ddddddg}dddd dd dd ddd g }|j|ƒ|ksct‚dS( NsH this_variable_is_very_long_because_we_try_to_test_line_break=1.0d0sI this_variable_is_very_long_because_we_try_to_test_line_break =1.0d0sK this_variable_is_very_long_because_we_try_to_test_line_break = 1.0d0sL this_variable_is_very_long_because_we_try_to_test_line_break = 1.0d0sM this_variable_is_very_long_because_we_try_to_test_line_break = 1.0d0sK this_variable_is_very_long_because_we_try_to_test_line_break = 10.0d0sD this_variable_is_very_long_because_we_try_to_test_line_break =s @ 1.0d0sE this_variable_is_very_long_because_we_try_to_test_line_break =sF this_variable_is_very_long_because_we_try_to_test_line_break =sG this_variable_is_very_long_because_we_try_to_test_line_break =s @ 10.0d0(R8R–RF(RBR™R’((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyttest_wrap_fortran_keep_d0s&   cCsttd„ƒdS(NcSsttdƒddƒS(Nitmethodtgarbage(R7R(((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyRNRO(R<t TypeError(((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyt test_settingsscCsDtdƒ\}}tt|ƒt|ƒddƒdks@t‚dS(Nsx,yRLRMssin(y) + cos(x)(R R7RRRF(R=RH((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyttest_free_form_code_linescCsZtdƒ\}}tt|ƒt|ƒdjƒddƒ}d}||ksVt‚dS(Nsx,yiRLRMsµsin(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(R R7RRR„RF(R=RHtresultR’((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyt test_free_form_continuation_line#s,cCsGtidd6ƒ}dg}ddg}|j|ƒ|ksCt‚dS(NRMRLs^! This is a long comment on a single line that must be wrapped properly to produce nice outputsG! This is a long comment on a single line that must be wrapped properlys! to produce nice output(R8R–RF(RBR™R’((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyttest_free_form_comment_line.s   c Csôtddtƒ\}}tdƒ}tdƒ}tdƒ}td|ƒ}td|ƒ}d}t|||f||d ||d d ƒ}||id d 6ksð||idd 6ksð||idd 6ksð||idd 6ksðt‚dS(Nsn,mR,R^R=RHR_tjs_do i = 1, m y(i) = 0 end do do i = 1, m do j = 1, n y(i) = %(rhs)s end do end doRZRLRMsy(i) + A(i, j)*x(j)trhssy(i) + x(j)*A(i, j)sx(j)*A(i, j) + y(i)sA(i, j)*x(j) + y(i)(R RPR9R:R7RF( RSRhR^R=RHR_R¨R’R‘((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyt test_loops7s    .cCsœtddtdtƒ\}}tdƒ}tdƒ}t||ƒ}di|jjd6|jd6}t||d ||d d ƒ}||ks˜t‚dS( Nsi mR,tclsR=RHsMdo i_%(icount)i = 1, m_%(mcount)i y(i_%(icount)i) = x(i_%(icount)i) end doticounttmcountRZRLRM( R RPRR9R:tlabelt dummy_indexR7RF(R_RhR=RHR’R‘((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyttest_dummy_loopsQs   cCsÊd}tdd|fƒ}tdd|fƒ}tdd|dfƒ}td|dƒ}t||||d||||d||ƒ}t|jd|jd tƒ}|jd ƒsÆt‚dS( NiRHtshapeR=tDyiR_RZtcontracts+Dy(i) = (y(i + 1) - y(i))/(x(i + 1) - x(i))( R9R:RR7R©tlhsR`tendswithRF(tlen_yRHR=R²R_R…tcode0((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyt2test_fcode_Indexed_without_looking_for_contraction_s7cCsSdtfd„ƒY}|ƒ}tdƒ}|jt|ƒdƒdksOt‚dS(NtMyFancyFCodePrintercBseZejjƒZRS((RDRER8t_default_settingstcopy(((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyR¹ksR=tborks bork = sin(x)(R8R RaRRF(R¹RBR=((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyttest_derived_classesjs  cCsDd}d}tidd6ƒ}|j|ƒ}||ks@t‚dS(Nsusubroutine 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 ssubroutine 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 RMRL(R8t indent_codeRF(t codelinesR’RrR¥((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyt test_indentss c Cs–tdƒ\}}}t||td||dkf|tfƒt|ƒgƒ}tdddƒ}t||ƒdks€t‚td|d|dkf|dtfƒt|dƒ|d}t|dd ƒd ksát‚td d dƒ}td ddƒ}tt|dƒdt|dƒg|d|d|dd gd|d|dt |dƒddggƒ}t||ƒdks’t‚dS(Nsx,y,ziiR^iis† 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)RKi_s? merge(2*A(3, 1), A(3, 1), x > 0) + sin(A(2, 1)) + A(1, 1)tqitMisí 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(ii(ii(ii(ii(ii(ii(ii(ii(ii(ii(ii(ii( R R5R RPRR6R7RFRR( R=RHRItmatR^RŽRÁRÂRh((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyttest_Matrix_printing²s> C#3 cCsbtdƒ\}}t|tdddƒt|||ƒgƒ}t|ƒ}|dks^t‚dS(Nsx yii is3 do x = 0, 10, 2 y = x*y end do(R RRR R7RF(R=RHtftsol((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyttest_fcode_ForÕs. c Cs d„}tddtƒ}tj|ƒ}t|ƒ}||dƒtdƒ\}}t|tdddthƒ}t|ƒ}||d ƒt|d tƒ}t|ƒ} || d ƒ|td ƒ|td ƒ|t d ditt 6ƒ|t d ditt 6ƒdS(Nc[s+t|dddd||ks't‚dS(NRKi_RLRM(R7RF(RŽtreftkwargs((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pytcheckàsR_R,sinteger*4 :: isx ytvaluei*tattrssreal*4, parameter :: x = 42ttypes logical :: ysreal*4sreal*8t type_aliases( R RPR#tdeducedR"R$R&R(R%R'( RÊR_tvar1tdcl1R=RHtvar2tdcl2btvar3tdcl3((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyttest_fcode_Declarationßs          cCs¨tdddƒ}tdddƒ}tdddƒ}t|d ƒdksRt‚td|d ƒdksrt‚|d j|||ƒ}t|ƒd ks¤t‚dS( NR^iitBtCis A(1, 1)s 3*A(1, 1)s (A - B)(1, 1)(ii(ii(ii(R6R7RFtsubs(R^R×RØtF((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyttest_MatrixElement_printingøs cCs:tdƒ}tt|ddƒddƒdks6t‚dS(NR=t+iRLRMs x = x + 1(R R7R!RF(R=((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyttest_aug_assigns cCsRtdƒ}ttt|ƒdkt|ddƒgƒddƒdksNt‚dS(NR=it-RLRMs)do while (abs(x) > 1) x = x - 1 end do(R R7R)tabsR!RF(R=((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyt test_While s 3csttdƒ}tddtƒ}t|dtƒ}t|dtƒ}ttd||gƒ‰tt‡fd†ƒdS(NR=RSR,RÍtpowercs tˆƒS(N(R7((tfp1(s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyRNRO(R RPR#R'R,R*R<RQ(R=RStvxtvn((Râs>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyttest_FunctionPrototype_prints  cs–tdƒ}tddtƒ}t|dtƒ}t|dtƒ}t|||ƒt|ƒg}ttd||g|ƒ‰tt ‡fd†ƒdS(NR=RSR,RÍRács tˆƒS(N(R7((tfd1(s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyRN'RO( R RPR#R'R,R R-R+R<RQ(R=RSRãRätbody((Ræs>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pyttest_FunctionDefinition_prints N(qtsympyRRRRRRRRRR R R R R RRRRRRRRRRRRRRRRRt sympy.codegenRR R!tsympy.codegen.astR"R#R$R%R&R'R(R)R*R+R,R-tsympy.core.compatibilityR.tsympy.core.relationalR/tsympy.logic.boolalgR0R1R2R3R4tsympy.matricesR5R6tsympy.printing.fcodeR7R8t sympy.tensorR9R:tsympy.utilities.lambdifyR;tsympy.utilities.pytestR<RGRRRTRURVRWRXRiRmRsRuRvRyRR€R‚R†RˆR‰RŠRR•RžRŸR£R¤R¦R§RªR°R¸R½RÀRÄRÇRÖRÛRÝRàRåRè(((s>lib/python2.7/site-packages/sympy/printing/tests/test_fcode.pytsfÄR(           1 c + <      ? #