Changeset 1125
- Timestamp:
- Jun 29, 2004, 10:28:47 AM (22 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/psLib/src/dataManip/psMatrixVectorArithmetic.c
r1124 r1125 29 29 * @author Ross Harman, MHPCC 30 30 * 31 * @version $Revision: 1.1 2$ $Name: not supported by cvs2svn $32 * @date $Date: 2004-06-29 20: 16:16$31 * @version $Revision: 1.13 $ $Name: not supported by cvs2svn $ 32 * @date $Date: 2004-06-29 20:28:47 $ 33 33 * 34 34 * Copyright 2004 Maui High Performance Computing Center, University of Hawaii … … 580 580 } else if(!strncmp(OP, "exp", 3)) { \ 581 581 if(PS_IS_PSELEMTYPE_COMPLEX(IN->type)) { \ 582 UNARY_TYPE(DIM,OUT,IN,cexp(*i1)); 583 } else { \ 584 UNARY_TYPE(DIM,OUT,IN,exp(*i1)); \582 UNARY_TYPE(DIM,OUT,IN,cexp(*i1)); \ 583 } else { \ 584 UNARY_TYPE(DIM,OUT,IN,exp(*i1)); \ 585 585 } \ 586 586 } else if(!strncmp(OP, "ln", 2)) { \ 587 587 if(PS_IS_PSELEMTYPE_COMPLEX(IN->type)) { \ 588 UNARY_TYPE(DIM,OUT,IN,clog(*i1)); 589 } else { \ 590 UNARY_TYPE(DIM,OUT,IN,log(*i1)); \588 UNARY_TYPE(DIM,OUT,IN,clog(*i1)); \ 589 } else { \ 590 UNARY_TYPE(DIM,OUT,IN,log(*i1)); \ 591 591 } \ 592 592 } else if(!strncmp(OP, "ten", 3)) { \ 593 593 if(PS_IS_PSELEMTYPE_COMPLEX(IN->type)) { \ 594 UNARY_TYPE(DIM,OUT,IN,cpow(*i1,10.0)); 595 } else { \ 596 UNARY_TYPE(DIM,OUT,IN,exp10(*i1)); \594 UNARY_TYPE(DIM,OUT,IN,cpow(*i1,10.0)); \ 595 } else { \ 596 UNARY_TYPE(DIM,OUT,IN,exp10(*i1)); \ 597 597 } \ 598 598 } else if(!strncmp(OP, "log", 3)) { \ 599 599 if(PS_IS_PSELEMTYPE_COMPLEX(IN->type)) { \ 600 UNARY_TYPE(DIM,OUT,IN,clog10(*i1)); 601 } else { \ 602 UNARY_TYPE(DIM,OUT,IN,log10(*i1)); \600 UNARY_TYPE(DIM,OUT,IN,clog10(*i1)); \ 601 } else { \ 602 UNARY_TYPE(DIM,OUT,IN,log10(*i1)); \ 603 603 } \ 604 604 } else if(!strncmp(OP, "sin", 3)) { \ 605 605 if(PS_IS_PSELEMTYPE_COMPLEX(IN->type)) { \ 606 UNARY_TYPE(DIM,OUT,IN,csin(*i1)); 607 } else { \ 608 UNARY_TYPE(DIM,OUT,IN,sin(*i1)); \606 UNARY_TYPE(DIM,OUT,IN,csin(*i1)); \ 607 } else { \ 608 UNARY_TYPE(DIM,OUT,IN,sin(*i1)); \ 609 609 } \ 610 610 } else if(!strncmp(OP, "dsin", 4)) { \ 611 611 if(PS_IS_PSELEMTYPE_COMPLEX(IN->type)) { \ 612 UNARY_TYPE(DIM,OUT,IN,csin(*i1*D2R)); 613 } else { \ 614 UNARY_TYPE(DIM,OUT,IN,sin(*i1*D2R)); \612 UNARY_TYPE(DIM,OUT,IN,csin(*i1*D2R)); \ 613 } else { \ 614 UNARY_TYPE(DIM,OUT,IN,sin(*i1*D2R)); \ 615 615 } \ 616 616 } else if(!strncmp(OP, "cos", 3)) { \ 617 617 if(PS_IS_PSELEMTYPE_COMPLEX(IN->type)) { \ 618 UNARY_TYPE(DIM,OUT,IN,ccos(*i1)); 619 } else { \ 620 UNARY_TYPE(DIM,OUT,IN,cos(*i1)); \618 UNARY_TYPE(DIM,OUT,IN,ccos(*i1)); \ 619 } else { \ 620 UNARY_TYPE(DIM,OUT,IN,cos(*i1)); \ 621 621 } \ 622 622 } else if(!strncmp(OP, "dcos", 4)) { \ 623 623 if(PS_IS_PSELEMTYPE_COMPLEX(IN->type)) { \ 624 UNARY_TYPE(DIM,OUT,IN,ccos(*i1*D2R)); 625 } else { \ 626 UNARY_TYPE(DIM,OUT,IN,cos(*i1*D2R)); \624 UNARY_TYPE(DIM,OUT,IN,ccos(*i1*D2R)); \ 625 } else { \ 626 UNARY_TYPE(DIM,OUT,IN,cos(*i1*D2R)); \ 627 627 } \ 628 628 } else if(!strncmp(OP, "tan", 3)) { \ 629 629 if(PS_IS_PSELEMTYPE_COMPLEX(IN->type)) { \ 630 UNARY_TYPE(DIM,OUT,IN,ctan(*i1)); 631 } else { \ 632 UNARY_TYPE(DIM,OUT,IN,tan(*i1)); \630 UNARY_TYPE(DIM,OUT,IN,ctan(*i1)); \ 631 } else { \ 632 UNARY_TYPE(DIM,OUT,IN,tan(*i1)); \ 633 633 } \ 634 634 } else if(!strncmp(OP, "dtan", 4)) { \ 635 635 if(PS_IS_PSELEMTYPE_COMPLEX(IN->type)) { \ 636 UNARY_TYPE(DIM,OUT,IN,ctan(*i1*D2R)); 637 } else { \ 638 UNARY_TYPE(DIM,OUT,IN,tan(*i1*D2R)); \636 UNARY_TYPE(DIM,OUT,IN,ctan(*i1*D2R)); \ 637 } else { \ 638 UNARY_TYPE(DIM,OUT,IN,tan(*i1*D2R)); \ 639 639 } \ 640 640 } else if(!strncmp(OP, "asin", 4)) { \ 641 641 if(PS_IS_PSELEMTYPE_COMPLEX(IN->type)) { \ 642 UNARY_TYPE(DIM,OUT,IN,casin(*i1)); 643 } else { \ 644 UNARY_TYPE(DIM,OUT,IN,asin(*i1)); \642 UNARY_TYPE(DIM,OUT,IN,casin(*i1)); \ 643 } else { \ 644 UNARY_TYPE(DIM,OUT,IN,asin(*i1)); \ 645 645 } \ 646 646 } else if(!strncmp(OP, "dasin", 5)) { \ 647 647 if(PS_IS_PSELEMTYPE_COMPLEX(IN->type)) { \ 648 UNARY_TYPE(DIM,OUT,IN,R2D*casin(*i1)); 649 } else { \ 650 UNARY_TYPE(DIM,OUT,IN,R2D*asin(*i1)); \648 UNARY_TYPE(DIM,OUT,IN,R2D*casin(*i1)); \ 649 } else { \ 650 UNARY_TYPE(DIM,OUT,IN,R2D*asin(*i1)); \ 651 651 } \ 652 652 } else if(!strncmp(OP, "acos", 4)) { \ 653 653 if(PS_IS_PSELEMTYPE_COMPLEX(IN->type)) { \ 654 UNARY_TYPE(DIM,OUT,IN,cacos(*i1)); 655 } else { \ 656 UNARY_TYPE(DIM,OUT,IN,acos(*i1)); \654 UNARY_TYPE(DIM,OUT,IN,cacos(*i1)); \ 655 } else { \ 656 UNARY_TYPE(DIM,OUT,IN,acos(*i1)); \ 657 657 } \ 658 658 } else if(!strncmp(OP, "dacos", 5)) { \ 659 659 if(PS_IS_PSELEMTYPE_COMPLEX(IN->type)) { \ 660 UNARY_TYPE(DIM,OUT,IN,R2D*cacos(*i1)); 661 } else { \ 662 UNARY_TYPE(DIM,OUT,IN,R2D*acos(*i1)); \660 UNARY_TYPE(DIM,OUT,IN,R2D*cacos(*i1)); \ 661 } else { \ 662 UNARY_TYPE(DIM,OUT,IN,R2D*acos(*i1)); \ 663 663 } \ 664 664 } else if(!strncmp(OP, "atan", 4)) { \ 665 665 if(PS_IS_PSELEMTYPE_COMPLEX(IN->type)) { \ 666 UNARY_TYPE(DIM,OUT,IN,catan(*i1)); 667 } else { \ 668 UNARY_TYPE(DIM,OUT,IN,atan(*i1)); \666 UNARY_TYPE(DIM,OUT,IN,catan(*i1)); \ 667 } else { \ 668 UNARY_TYPE(DIM,OUT,IN,atan(*i1)); \ 669 669 } \ 670 670 } else if(!strncmp(OP, "datan", 5)) { \ 671 671 if(PS_IS_PSELEMTYPE_COMPLEX(IN->type)) { \ 672 UNARY_TYPE(DIM,OUT,IN,R2D*catan(*i1)); 673 } else { \ 674 UNARY_TYPE(DIM,OUT,IN,R2D*atan(*i1)); \672 UNARY_TYPE(DIM,OUT,IN,R2D*catan(*i1)); \ 673 } else { \ 674 UNARY_TYPE(DIM,OUT,IN,R2D*atan(*i1)); \ 675 675 } \ 676 676 } else { \
Note:
See TracChangeset
for help on using the changeset viewer.
