1 |
#include "EXTERN.h" |
2 |
#include "perl.h" |
3 |
#include "XSUB.h" |
4 |
|
5 |
#include <string.h> |
6 |
#include <setjmp.h> |
7 |
#include <math.h> |
8 |
|
9 |
#include <magic.h> |
10 |
|
11 |
#include <jpeglib.h> |
12 |
#include <jerror.h> |
13 |
|
14 |
#include <glib.h> |
15 |
#include <gtk/gtk.h> |
16 |
#include <gdk/gdkx.h> |
17 |
#include <gdk-pixbuf/gdk-pixbuf.h> |
18 |
|
19 |
#include <gperl.h> |
20 |
#include <gtk2perl.h> |
21 |
|
22 |
#include <assert.h> |
23 |
|
24 |
#if WEBP |
25 |
#include <webp/demux.h> |
26 |
#include <webp/decode.h> |
27 |
#endif |
28 |
|
29 |
#if JXL |
30 |
#include <jxl/decode.h> |
31 |
#include "jxl/thread_parallel_runner.h" |
32 |
#endif |
33 |
|
34 |
#include "perlmulticore.h" |
35 |
|
36 |
#define IW 80 /* MUST match Schnauzer.pm! */ |
37 |
#define IH 60 /* MUST match Schnauzer.pm! */ |
38 |
|
39 |
#define LINELENGTH 240 |
40 |
|
41 |
#define ELLIPSIS "\xe2\x80\xa6" |
42 |
|
43 |
typedef char *octet_string; |
44 |
|
45 |
static magic_t magic_cookie[2]; /* !mime, mime */ |
46 |
|
47 |
struct jpg_err_mgr |
48 |
{ |
49 |
struct jpeg_error_mgr err; |
50 |
jmp_buf setjmp_buffer; |
51 |
}; |
52 |
|
53 |
static void |
54 |
cv_error_exit (j_common_ptr cinfo) |
55 |
{ |
56 |
cinfo->err->output_message (cinfo); |
57 |
longjmp (((struct jpg_err_mgr *)cinfo->err)->setjmp_buffer, 99); |
58 |
} |
59 |
|
60 |
static void |
61 |
cv_error_output (j_common_ptr cinfo) |
62 |
{ |
63 |
char msg[JMSG_LENGTH_MAX]; |
64 |
|
65 |
cinfo->err->format_message (cinfo, msg); |
66 |
|
67 |
fprintf (stderr, "JPEG decoding error: %s\n", msg); |
68 |
return; |
69 |
} |
70 |
|
71 |
static void |
72 |
rgb_to_hsv (unsigned int r, unsigned int g, unsigned int b, |
73 |
unsigned int *h, unsigned int *s, unsigned int *v) |
74 |
{ |
75 |
unsigned int mx = r; if (g > mx) mx = g; if (b > mx) mx = b; |
76 |
unsigned int mn = r; if (g < mn) mn = g; if (b < mn) mn = b; |
77 |
unsigned int delta = mx - mn; |
78 |
|
79 |
*v = mx; |
80 |
|
81 |
*s = mx ? delta * 255 / mx : 0; |
82 |
|
83 |
if (delta == 0) |
84 |
*h = 0; |
85 |
else |
86 |
{ |
87 |
if (r == mx) |
88 |
*h = ((int)g - (int)b) * 255 / (int)(delta * 3); |
89 |
else if (g == mx) |
90 |
*h = ((int)b - (int)r) * 255 / (int)(delta * 3) + 52; |
91 |
else if (b == mx) |
92 |
*h = ((int)r - (int)g) * 255 / (int)(delta * 3) + 103; |
93 |
|
94 |
*h &= 255; |
95 |
} |
96 |
} |
97 |
|
98 |
struct feature { |
99 |
float v1, v2, v3; // mean, square, cube |
100 |
int n; |
101 |
}; |
102 |
|
103 |
static void |
104 |
feature_init (struct feature *f) |
105 |
{ |
106 |
f->v1 = 0.; |
107 |
f->v2 = 0.; |
108 |
f->v3 = 0.; |
109 |
f->n = 0; |
110 |
} |
111 |
|
112 |
// didn't find an algorithm to neatly do mean, variance and skew in one pass. |
113 |
// elmex ist schuld. |
114 |
static void |
115 |
feature_update_pass_1 (struct feature *f, unsigned int v) |
116 |
{ |
117 |
f->v1 += v; |
118 |
f->n += 1; |
119 |
} |
120 |
|
121 |
static void |
122 |
feature_finish_pass_1 (struct feature *f) |
123 |
{ |
124 |
if (f->n < 1) |
125 |
return; |
126 |
|
127 |
f->v1 /= f->n; |
128 |
} |
129 |
|
130 |
static void |
131 |
feature_update_pass_2 (struct feature *f, unsigned int v) |
132 |
{ |
133 |
float d = v - f->v1; |
134 |
|
135 |
f->v2 += d * d; |
136 |
f->v3 += d * d * d; |
137 |
} |
138 |
|
139 |
static void |
140 |
feature_finish_pass_2 (struct feature *f) |
141 |
{ |
142 |
if (f->n < 1) |
143 |
return; |
144 |
|
145 |
f->v2 /= f->n; |
146 |
f->v3 /= f->n; |
147 |
|
148 |
f->v1 /= 255.; |
149 |
f->v2 /= 255. * 255.; f->v2 = sqrtf (f->v2); |
150 |
f->v3 /= 255. * 255. * 255.; f->v3 = powf (fabsf (f->v3), 1./3.); |
151 |
} |
152 |
|
153 |
static guint32 a85_val; |
154 |
static guint a85_cnt; |
155 |
static guchar a85_buf[LINELENGTH], *a85_ptr; |
156 |
|
157 |
static void |
158 |
a85_init (void) |
159 |
{ |
160 |
a85_cnt = 4; |
161 |
a85_ptr = a85_buf; |
162 |
} |
163 |
|
164 |
static void |
165 |
a85_push (PerlIO *fp, guchar c) |
166 |
{ |
167 |
a85_val = a85_val << 8 | c; |
168 |
|
169 |
if (!--a85_cnt) |
170 |
{ |
171 |
a85_cnt = 4; |
172 |
if (a85_val) |
173 |
{ |
174 |
a85_ptr[4] = (a85_val % 85) + 33; a85_val /= 85; |
175 |
a85_ptr[3] = (a85_val % 85) + 33; a85_val /= 85; |
176 |
a85_ptr[2] = (a85_val % 85) + 33; a85_val /= 85; |
177 |
a85_ptr[1] = (a85_val % 85) + 33; a85_val /= 85; |
178 |
a85_ptr[0] = (a85_val ) + 33; |
179 |
|
180 |
a85_ptr += 5; |
181 |
} |
182 |
else |
183 |
*a85_ptr++ = 'z'; |
184 |
|
185 |
if (a85_ptr >= a85_buf + sizeof (a85_buf) - 7) |
186 |
{ |
187 |
*a85_ptr++ = '\n'; |
188 |
PerlIO_write (fp, a85_buf, a85_ptr - a85_buf); |
189 |
a85_ptr = a85_buf; |
190 |
} |
191 |
} |
192 |
|
193 |
} |
194 |
|
195 |
static void |
196 |
a85_finish (PerlIO *fp) |
197 |
{ |
198 |
while (a85_cnt != 4) |
199 |
a85_push (fp, 0); |
200 |
|
201 |
*a85_ptr++ = '~'; // probably buggy end-marker |
202 |
*a85_ptr++ = '>'; // probably buggy end-marker |
203 |
*a85_ptr++ = '\n'; |
204 |
|
205 |
PerlIO_write (fp, a85_buf, a85_ptr - a85_buf); |
206 |
} |
207 |
|
208 |
///////////////////////////////////////////////////////////////////////////// |
209 |
// memory source for libjpeg |
210 |
|
211 |
static void cv_ms_init (j_decompress_ptr cinfo) |
212 |
{ |
213 |
} |
214 |
|
215 |
static void cv_ms_term (j_decompress_ptr cinfo) |
216 |
{ |
217 |
} |
218 |
|
219 |
static boolean cv_ms_fill (j_decompress_ptr cinfo) |
220 |
{ |
221 |
// unexpected EOF, warn and generate fake EOI marker |
222 |
|
223 |
WARNMS (cinfo, JWRN_JPEG_EOF); |
224 |
|
225 |
struct jpeg_source_mgr *src = (struct jpeg_source_mgr *)cinfo->src; |
226 |
|
227 |
static const JOCTET eoi[] = { 0xFF, JPEG_EOI }; |
228 |
|
229 |
src->next_input_byte = eoi; |
230 |
src->bytes_in_buffer = sizeof (eoi); |
231 |
|
232 |
return TRUE; |
233 |
} |
234 |
|
235 |
static void cv_ms_skip (j_decompress_ptr cinfo, long num_bytes) |
236 |
{ |
237 |
struct jpeg_source_mgr *src = (struct jpeg_source_mgr *)cinfo->src; |
238 |
|
239 |
src->next_input_byte += num_bytes; |
240 |
src->bytes_in_buffer -= num_bytes; |
241 |
} |
242 |
|
243 |
static void cv_jpeg_mem_src (j_decompress_ptr cinfo, void *buf, size_t buflen) |
244 |
{ |
245 |
struct jpeg_source_mgr *src; |
246 |
|
247 |
if (!cinfo->src) |
248 |
cinfo->src = (struct jpeg_source_mgr *) |
249 |
(*cinfo->mem->alloc_small) ( |
250 |
(j_common_ptr) cinfo, JPOOL_PERMANENT, sizeof (struct jpeg_source_mgr) |
251 |
); |
252 |
|
253 |
src = (struct jpeg_source_mgr *)cinfo->src; |
254 |
src->init_source = cv_ms_init; |
255 |
src->fill_input_buffer = cv_ms_fill; |
256 |
src->skip_input_data = cv_ms_skip; |
257 |
src->resync_to_restart = jpeg_resync_to_restart; |
258 |
src->term_source = cv_ms_term; |
259 |
src->next_input_byte = (JOCTET *)buf; |
260 |
src->bytes_in_buffer = buflen; |
261 |
} |
262 |
|
263 |
///////////////////////////////////////////////////////////////////////////// |
264 |
|
265 |
/* great, the jpeg-xl reference implementaton requires us to parse bmff files */ |
266 |
|
267 |
struct bmff_box |
268 |
{ |
269 |
char type[4]; |
270 |
const uint8_t *ptr; |
271 |
size_t size; |
272 |
}; |
273 |
|
274 |
static int |
275 |
bmff_parse_box (struct bmff_box *box, const uint8_t **next_in, size_t *avail_in) |
276 |
{ |
277 |
if (*avail_in < 8) |
278 |
return 0; |
279 |
|
280 |
box->size = ((*next_in)[0] << 24) |
281 |
| ((*next_in)[1] << 16) |
282 |
| ((*next_in)[2] << 8) |
283 |
| ((*next_in)[3] ); |
284 |
|
285 |
if (box->size < 8) |
286 |
return 0; |
287 |
|
288 |
if (*avail_in < box->size) |
289 |
return 0; |
290 |
|
291 |
memcpy (box->type, *next_in + 4, 4); |
292 |
box->ptr = *next_in + 8; |
293 |
|
294 |
*next_in += box->size; |
295 |
*avail_in -= box->size; |
296 |
|
297 |
box->size -= 8; |
298 |
|
299 |
return 1; |
300 |
} |
301 |
|
302 |
///////////////////////////////////////////////////////////////////////////// |
303 |
|
304 |
MODULE = Gtk2::CV PACKAGE = Gtk2::CV |
305 |
|
306 |
PROTOTYPES: ENABLE |
307 |
|
308 |
void |
309 |
_exit (int code) |
310 |
|
311 |
# calculate the common prefix length of two strings |
312 |
# missing function in perl. really :) |
313 |
int |
314 |
common_prefix_length (a, b) |
315 |
unsigned char *a = (unsigned char *)SvPVutf8_nolen ($arg); |
316 |
unsigned char *b = (unsigned char *)SvPVutf8_nolen ($arg); |
317 |
CODE: |
318 |
RETVAL = 0; |
319 |
|
320 |
while (*a == *b && *a) |
321 |
{ |
322 |
RETVAL += (*a & 0xc0) != 0x80; |
323 |
a++, b++; |
324 |
} |
325 |
|
326 |
OUTPUT: |
327 |
RETVAL |
328 |
|
329 |
int |
330 |
common_prefix_length_byte (a, b) |
331 |
unsigned char *a = (unsigned char *)SvPVbyte_nolen ($arg); |
332 |
unsigned char *b = (unsigned char *)SvPVbyte_nolen ($arg); |
333 |
CODE: |
334 |
RETVAL = 0; |
335 |
while (*a == *b && *a) |
336 |
a++, b++, RETVAL++; |
337 |
OUTPUT: |
338 |
RETVAL |
339 |
|
340 |
const char * |
341 |
magic (SV *path_or_data) |
342 |
ALIAS: |
343 |
magic = 0 |
344 |
magic_mime = 1 |
345 |
magic_buffer = 2 |
346 |
magic_buffer_mime = 3 |
347 |
CODE: |
348 |
{ |
349 |
STRLEN len; |
350 |
char *data = SvPVbyte (path_or_data, len); |
351 |
|
352 |
if (!magic_cookie[0]) |
353 |
{ |
354 |
magic_cookie[0] = magic_open (MAGIC_SYMLINK); |
355 |
magic_cookie[1] = magic_open (MAGIC_SYMLINK | MAGIC_MIME_TYPE); |
356 |
magic_load (magic_cookie[0], 0); |
357 |
magic_load (magic_cookie[1], 0); |
358 |
} |
359 |
|
360 |
perlinterp_release (); |
361 |
|
362 |
RETVAL = ix & 2 |
363 |
? magic_buffer (magic_cookie[ix & 1], data, len) |
364 |
: magic_file (magic_cookie[ix & 1], data); |
365 |
|
366 |
perlinterp_acquire (); |
367 |
} |
368 |
OUTPUT: |
369 |
RETVAL |
370 |
|
371 |
# missing/broken in Gtk2 perl module |
372 |
|
373 |
void |
374 |
gdk_window_clear_hints (GdkWindow *window) |
375 |
CODE: |
376 |
gdk_window_set_geometry_hints (window, 0, 0); |
377 |
|
378 |
gboolean |
379 |
gdk_net_wm_supports (GdkAtom property) |
380 |
CODE: |
381 |
#if defined(GDK_WINDOWING_X11) && !defined(GDK_MULTIHEAD_SAFE) |
382 |
RETVAL = gdk_net_wm_supports (property); |
383 |
#else |
384 |
RETVAL = 0; |
385 |
#endif |
386 |
OUTPUT: |
387 |
RETVAL |
388 |
|
389 |
GdkPixbuf_noinc * |
390 |
dealpha_expose (GdkPixbuf *pb) |
391 |
CODE: |
392 |
perlinterp_release (); |
393 |
{ |
394 |
int w = gdk_pixbuf_get_width (pb); |
395 |
int h = gdk_pixbuf_get_height (pb); |
396 |
int bpp = gdk_pixbuf_get_n_channels (pb); |
397 |
int x, y, i; |
398 |
guchar *src = gdk_pixbuf_get_pixels (pb), *dst; |
399 |
int sstr = gdk_pixbuf_get_rowstride (pb), dstr; |
400 |
|
401 |
RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, 0, 8, w, h); |
402 |
|
403 |
dst = gdk_pixbuf_get_pixels (RETVAL); |
404 |
dstr = gdk_pixbuf_get_rowstride (RETVAL); |
405 |
|
406 |
for (x = 0; x < w; x++) |
407 |
for (y = 0; y < h; y++) |
408 |
for (i = 0; i < 3; i++) |
409 |
dst[x * 3 + y * dstr + i] = src[x * bpp + y * sstr + i]; |
410 |
} |
411 |
perlinterp_acquire (); |
412 |
OUTPUT: |
413 |
RETVAL |
414 |
|
415 |
GdkPixbuf_noinc * |
416 |
rotate (GdkPixbuf *pb, int angle) |
417 |
CODE: |
418 |
perlinterp_release (); |
419 |
if (angle < 0) |
420 |
angle += 360; |
421 |
RETVAL = gdk_pixbuf_rotate_simple (pb, angle == 0 ? GDK_PIXBUF_ROTATE_NONE |
422 |
: angle == 90 ? GDK_PIXBUF_ROTATE_COUNTERCLOCKWISE |
423 |
: angle == 180 ? GDK_PIXBUF_ROTATE_UPSIDEDOWN |
424 |
: angle == 270 ? GDK_PIXBUF_ROTATE_CLOCKWISE |
425 |
: angle); |
426 |
perlinterp_acquire (); |
427 |
OUTPUT: |
428 |
RETVAL |
429 |
|
430 |
const char * |
431 |
filetype (SV *image_data) |
432 |
CODE: |
433 |
{ |
434 |
STRLEN data_len; |
435 |
U8 *data = SvPVbyte (image_data, data_len); |
436 |
static const unsigned char jxl_header[] = { |
437 |
0, 0, 0, 0x0c, 0x4a, 0x58, 0x4c, 0x20, 0x0d, 0xa, 0x87, 0x0a |
438 |
}; |
439 |
|
440 |
if (data_len >= 20 |
441 |
&& data[0] == 0xff |
442 |
&& data[1] == 0xd8 |
443 |
&& data[2] == 0xff) |
444 |
RETVAL = "image/jpeg"; |
445 |
else if (data_len >= 12 |
446 |
&& data[ 0] == (U8)'R' |
447 |
&& data[ 1] == (U8)'I' |
448 |
&& data[ 2] == (U8)'F' |
449 |
&& data[ 3] == (U8)'F' |
450 |
&& data[ 8] == (U8)'W' |
451 |
&& data[ 9] == (U8)'E' |
452 |
&& data[10] == (U8)'B' |
453 |
&& data[11] == (U8)'P') |
454 |
RETVAL = "image/webp"; |
455 |
else if (data_len >= 16 |
456 |
&& data[ 0] == 0x89 |
457 |
&& data[ 1] == (U8)'P' |
458 |
&& data[ 2] == (U8)'N' |
459 |
&& data[ 3] == (U8)'G' |
460 |
&& data[ 4] == 0x0d |
461 |
&& data[ 5] == 0x0a |
462 |
&& data[ 6] == 0x1a |
463 |
&& data[ 7] == 0x0a) |
464 |
RETVAL = "image/png"; |
465 |
else if (data_len >= sizeof (jxl_header) && memcmp (data, jxl_header, sizeof (jxl_header)) == 0) |
466 |
RETVAL = "image/jxl"; // todo: might want to use JxlSignatureCheck |
467 |
else if (data_len >= 2 |
468 |
&& data[0] == 0xff |
469 |
&& data[1] == 0x0a) |
470 |
RETVAL = "image/jxl"; |
471 |
else if (data_len >= 13 |
472 |
&& data[0] == 'G' |
473 |
&& data[1] == 'I' |
474 |
&& data[2] == 'F' |
475 |
&& data[3] == '8' |
476 |
//&& (data[4] == '7' || data[4] == '9') |
477 |
&& data[5] == 'a') |
478 |
{ |
479 |
RETVAL = "image/gif"; |
480 |
|
481 |
// now see if its animated - we require the netscape application header for this |
482 |
int ofs = 13; |
483 |
|
484 |
if (data[10] & 0x80) |
485 |
ofs += (1 << ((data[10] & 7) + 1)) * 3; |
486 |
|
487 |
if (data_len >= ofs + 2 + 1 + 11) |
488 |
{ |
489 |
|
490 |
// skip a graphic control extension block. we assume |
491 |
// there is at most one such block - while the NAB |
492 |
// has to come firstz, some files do not obey this |
493 |
if (data[ofs] == 0x21 && data[ofs + 1] == 0xf9) |
494 |
ofs += 3 + data[ofs + 2] + 1; |
495 |
|
496 |
if (data_len >= ofs + 2 + 1 + 11) |
497 |
if (!memcmp (data + ofs, "\x21\xff\x0bNETSCAPE2.0", sizeof ("\x21\xff\x0bNETSCAPE2.0") - 1)) |
498 |
RETVAL = "video/gif"; |
499 |
} |
500 |
} |
501 |
else if (data_len >= 0x8000 + 6 |
502 |
&& data[0x8000+1] == (U8)'B' |
503 |
&& data[0x8000+2] == (U8)'E' |
504 |
&& data[0x8000+3] == (U8)'A' |
505 |
&& data[0x8000+4] == (U8)'0' |
506 |
&& data[0x8000+5] == (U8)'1') |
507 |
RETVAL = "video/iso-bluray"; |
508 |
else if (data_len >= 0x8000 + 6 |
509 |
&& data[0x8000+1] == (U8)'C' |
510 |
&& data[0x8000+2] == (U8)'D' |
511 |
&& data[0x8000+3] == (U8)'0' |
512 |
&& data[0x8000+4] == (U8)'0' |
513 |
&& data[0x8000+5] == (U8)'1') |
514 |
RETVAL = "video/iso9660"; |
515 |
|
516 |
else |
517 |
XSRETURN_UNDEF; |
518 |
} |
519 |
OUTPUT: |
520 |
RETVAL |
521 |
|
522 |
GdkPixbuf_noinc * |
523 |
decode_webp (SV *image_data, int thumbnail = 0, int iw = 0, int ih = 0) |
524 |
CODE: |
525 |
{ |
526 |
#if WEBP |
527 |
STRLEN data_size; |
528 |
int alpha; |
529 |
WebPData data; |
530 |
WebPDemuxer *demux; |
531 |
WebPIterator iter; |
532 |
WebPDecoderConfig config; |
533 |
int inw, inh; |
534 |
|
535 |
data.bytes = (uint8_t *)SvPVbyte (image_data, data_size); |
536 |
data.size = data_size; |
537 |
|
538 |
perlinterp_release (); |
539 |
|
540 |
RETVAL = 0; |
541 |
|
542 |
if (!(demux = WebPDemux (&data))) |
543 |
goto err_demux; |
544 |
|
545 |
if (!WebPDemuxGetFrame (demux, 1, &iter)) |
546 |
goto err_iter; |
547 |
|
548 |
if (!WebPInitDecoderConfig (&config)) |
549 |
goto err_iter; |
550 |
|
551 |
config.options.use_threads = 1; |
552 |
|
553 |
if (WebPGetFeatures (iter.fragment.bytes, iter.fragment.size, &config.input) != VP8_STATUS_OK) |
554 |
goto err_iter; |
555 |
|
556 |
inw = config.input.width; |
557 |
inh = config.input.height; |
558 |
|
559 |
if (thumbnail) |
560 |
{ |
561 |
if (inw * ih > inh * iw) |
562 |
ih = (iw * inh + inw - 1) / inw; |
563 |
else |
564 |
iw = (ih * inw + inh - 1) / inh; |
565 |
|
566 |
config.options.bypass_filtering = 1; |
567 |
config.options.no_fancy_upsampling = 1; |
568 |
|
569 |
config.options.use_scaling = 1; |
570 |
config.options.scaled_width = iw; |
571 |
config.options.scaled_height = ih; |
572 |
} |
573 |
else |
574 |
{ |
575 |
iw = inw; |
576 |
ih = inh; |
577 |
} |
578 |
|
579 |
alpha = !!config.input.has_alpha; |
580 |
|
581 |
RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, alpha, 8, iw, ih); |
582 |
if (!RETVAL) |
583 |
goto err_iter; |
584 |
|
585 |
config.output.colorspace = alpha ? MODE_RGBA : MODE_RGB; |
586 |
config.output.u.RGBA.rgba = gdk_pixbuf_get_pixels (RETVAL); |
587 |
config.output.u.RGBA.stride = gdk_pixbuf_get_rowstride (RETVAL); |
588 |
config.output.u.RGBA.size = gdk_pixbuf_get_byte_length (RETVAL); |
589 |
config.output.is_external_memory = 1; |
590 |
|
591 |
if (WebPDecode (iter.fragment.bytes, iter.fragment.size, &config) != VP8_STATUS_OK) |
592 |
{ |
593 |
g_object_unref (RETVAL); |
594 |
RETVAL = 0; |
595 |
goto err_iter; |
596 |
} |
597 |
|
598 |
err_iter: |
599 |
WebPDemuxReleaseIterator (&iter); |
600 |
err_demux: |
601 |
WebPDemuxDelete (demux); |
602 |
|
603 |
perlinterp_acquire (); |
604 |
#else |
605 |
croak ("load_webp: webp not enabled at compile time"); |
606 |
#endif |
607 |
} |
608 |
OUTPUT: |
609 |
RETVAL |
610 |
|
611 |
GdkPixbuf_noinc * |
612 |
decode_jxl (SV *image_data, int thumbnail = 0, int iw = 0, int ih = 0) |
613 |
CODE: |
614 |
{ |
615 |
#if JXL |
616 |
JxlDecoder *dec; |
617 |
JxlBasicInfo info; |
618 |
const uint8_t *next_in = (uint8_t *)SvPVbyte_nolen (image_data); |
619 |
size_t avail_in = SvCUR (image_data); |
620 |
const char *error = 0; |
621 |
JxlDecoderStatus status; |
622 |
static void *runner_cache; |
623 |
void *runner = 0; |
624 |
|
625 |
RETVAL = 0; |
626 |
|
627 |
if (runner_cache) |
628 |
runner = runner_cache; |
629 |
else |
630 |
runner = JxlThreadParallelRunnerCreate (0, JxlThreadParallelRunnerDefaultNumWorkerThreads ()); |
631 |
|
632 |
runner_cache = 0; |
633 |
|
634 |
perlinterp_release (); |
635 |
|
636 |
dec = JxlDecoderCreate (0); |
637 |
|
638 |
error = "JxlDecoderCreate failed"; |
639 |
if (!dec) |
640 |
goto done; |
641 |
|
642 |
status = JxlDecoderSetParallelRunner (dec, JxlThreadParallelRunner, runner); |
643 |
error = "JxlDecoderSetParallelRunner failed"; |
644 |
if (status != JXL_DEC_SUCCESS) |
645 |
goto done; |
646 |
|
647 |
error = "JxlDecoderSubscribeEvents failed"; |
648 |
status = JxlDecoderSubscribeEvents (dec, JXL_DEC_BASIC_INFO | JXL_DEC_FULL_IMAGE); |
649 |
if (status != JXL_DEC_SUCCESS) |
650 |
goto done; |
651 |
|
652 |
status = JxlDecoderSetInput (dec, next_in, avail_in); |
653 |
error = "JxlDecoderSetInput failed"; |
654 |
if (status != JXL_DEC_SUCCESS) |
655 |
goto done; |
656 |
|
657 |
for (;;) |
658 |
{ |
659 |
status = JxlDecoderProcessInput (dec); |
660 |
|
661 |
switch (status) |
662 |
{ |
663 |
case JXL_DEC_FULL_IMAGE: |
664 |
error = 0; |
665 |
goto done; |
666 |
|
667 |
case JXL_DEC_ERROR: |
668 |
error = "JxlDecoderProcessInput failed"; |
669 |
goto done; |
670 |
|
671 |
case JXL_DEC_NEED_MORE_INPUT: |
672 |
error = "incomplete file"; |
673 |
goto done; |
674 |
|
675 |
case JXL_DEC_SUCCESS: |
676 |
error = "incomplete decode"; |
677 |
goto done; |
678 |
|
679 |
case JXL_DEC_BASIC_INFO: |
680 |
{ |
681 |
status = JxlDecoderGetBasicInfo (dec, &info); |
682 |
error = "JxlDecoderGetBasicInfo failed"; |
683 |
if (status != JXL_DEC_SUCCESS) |
684 |
goto done; |
685 |
|
686 |
RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, !!info.alpha_bits, 8, info.xsize, info.ysize); |
687 |
error = "unable to allocate pixbuf"; |
688 |
if (!RETVAL) |
689 |
goto done; |
690 |
|
691 |
JxlPixelFormat format = { |
692 |
info.alpha_bits ? 4 : 3, |
693 |
JXL_TYPE_UINT8, |
694 |
JXL_NATIVE_ENDIAN, |
695 |
gdk_pixbuf_get_rowstride (RETVAL) |
696 |
}; |
697 |
|
698 |
// cannot use gdk_pixbuf_get_byte_length because that does |
699 |
// not return the size of the buffer, but the size of the buffer without |
700 |
// the last padding bytes. the internal buffer is rowstride * ysize, |
701 |
// and this is what the jxl decoder needs. none of this is documented |
702 |
// in either library, of course. |
703 |
|
704 |
status = JxlDecoderSetImageOutBuffer ( |
705 |
dec, |
706 |
&format, |
707 |
gdk_pixbuf_get_pixels (RETVAL), |
708 |
gdk_pixbuf_get_rowstride (RETVAL) * info.ysize |
709 |
); |
710 |
error = "JxlDecoderSetImageOutBuffer failed"; |
711 |
if (status != JXL_DEC_SUCCESS) |
712 |
goto done; |
713 |
} |
714 |
break; |
715 |
|
716 |
default: |
717 |
error = "unexpected event"; |
718 |
goto done; |
719 |
} |
720 |
} |
721 |
|
722 |
done: |
723 |
if (dec) |
724 |
JxlDecoderDestroy (dec); |
725 |
|
726 |
perlinterp_acquire (); |
727 |
|
728 |
if (runner_cache) |
729 |
JxlThreadParallelRunnerDestroy (runner); |
730 |
|
731 |
runner_cache = runner; |
732 |
|
733 |
if (error) |
734 |
{ |
735 |
if (RETVAL) |
736 |
g_object_unref (RETVAL); |
737 |
|
738 |
croak ("load_jxl: %s (status %d)", error, status); |
739 |
} |
740 |
#else |
741 |
croak ("load_jxl: jpeg-xl not enabled at compile time"); |
742 |
#endif |
743 |
} |
744 |
OUTPUT: |
745 |
RETVAL |
746 |
|
747 |
GdkPixbuf_noinc * |
748 |
decode_jpeg (SV *image_data, int thumbnail = 0, int iw = 0, int ih = 0) |
749 |
CODE: |
750 |
{ |
751 |
struct jpeg_decompress_struct cinfo; |
752 |
struct jpg_err_mgr jerr; |
753 |
int rs; |
754 |
volatile GdkPixbuf *pb = 0; |
755 |
STRLEN data_len; |
756 |
guchar *data = SvPVbyte (image_data, data_len); |
757 |
|
758 |
RETVAL = 0; |
759 |
|
760 |
perlinterp_release (); |
761 |
|
762 |
cinfo.err = jpeg_std_error (&jerr.err); |
763 |
|
764 |
jerr.err.error_exit = cv_error_exit; |
765 |
jerr.err.output_message = cv_error_output; |
766 |
|
767 |
if ((rs = setjmp (jerr.setjmp_buffer))) |
768 |
{ |
769 |
jpeg_destroy_decompress (&cinfo); |
770 |
|
771 |
if (pb) |
772 |
g_object_unref ((gpointer)pb); |
773 |
|
774 |
perlinterp_acquire (); |
775 |
XSRETURN_UNDEF; |
776 |
} |
777 |
|
778 |
if (!data_len) |
779 |
longjmp (jerr.setjmp_buffer, 4); |
780 |
|
781 |
jpeg_create_decompress (&cinfo); |
782 |
cv_jpeg_mem_src (&cinfo, data, data_len); |
783 |
|
784 |
jpeg_read_header (&cinfo, TRUE); |
785 |
|
786 |
cinfo.dct_method = JDCT_DEFAULT; |
787 |
cinfo.do_fancy_upsampling = FALSE; /* worse quality, but nobody compained so far, and gdk-pixbuf does the same */ |
788 |
cinfo.do_block_smoothing = FALSE; |
789 |
cinfo.out_color_space = JCS_RGB; |
790 |
cinfo.quantize_colors = FALSE; |
791 |
|
792 |
cinfo.scale_num = 1; |
793 |
cinfo.scale_denom = 1; |
794 |
|
795 |
jpeg_calc_output_dimensions (&cinfo); |
796 |
|
797 |
if (thumbnail) |
798 |
{ |
799 |
cinfo.dct_method = JDCT_FASTEST; |
800 |
cinfo.do_fancy_upsampling = FALSE; |
801 |
|
802 |
while (cinfo.scale_denom < 8 |
803 |
&& cinfo.output_width >= iw*4 |
804 |
&& cinfo.output_height >= ih*4) |
805 |
{ |
806 |
cinfo.scale_denom <<= 1; |
807 |
jpeg_calc_output_dimensions (&cinfo); |
808 |
} |
809 |
} |
810 |
|
811 |
if (cinfo.output_components != 3) |
812 |
longjmp (jerr.setjmp_buffer, 3); |
813 |
|
814 |
if (cinfo.jpeg_color_space == JCS_YCCK || cinfo.jpeg_color_space == JCS_CMYK) |
815 |
{ |
816 |
cinfo.out_color_space = JCS_CMYK; |
817 |
cinfo.output_components = 4; |
818 |
} |
819 |
|
820 |
pb = RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, cinfo.output_components == 4, 8, cinfo.output_width, cinfo.output_height); |
821 |
if (!RETVAL) |
822 |
longjmp (jerr.setjmp_buffer, 2); |
823 |
|
824 |
data = gdk_pixbuf_get_pixels (RETVAL); |
825 |
rs = gdk_pixbuf_get_rowstride (RETVAL); |
826 |
|
827 |
jpeg_start_decompress (&cinfo); |
828 |
|
829 |
while (cinfo.output_scanline < cinfo.output_height) |
830 |
{ |
831 |
int remaining = cinfo.output_height - cinfo.output_scanline; |
832 |
JSAMPROW rp[4]; |
833 |
|
834 |
rp [0] = data + cinfo.output_scanline * rs; |
835 |
rp [1] = (guchar *)rp [0] + rs; |
836 |
rp [2] = (guchar *)rp [1] + rs; |
837 |
rp [3] = (guchar *)rp [2] + rs; |
838 |
|
839 |
jpeg_read_scanlines (&cinfo, rp, remaining < 4 ? remaining : 4); |
840 |
} |
841 |
|
842 |
if (cinfo.out_color_space == JCS_CMYK) |
843 |
{ |
844 |
guchar *end = data + cinfo.output_height * rs; |
845 |
|
846 |
while (data < end) |
847 |
{ |
848 |
U32 c = data [0]; |
849 |
U32 m = data [1]; |
850 |
U32 y = data [2]; |
851 |
U32 k = data [3]; |
852 |
|
853 |
if (0) |
854 |
if (cinfo.Adobe_transform == 2) |
855 |
{ |
856 |
c ^= 0xff; |
857 |
m ^= 0xff; |
858 |
y ^= 0xff; |
859 |
k ^= 0xff; |
860 |
} |
861 |
|
862 |
data [0] = (c * k + 0x80) / 0xff; |
863 |
data [1] = (m * k + 0x80) / 0xff; |
864 |
data [2] = (y * k + 0x80) / 0xff; |
865 |
data [3] = 0xff; |
866 |
|
867 |
data += 4; |
868 |
} |
869 |
} |
870 |
|
871 |
jpeg_finish_decompress (&cinfo); |
872 |
jpeg_destroy_decompress (&cinfo); |
873 |
perlinterp_acquire (); |
874 |
} |
875 |
OUTPUT: |
876 |
RETVAL |
877 |
|
878 |
void |
879 |
compare (GdkPixbuf *a, GdkPixbuf *b) |
880 |
PPCODE: |
881 |
perlinterp_release (); |
882 |
{ |
883 |
int w = gdk_pixbuf_get_width (a); |
884 |
int h = gdk_pixbuf_get_height (a); |
885 |
|
886 |
int sa = gdk_pixbuf_get_rowstride (a); |
887 |
int sb = gdk_pixbuf_get_rowstride (b); |
888 |
|
889 |
guchar *pa = gdk_pixbuf_get_pixels (a); |
890 |
guchar *pb = gdk_pixbuf_get_pixels (b); |
891 |
|
892 |
int x, y; |
893 |
|
894 |
assert (w == gdk_pixbuf_get_width (b)); |
895 |
assert (h == gdk_pixbuf_get_height (b)); |
896 |
|
897 |
assert (gdk_pixbuf_get_n_channels (a) == 3); |
898 |
assert (gdk_pixbuf_get_n_channels (b) == 3); |
899 |
|
900 |
double diff = 0.; |
901 |
int peak = 0; |
902 |
|
903 |
if (w && h) |
904 |
for (y = 0; y < h; y++) |
905 |
{ |
906 |
guchar *pa_ = pa + y * sa; |
907 |
guchar *pb_ = pb + y * sb; |
908 |
|
909 |
for (x = 0; x < w; x++) |
910 |
{ |
911 |
int d; |
912 |
|
913 |
d = ((int)*pa_++) - ((int)*pb_++); diff += d*d; peak = MAX (peak, abs (d)); |
914 |
d = ((int)*pa_++) - ((int)*pb_++); diff += d*d; peak = MAX (peak, abs (d)); |
915 |
d = ((int)*pa_++) - ((int)*pb_++); diff += d*d; peak = MAX (peak, abs (d)); |
916 |
} |
917 |
} |
918 |
|
919 |
perlinterp_acquire (); |
920 |
|
921 |
EXTEND (SP, 2); |
922 |
PUSHs (sv_2mortal (newSVnv (sqrt (diff / (w * h * 3. * 255. * 255.))))); |
923 |
PUSHs (sv_2mortal (newSVnv (peak / 255.))); |
924 |
} |
925 |
|
926 |
############################################################################# |
927 |
|
928 |
MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer |
929 |
|
930 |
# currently only works for filenames (octet strings) |
931 |
|
932 |
SV * |
933 |
foldcase (SV *pathsv) |
934 |
PROTOTYPE: $ |
935 |
CODE: |
936 |
{ |
937 |
STRLEN plen; |
938 |
U8 *path = (U8 *)SvPV (pathsv, plen); |
939 |
U8 *pend = path + plen; |
940 |
U8 dst [plen * 8 * 3], *dstp = dst; |
941 |
|
942 |
while (path < pend) |
943 |
{ |
944 |
U8 ch = *path; |
945 |
|
946 |
if (ch >= 'a' && ch <= 'z') |
947 |
*dstp++ = *path++; |
948 |
else if (ch >= 'A' && ch <= 'Z') |
949 |
*dstp++ = *path++ + ('a' - 'A'); |
950 |
else if (ch >= '0' && ch <= '9') |
951 |
{ |
952 |
/* version sort, up to 8 digits */ |
953 |
STRLEN el, nl = 0; |
954 |
while (*path >= '0' && *path <= '9' && path < pend) |
955 |
path++, nl++; |
956 |
|
957 |
for (el = nl; el < 8; el++) |
958 |
*dstp++ = '0'; |
959 |
|
960 |
memcpy (dstp, path - nl, nl); |
961 |
dstp += nl; |
962 |
} |
963 |
else |
964 |
*dstp++ = *path++; |
965 |
#if 0 |
966 |
else |
967 |
{ |
968 |
STRLEN cl; |
969 |
to_utf8_fold (path, dstp, &cl); |
970 |
dstp += cl; |
971 |
path += is_utf8_char (path); |
972 |
} |
973 |
#endif |
974 |
} |
975 |
|
976 |
RETVAL = newSVpvn ((const char *)dst, dstp - dst); |
977 |
} |
978 |
OUTPUT: |
979 |
RETVAL |
980 |
|
981 |
GdkPixbuf_noinc * |
982 |
p7_to_pb (int w, int h, SV *src_sv) |
983 |
PROTOTYPE: @ |
984 |
CODE: |
985 |
{ |
986 |
int x, y; |
987 |
guchar *dst, *d; |
988 |
int dstr; |
989 |
guchar *src = (guchar *)SvPVbyte_nolen (src_sv); |
990 |
|
991 |
RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, 0, 8, w, h); |
992 |
dst = gdk_pixbuf_get_pixels (RETVAL); |
993 |
dstr = gdk_pixbuf_get_rowstride (RETVAL); |
994 |
|
995 |
for (y = 0; y < h; y++) |
996 |
for (d = dst + y * dstr, x = 0; x < w; x++) |
997 |
{ |
998 |
*d++ = (((*src >> 5) & 7) * 255 + 4) / 7; |
999 |
*d++ = (((*src >> 2) & 7) * 255 + 4) / 7; |
1000 |
*d++ = (((*src >> 0) & 3) * 255 + 2) / 3; |
1001 |
|
1002 |
src++; |
1003 |
} |
1004 |
} |
1005 |
OUTPUT: |
1006 |
RETVAL |
1007 |
|
1008 |
############################################################################# |
1009 |
|
1010 |
MODULE = Gtk2::CV PACKAGE = Gtk2::CV::PostScript |
1011 |
|
1012 |
void |
1013 |
dump_ascii85 (PerlIO *fp, GdkPixbuf *pb) |
1014 |
CODE: |
1015 |
{ |
1016 |
int w = gdk_pixbuf_get_width (pb); |
1017 |
int h = gdk_pixbuf_get_height (pb); |
1018 |
int x, y, i; |
1019 |
guchar *dst; |
1020 |
int bpp = gdk_pixbuf_get_n_channels (pb); |
1021 |
guchar *src = gdk_pixbuf_get_pixels (pb); |
1022 |
int sstr = gdk_pixbuf_get_rowstride (pb); |
1023 |
|
1024 |
a85_init (); |
1025 |
|
1026 |
for (y = 0; y < h; y++) |
1027 |
for (x = 0; x < w; x++) |
1028 |
for (i = 0; i < (bpp < 3 ? 1 : 3); i++) |
1029 |
a85_push (fp, src [x * bpp + y * sstr + i]); |
1030 |
|
1031 |
a85_finish (fp); |
1032 |
} |
1033 |
|
1034 |
void |
1035 |
dump_binary (PerlIO *fp, GdkPixbuf *pb) |
1036 |
CODE: |
1037 |
{ |
1038 |
int w = gdk_pixbuf_get_width (pb); |
1039 |
int h = gdk_pixbuf_get_height (pb); |
1040 |
int x, y, i; |
1041 |
guchar *dst; |
1042 |
int bpp = gdk_pixbuf_get_n_channels (pb); |
1043 |
guchar *src = gdk_pixbuf_get_pixels (pb); |
1044 |
int sstr = gdk_pixbuf_get_rowstride (pb); |
1045 |
|
1046 |
for (y = 0; y < h; y++) |
1047 |
for (x = 0; x < w; x++) |
1048 |
for (i = 0; i < (bpp < 3 ? 1 : 3); i++) |
1049 |
PerlIO_putc (fp, src [x * bpp + y * sstr + i]); |
1050 |
} |
1051 |
|
1052 |
############################################################################# |
1053 |
|
1054 |
MODULE = Gtk2::CV PACKAGE = Gtk2::CV |
1055 |
|
1056 |
SV * |
1057 |
pb_to_hv84 (GdkPixbuf *pb) |
1058 |
CODE: |
1059 |
{ |
1060 |
int w = gdk_pixbuf_get_width (pb); |
1061 |
int h = gdk_pixbuf_get_height (pb); |
1062 |
int x, y; |
1063 |
guchar *dst; |
1064 |
int bpp = gdk_pixbuf_get_n_channels (pb); |
1065 |
guchar *src = gdk_pixbuf_get_pixels (pb); |
1066 |
int sstr = gdk_pixbuf_get_rowstride (pb); |
1067 |
|
1068 |
RETVAL = newSV (6 * 8 * 12 / 8); |
1069 |
SvPOK_only (RETVAL); |
1070 |
SvCUR_set (RETVAL, 6 * 8 * 12 / 8); |
1071 |
|
1072 |
dst = (guchar *)SvPVX (RETVAL); |
1073 |
|
1074 |
/* some primitive error distribution + random dithering */ |
1075 |
|
1076 |
for (y = 0; y < h; y++) |
1077 |
{ |
1078 |
guchar *p = src + y * sstr; |
1079 |
|
1080 |
for (x = 0; x < w; x += 2) |
1081 |
{ |
1082 |
unsigned int r, g, b, h, s, v, H, V1, V2; |
1083 |
|
1084 |
if (bpp == 3) |
1085 |
r = *p++, g = *p++, b = *p++; |
1086 |
else if (bpp == 1) |
1087 |
r = g = b = *p++; |
1088 |
else |
1089 |
abort (); |
1090 |
|
1091 |
rgb_to_hsv (r, g, b, &h, &s, &v); |
1092 |
|
1093 |
H = (h * 15 / 255) << 4; |
1094 |
V1 = v; |
1095 |
|
1096 |
if (bpp == 3) |
1097 |
r = *p++, g = *p++, b = *p++; |
1098 |
else if (bpp == 1) |
1099 |
r = g = b = *p++; |
1100 |
else |
1101 |
abort (); |
1102 |
|
1103 |
rgb_to_hsv (r, g, b, &h, &s, &v); |
1104 |
|
1105 |
H |= h * 15 / 255; |
1106 |
V2 = v; |
1107 |
|
1108 |
*dst++ = H; |
1109 |
*dst++ = V1; |
1110 |
*dst++ = V2; |
1111 |
} |
1112 |
} |
1113 |
} |
1114 |
OUTPUT: |
1115 |
RETVAL |
1116 |
|
1117 |
SV * |
1118 |
hv84_to_av (unsigned char *hv84) |
1119 |
CODE: |
1120 |
{ |
1121 |
int i = 72 / 3; |
1122 |
AV *av = newAV (); |
1123 |
|
1124 |
RETVAL = (SV *)newRV_noinc ((SV *)av); |
1125 |
while (i--) |
1126 |
{ |
1127 |
int h = *hv84++; |
1128 |
int v1 = *hv84++; |
1129 |
int v2 = *hv84++; |
1130 |
|
1131 |
av_push (av, newSViv (v1)); |
1132 |
av_push (av, newSViv ((h >> 4) * 255 / 15)); |
1133 |
av_push (av, newSViv (v2)); |
1134 |
av_push (av, newSViv ((h & 15) * 255 / 15)); |
1135 |
} |
1136 |
} |
1137 |
OUTPUT: |
1138 |
RETVAL |
1139 |
|
1140 |
############################################################################# |
1141 |
|
1142 |
MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Plugin::RCluster |
1143 |
|
1144 |
SV * |
1145 |
extract_features (SV *ar) |
1146 |
CODE: |
1147 |
{ |
1148 |
int i; |
1149 |
AV *av, *result; |
1150 |
|
1151 |
if (!SvROK (ar) || SvTYPE (SvRV (ar)) != SVt_PVAV) |
1152 |
croak ("Not an array ref as first argument to extract_features"); |
1153 |
|
1154 |
av = (AV *) SvRV (ar); |
1155 |
result = newAV (); |
1156 |
|
1157 |
for (i = 0; i <= av_len (av); ++i) |
1158 |
{ |
1159 |
SV *sv = *av_fetch (av, i, 1); |
1160 |
SV *histsv = newSV (9 * sizeof (float) + 1); |
1161 |
|
1162 |
SvPOK_on (histsv); |
1163 |
SvCUR_set (histsv, 9 * sizeof (float)); |
1164 |
float *hist = (float *)SvPVX (histsv); |
1165 |
|
1166 |
struct feature f_h, f_s, f_v; |
1167 |
feature_init (&f_h); |
1168 |
feature_init (&f_s); |
1169 |
feature_init (&f_v); |
1170 |
|
1171 |
{ |
1172 |
STRLEN len; |
1173 |
unsigned char *buf = (unsigned char *)SvPVbyte (sv, len); |
1174 |
while (len >= 3) |
1175 |
{ |
1176 |
unsigned int r, g, b, h, s, v; |
1177 |
r = *buf++; g = *buf++; b = *buf++; |
1178 |
rgb_to_hsv (r, g, b, &h, &s, &v); |
1179 |
|
1180 |
feature_update_pass_1 (&f_h, h); |
1181 |
feature_update_pass_1 (&f_s, s); |
1182 |
feature_update_pass_1 (&f_v, v); |
1183 |
|
1184 |
len -= 3; |
1185 |
} |
1186 |
|
1187 |
feature_finish_pass_1 (&f_h); |
1188 |
feature_finish_pass_1 (&f_s); |
1189 |
feature_finish_pass_1 (&f_v); |
1190 |
} |
1191 |
|
1192 |
{ |
1193 |
STRLEN len; |
1194 |
unsigned char *buf = (unsigned char *)SvPVbyte (sv, len); |
1195 |
while (len >= 3) |
1196 |
{ |
1197 |
unsigned int r, g, b, h, s, v; |
1198 |
r = *buf++; g = *buf++; b = *buf++; |
1199 |
rgb_to_hsv (r, g, b, &h, &s, &v); |
1200 |
|
1201 |
feature_update_pass_2 (&f_h, h); |
1202 |
feature_update_pass_2 (&f_s, s); |
1203 |
feature_update_pass_2 (&f_v, v); |
1204 |
|
1205 |
len -= 3; |
1206 |
} |
1207 |
|
1208 |
feature_finish_pass_2 (&f_h); |
1209 |
feature_finish_pass_2 (&f_s); |
1210 |
feature_finish_pass_2 (&f_v); |
1211 |
} |
1212 |
|
1213 |
hist [0] = f_h.v1 * 2.; hist [1] = f_h.v2 * 2.; hist [2] = f_h.v3 * 2.; |
1214 |
hist [3] = f_s.v1 ; hist [4] = f_s.v2 ; hist [5] = f_s.v3 ; |
1215 |
hist [6] = f_v.v1 * .5; hist [7] = f_v.v2 * .5; hist [8] = f_v.v3 * .5; |
1216 |
|
1217 |
av_push (result, histsv); |
1218 |
} |
1219 |
|
1220 |
RETVAL = newRV_noinc ((SV *)result); |
1221 |
} |
1222 |
OUTPUT: |
1223 |
RETVAL |
1224 |
|