The R Project SVN R

Rev

Rev 60595 | Rev 64401 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 60595 Rev 62908
Line 657... Line 657...
657 
    Rcomplex *work, tmp;
 657 
    Rcomplex *work, tmp;
658 
 
 658 
 
659 
    k = LENGTH(tau);
 659 
    k = LENGTH(tau);
660 
    if (!isMatrix(Bin)) error(_("'b' must be a complex matrix"));
 660 
    if (!isMatrix(Bin)) error(_("'b' must be a complex matrix"));
661 
 
 661 
 
662 
    if (!isReal(Bin)) B = PROTECT(coerceVector(Bin, CPLXSXP));
 662 
    if (!isComplex(Bin)) B = PROTECT(coerceVector(Bin, CPLXSXP));
663 
    else B = PROTECT(duplicate(Bin));
 663 
    else B = PROTECT(duplicate(Bin));
664 
 
 664 
 
665 
    n = INTEGER(coerceVector(getAttrib(qr, R_DimSymbol), INTSXP))[0];
 665 
    n = INTEGER(coerceVector(getAttrib(qr, R_DimSymbol), INTSXP))[0];
666 
    Bdims = INTEGER(coerceVector(getAttrib(Bin, R_DimSymbol), INTSXP));
 666 
    Bdims = INTEGER(coerceVector(getAttrib(Bin, R_DimSymbol), INTSXP));
667 
    if(Bdims[0] != n)
 667 
    if(Bdims[0] != n)