… | |
… | |
648 | { |
648 | { |
649 | int i; |
649 | int i; |
650 | AV *av, *result; |
650 | AV *av, *result; |
651 | |
651 | |
652 | if (!SvROK (ar) || SvTYPE (SvRV (ar)) != SVt_PVAV) |
652 | if (!SvROK (ar) || SvTYPE (SvRV (ar)) != SVt_PVAV) |
653 | croak ("Not an array ref as first argument to make_histogram"); |
653 | croak ("Not an array ref as first argument to extract_features"); |
654 | |
654 | |
655 | av = (AV *) SvRV (ar); |
655 | av = (AV *) SvRV (ar); |
656 | result = newAV (); |
656 | result = newAV (); |
657 | |
657 | |
658 | for (i = 0; i <= av_len (av); ++i) |
658 | for (i = 0; i <= av_len (av); ++i) |
… | |
… | |
721 | RETVAL = newRV_noinc ((SV *)result); |
721 | RETVAL = newRV_noinc ((SV *)result); |
722 | } |
722 | } |
723 | OUTPUT: |
723 | OUTPUT: |
724 | RETVAL |
724 | RETVAL |
725 | |
725 | |
726 | SV * |
|
|
727 | make_histograms (SV *ar) |
|
|
728 | CODE: |
|
|
729 | { |
|
|
730 | int i; |
|
|
731 | AV *av, *result; |
|
|
732 | |
|
|
733 | if (!SvROK (ar) || SvTYPE (SvRV (ar)) != SVt_PVAV) |
|
|
734 | croak ("Not an array ref as first argument to make_histogram"); |
|
|
735 | |
|
|
736 | av = (AV *) SvRV (ar); |
|
|
737 | result = newAV (); |
|
|
738 | |
|
|
739 | for (i = 0; i <= av_len (av); ++i) |
|
|
740 | { |
|
|
741 | const int HISTSIZE = 64; |
|
|
742 | |
|
|
743 | int j; |
|
|
744 | SV *sv = *av_fetch (av, i, 1); |
|
|
745 | STRLEN len; |
|
|
746 | char *buf = SvPVbyte (sv, len); |
|
|
747 | |
|
|
748 | int tmphist[HISTSIZE]; |
|
|
749 | float *hist; |
|
|
750 | |
|
|
751 | SV *histsv = newSV (HISTSIZE * sizeof (float) + 1); |
|
|
752 | SvPOK_on (histsv); |
|
|
753 | SvCUR_set (histsv, HISTSIZE * sizeof (float)); |
|
|
754 | hist = (float *)SvPVX (histsv); |
|
|
755 | |
|
|
756 | Zero (tmphist, sizeof (tmphist), char); |
|
|
757 | |
|
|
758 | for (j = len; j--; ) |
|
|
759 | { |
|
|
760 | unsigned int idx |
|
|
761 | = ((*buf & 0xc0) >> 2) |
|
|
762 | | ((*buf & 0x18) >> 1) |
|
|
763 | | (*buf & 0x03); |
|
|
764 | |
|
|
765 | ++tmphist[idx]; |
|
|
766 | ++buf; |
|
|
767 | } |
|
|
768 | |
|
|
769 | for (j = 0; j < HISTSIZE; ++j) |
|
|
770 | hist[j] = (float)tmphist[j] / (len + 1e-30); |
|
|
771 | |
|
|
772 | av_push (result, histsv); |
|
|
773 | } |
|
|
774 | |
|
|
775 | RETVAL = newRV_noinc ((SV *)result); |
|
|
776 | } |
|
|
777 | OUTPUT: |
|
|
778 | RETVAL |
|
|
779 | |
|
|
780 | |
|
|