PDL::Ops PDL::Ops Perl module contains fundamental mathematical operators.
Download

PDL::Ops Ranking & Summary

PDL::Ops Tags

mathematical Perl module mathematical operators fundamental operators fundamental

PDL::Ops Description

PDL::Ops Perl module contains fundamental mathematical operators. PDL::Ops Perl module contains fundamental mathematical operators.This module provides the functions used by PDL to overload the basic mathematical operators (+ - / * etc.) and functions (sin sqrt etc.)It also includes the function log10, which should be a perl function so that we can overload it!Matrix multiplication (the operator x) is handled by the module PDL::Primitive.$doc $c = $name $a, $b, 0; # explicit call with trailing 0 $c = $a $op $b; # overloaded call $a->inplace->$name($b,0); # modify $a inplaceIt can be made to work inplace with the $a->inplace syntax. This function is used to overload the binary $optxt operator. Note that when calling this function explicitly you need to supply a third argument that should generally be zero (see first example). This restriction is expected to go away in future releases. EOD } # sub: biop()#simple binary functions sub bifunc { my ($name,$func,$swap,$doc,%extra) = @_; my $funcov = ref $func eq 'ARRAY' ? $func-> : $func; my $isop=0; if ($funcov =~ s/^op//) { $isop = 1; } my $funcovp = protect_chars $funcov; $func = $func-> if ref $func eq 'ARRAY'; if ($swap) { $extra{HdrCode} .= 1, NoBadifNaN => 1, Pars => 'a(); b(); c();', OtherPars => 'int swap', Inplace => , # quick and dirty solution to get ->inplace do its job Code => "$c() = $func($a(),$b());", BadCode => 'if ( $ISBAD(a()) || $ISBAD(b()) ) $SETBAD(c()); else' . "n $c() = $func($a(),$b());n", CopyBadStatusCode => 'if ( $BADFLAGCACHE() ) { if ( a == c && $ISPDLSTATEGOOD(a) ) { PDL->propogate_badflag( c, 1 ); /* have inplace op AND badflag has changed */ } $SETPDLSTATEBAD(c); }', %extra, Doc => $name($b,0); # explicit function call $ovcall $a->inplace->$name($b,0); # modify $a inplaceIt can be made to work inplace with the $a->inplace syntax. This function is used to overload the binary $funcovp function. Note that when calling this function explicitly you need to supply a third argument that should generally be zero (see first example). This restriction is expected to go away in future releases. EOD } # sub: bifunc()# simple unary functions and operators sub ufunc { my ($name,$func,$doc,%extra) = @_; my $funcov = ref $func eq 'ARRAY' ? $func-> : $func; my $funcovp = protect_chars $funcov; $func = $func-> if ref $func eq 'ARRAY'; # handle exceptions my $badcode = '$ISBAD(a())'; if ( exists $extra{Exception} ) {# $badcode .= " || $extra{Exception}";# print "Warning: ignored exception for $namen"; delete $exists{Exception}; } # do not have to worry about propogation of the badflag when # inplace since only input piddle is a, hence its badflag # won't change # UNLESS an exception occurs... pp_def($name, Pars => 'a(); b()', HandleBad => 1, NoBadifNaN => 1, Inplace => 1, Code => "$b() = $func($a());", BadCode => 'if ( ' . $badcode . ' ) $SETBAD(b()); else' . "n $b() = $func($a());n", %extra, Doc => inplace->$name; # modify $a inplaceIt can be made to work inplace with the $a->inplace syntax. This function is used to overload the unary $funcovp operator/function. EOD } # sub: ufunc()####################################################################### we trap some illegal operations here -- see the Exception option # note, for the ufunc()'s, the checks do not work too well # for unsigned integer types (ie < 0) # # XXX needs thinking about # - have to integrate into Code section as well (so # 12/pdl(2,4,0,3) is trapped and flagged bad) # --> complicated # - perhaps could use type %{ %} ? # # ==> currently have commented out the exception code, since # want to see if can use NaN/Inf for bad values # (would solve many problems for F,D types) # # there is an issue over how we handle comparison operators # - see Primitive/primitive.pd/zcover() for more discussion ### arithmetic ops # no swap biop('plus','+',0,'add two piddles'); biop('mult','*',0,'multiply two piddles');# all those need swapping biop('minus','-',1,'subtract two piddles'); biop('divide','/',1,'divide two piddles', Exception => '$b() == 0' );## note: divide should perhaps trap division by zero as well## comparison ops # need swapping biop('gt','>',1,'the binary > (greater than) operation'); biop('lt','',1,'leftshift a$ by $b',GenericTypes => $T); biop('or2','|',0,'binary or of two piddles',GenericTypes => $T); biop('and2','&',0,'binary and of two piddles',GenericTypes => $T); biop('xor','^',0,'binary exclusive or of two piddles',GenericTypes => $T);# really an ufunc ufunc('bitnot','~','unary bit negation',GenericTypes => $T);# some standard binary functions bifunc('power',,1,'raise piddle $a to the power b',GenericTypes => ); bifunc('atan2','atan2',1,'elementwise atan2 of two piddles',GenericTypes => ); bifunc('modulo',,1,'elementwise modulo operation'); bifunc('spaceship',,1,'elementwise ~ operation');# some standard unary functions ufunc('sqrt','sqrt','elementwise square root', Exception => '$a() < 0' ); ufunc('abs',,'elementwise absolute value',GenericTypes => ); ufunc('sin','sin','the sin function'); ufunc('cos','cos','the cos function'); ufunc('not','!','the elementwise not operation'); ufunc('exp','exp','the exponential function',GenericTypes => ); ufunc('log','log','the natural logarithm',GenericTypes => , Exception => '$a() is_inplace ) { $x->set_inplace(0); $y = $x; } elsif( ref($x) eq "PDL"){ #PDL Objects, use nullcreate: $y = PDL->nullcreate($x); }else{ #PDL-Derived Object, use copy: (Consistent with # Auto-creation docs in Objects.pod) $y = $x->copy; } &PDL::_log10_int( $x, $y ); return $y; }; ' );# note: the extra code that adding 'HandleBad => 1' creates is # unneeded here. Things could be made clever enough to work this out, # but it's very low priority. # It does add doc information though, and lets people know it's been # looked at for bad value support # DJB adds: not completely sure about this now that I have added code # to avoid a valgrind-reported error (see the CacheBadFlagInit rule # in PP.pm) # # Can't this be handled in Core.pm when '.=' is overloaded ? # pp_def( 'assgn', # HandleBad => 1, Pars => 'a(); b();', Code => '$b() = $a();', # BadCode => # 'if ( $ISBAD(a()) ) { $SETBAD(b()); } else { $b() = $a(); }', Doc => 'Plain numerical assignment. This is used to implement the ".=" operator', ); # pp_def assgn#pp_export_nothing();pp_done();Requirements:· Perl Requirements: · Perl

PDL::Ops Related Software