B < [OO @srddlmZdgZddlZddlmZddlZddlZddl m Z ddl m Z ddl mZmZmZmZmZddlmZmZmZdd lmZGd ddeZd d Zd dZedddedddedddedddedddedddedddedddgZddgZddZddZ d d!Z!d"d#Z"Gd$d%d%eZ#d&d'Z$d(d)Z%d*d+Z&d,d-d.d/d0d1d2d3d4d5g Z'd6d7Z(dS)8)print_functionLinearConstraintN)Mapping) PatsyError)Origin)atleast_2d_column_defaultrepr_pretty_delegaterepr_pretty_impl no_picklingassert_no_pickling)TokenOperator infix_parse)_parsing_error_testc@s6eZdZdZd ddZeZddZeZ e ddZ dS) rapA linear constraint in matrix form. This object represents a linear constraint of the form `Ax = b`. Usually you won't be constructing these by hand, but instead get them as the return value from :meth:`DesignInfo.linear_constraint`. .. attribute:: coefs A 2-dimensional ndarray with float dtype, representing `A`. .. attribute:: constants A 2-dimensional single-column ndarray with float dtype, representing `b`. .. attribute:: variable_names A list of strings giving the names of the variables being constrained. (Used only for consistency checking.) NcCst||_ttj|td|_|dkr>tj|jjdtd}tj|td}t ||_ |j j dksr|j jddkrzt d|jj dks|jjdt |krt d|jjddkrt d|jjd|j jdkrt ddS) N)dtyperz1constants is not (convertible to) a column matrixzwrong shape for coefsz/must have at least one row in constraint matrixz*shape mismatch between coefs and constants)listvariable_namesnpZ atleast_2dasarrayfloatcoefszerosshaper constantsndim ValueErrorlen)selfrrrr /lib/python3.7/site-packages/patsy/constraint.py__init__/s   zLinearConstraint.__init__cCs |rtt|||j|j|jgS)N)AssertionErrorr rrr)rpcycler r r! _repr_pretty_@szLinearConstraint._repr_pretty_cCsj|s td|dj}x|D]}|j|krtdqWtdd|D}tdd|D}||||S)aCreate a new LinearConstraint by ANDing together several existing LinearConstraints. :arg constraints: An iterable of LinearConstraint objects. Their :attr:`variable_names` attributes must all match. :returns: A new LinearConstraint object. zno constraints specifiedrzvariable names don't matchcSsg|] }|jqSr )r).0cr r r! Vsz,LinearConstraint.combine..cSsg|] }|jqSr )r)r'r(r r r!r)Ws)rrrZ row_stack)cls constraintsr constraintrrr r r!combineGs     zLinearConstraint.combine)N) __name__ __module__ __qualname____doc__r"r__repr__r&r __getstate__ classmethodr-r r r r!rs  cCsyddlm}Wn tk r0ddlm}YnXtddgddg}|jddgksVt||jddgg||jdggtddgddgddggdd g}||jddgddgg||jdgd gg|jj t t kst|jj t t ksttd gdgg}||jdggdd l m }|ttd gddgg|ttd gdggg|ttd gddggdd g|ttd d gddggdd g|ttd gdgggg|ttd d gg|ttd d gt dt|dS)Nr) assert_equalZfooZbarrr a) assert_raisesb)rr) numpy.testingr5 ImportErrornumpy.testing.utilsrrr#rrrrr nose.toolsr:rrr )r5Zlcr:r r r!test_LinearConstraintZs2    rAcCsttddgddgtddgddgdgg}|jddgks@tyddlm}Wn tk rpddlm}YnX||jddgddgg||j dgdggddl m }|t tjg|t tjtdgdgtdgdggdS)Nr9r<rr)r5)r:) rr-rr#r=r5r>r?rrr@r:r)Zcombr5r:r r r!test_LinearConstraint_combines rB,ri=+rd-*/NUMBERVARIABLEcsfdd}|S)Ncs.dkr|}n}t|tf|j|S)N__OP__)r rmatchspan)scannerZ token_stringZ actual_type)stringtyper r! make_tokens z _token_maker..make_tokenr )rRrQrSr )rQrRr! _token_makersrTc Csd}d}dddtD}d}d}t|tdd }dd d|D}|ttj|f|ttj|f|td |f|td |f|td |f|dfg}t |} | |\} } | rt|t| } t dt || | d| S)Nz\(z\)|cSsg|]}t|jqSr )reescape token_type)r'opr r r!r)sz(_tokenize_constraint..z&[-+]?[0-9]*\.?[0-9]+([eE][-+]?[0-9]+)?z\s+T)keyreversecSsg|]}t|qSr )rVrW)r'nr r r!r)srMrLrKz unrecognized token in constraintr) join_opssortedrrTr LPARENRPARENrVScannerscanrr) rQrZ lparen_reZ rparen_reZop_reZnum_reZ whitespace_reZ variable_relexiconrPtokensZleftoveroffsetr r r!_tokenize_constraints(     rgc Csd}t|dddg}ddtjddd fd d d tjd ddfddg }xft||D]X\}}t|ts`t|j|dksrt|jt ||d|dkst|j |dksJtqJWddl m }|t tddgtddgxNddgddgfD]:}td|}t|dkstdd|DdddgkstqWtddg}t|dks.z 2 * a[1,1],za[1,1]cSsg|] }|jqSr )rR)r'rwr r r!r)srKrHrLrCcSsg|] }|jqSr )rv)r'rwr r r!r)sri)rgr r`razip isinstancer#rRoriginrrvr@r:rr)codereZ expectedsgotexpectedr:namesr r r!test__tokenize_constraints6        rcCstt||ttS)N)rrgr^_atomic)rQrr r r!parse_constraints rc@sveZdZddZddZddZddZd d Zd d Zd dZ ddZ ddZ ddZ ddZ ddZdddZdS)_EvalConstraintc CsD||_t||_|j|j|j|j|j|j|j |j |j |j d |_ dS)N) )rLr)rKr)rEr)rGr)rEr)rGr)rHr)rJr)rDr)rCr)_variable_namesr_N_eval_variable _eval_number_eval_unary_plus_eval_unary_minus_eval_binary_plus_eval_binary_minus_eval_binary_multiply_eval_binary_div_eval_binary_eq_eval_binary_comma _dispatch)rrr r r!r"s z_EvalConstraint.__init__cCst|d|jdkS)Nr)rallr)rrr r r! is_constant sz_EvalConstraint.is_constantcCs2|jj}tj|jdftd}d||j|<|S)Nr)r)tokenrvrrrrrindex)rtreevarrr r r!r sz_EvalConstraint._eval_variablecCs*tj|jdftd}t|jj|d<|S)Nr)r)rrrrrrv)rrrr r r!rsz_EvalConstraint._eval_numbercCs||jdS)Nr)evalargs)rrr r r!rsz _EvalConstraint._eval_unary_pluscCsd||jdS)Nrr)rr)rrr r r!rsz!_EvalConstraint._eval_unary_minuscCs ||jd||jdS)Nrr)rr)rrr r r!rsz!_EvalConstraint._eval_binary_pluscCs ||jd||jdS)Nrr)rr)rrr r r!r sz"_EvalConstraint._eval_binary_minuscCsF||jd}||jd}||s:td|jd||dS)Nrrz1Can't divide by a variable in a linear constraintr)rrrr)rrleftrightr r r!r#s   z _EvalConstraint._eval_binary_divcCsZ||jd}||jd}||r6|d|S||rL||dStd|dS)Nrrrz=Can't multiply one variable by another in a linear constraint)rrrr)rrrrr r r!r+s    z%_EvalConstraint._eval_binary_multiplyc Cst|j}g}xDt|D]8\}}|jdkr||j|dd|jd|||<qW||d}||d}|d|j|d|j}t|dkrt d||d |d} t |j || } || t |S)NrDT)r,rrz!no variables appear in constraintr) rr enumeraterRappendrrrrrrrr-) rrrr+iargrrrZconstantr,r r r!r6s    z_EvalConstraint._eval_binary_eqcCs6|j|jddd}|j|jddd}t||gS)NrT)r,r)rrrr-)rrrrr r r!rJsz"_EvalConstraint._eval_binary_commaFcCs|jt|jf}||jkst|j||}|rt|tr>|S|j|jdksRtt |d|jdkrtt d|t|j |d|j|d Snt|trt d||SdS)Nrrz#term is constant, with no variablesrzunexpected constraint object) rRrrrr#ryrsizerrrrr)rrr,rZvalr r r!rOs     z_EvalConstraint.evalN)F)r.r/r0r"rrrrrrrrrrrrr r r r!rs rcCst|tr*|j|kr&td|j|f|St|trtjt|t|ftd}tt|}t }xt t |D]x\}\}}||kr| |}n t|t jr|}ntd|f||krtd||f||d|||f<|||<qrWt|||St|tr |g}t|tr|rt|dtrg} xN|D]F} t| tsVtd| ft| |} t|} | | j| dd q8Wt| St|trt|d krtd |\} }t|| |Stj|td}t||S) zaThis is the internal interface implementing DesignInfo.linear_constraint, see there for docs.z?LinearConstraint has wrong variable_names (got %r, expected %r))rz#unrecognized variable name/index %rzduplicated constraint on %rrrzexpected a string, not %rT)r,rz#constraint tuple must have length 2)ryrrrrrrrrsetrsixZ iteritemsrZ integer_typesaddstrrrrrrr-tupler)Zconstraint_likerrrZusedrnamevalueidxr+r{rZ evaluatorZcoefr r r!linear_constraintesZ                   rcCsyddlm}Wn tk r0ddlm}YnXt||}td|t|||}td|||j|j||j|j||j |j ||jj t t ||j j t t dS)Nr)r5r|r}) r=r5r>r?rprintrrrrrrr)inputZvarnamesrrr5r|r}r r r! _check_lincons    rcCsddlm}ddlm}t}|tddgddgddgddggdgg|tttddgddgddg|ddiddgddggdgg||d d gddgddgddggdgdgg||d d gddgddgddggdgdgg|ddiddgddggdgg||d d gddgddgddggdgdgg||d d gddgddgddggdgdgg|ttd diddg|ttdddddg|t ddgddgddggdgg|t ddgddggddgddgddggdgdgg|dddgddggdgg|dddgddggdgg|dddgddggdgg|dddgddggdgg|dddgddgddggdgdgg|dddgddgddggdgdgg|ddgddgddgddggdgdgg|ttdddigddg|dddgddggdgg|dddgddggdgg|dddgddgddggdgdgg|dddgddgddgddggdgdgdgg|d ddgddggdgg|d!ddgddggdgg|d"ddgddggdgg|d#ddgddgddggdgdgg|d$ddgddgddggdgdgg|d%ddd&gdddgdddgdddggdgdgdgg|d'dd(gddggdgg|d)d*gd+gfddgd)d*ggd+gg|d)d*gd*d,ggd+gd-ggfddgd)d*gd*d,ggd+gd-gg|ttddgdgdgfddg|ttddgfddg|d)d*gddgd)d*ggdgg|d)d*gd*d,ggddgd)d*gd*d,ggdgdgg|t d)d*gddgd)d*ggdgg|t d)d*gd*d,ggddgd)d*gd*d,ggdgdgg|ttdddgdS).Nr)r:) OrderedDictr9r<rr6r)r9r)r<r6)rr)rr6rh)r9rr;rjza = 2za - 2z a + 1 = 3z a + b = 3z a = 2, b = 3z b = 3, a = 2zb = 3z#2 * (a + b/3) + b + 2*3/4 = 1 + 2*3g?g@z+2 * -aza - b, a + b = 2rza = 1, a = 2, a = 3za * 2z-a = 1z(2 + a - a) * bz a = 1 = bz a = (1 = b)za = 1, a = b = cr(z a + 1 = 2za + 1r7r8(#) r@r:Z patsy.compatrrrrrrZarray)r:rrwr r r!test_linear_constraintsl  *   " " . "&8&&* & &&& ,&" 2&rz a + ooza = 1, <1 = 1>, b = 1za = 1, za = 1, <1>, b = 2za = 1, <2 * b = b + b>, cza + + cza + 2 / + cza = 1, 2 * <(a = b)>, cza = 1, a + <(a = b)>, cza = 1, <(a, b)> + 2, ccCsdd}t|tdS)NcSst|dddgS)Nr9r<r()r)Zbad_coder r r!doitsztest_eval_errors..doit)r_parse_eval_error_tests)rr r r!test_eval_errorssr))Z __future__r__all__rV collectionsrrZnumpyrZpatsyrZ patsy.originrZ patsy.utilrrr r r Zpatsy.infix_parserr r rZpatsy.parse_formularobjectrrArBr^rrTrgrrrrrrrrr r r r!sR     A%        $u>N