Here is a little test script: I encountered two problems.First, changing the interpolation method via ->set breaks interpolate, apparently unsetting x and y.
Second, PDL::Primitive::interpolate breaks when expected to thread over y. Best Ingo On 04/17/2012 06:55 PM, Ingo Schmid wrote:
Hi,I willl provide an example tomorrow, it was very busy today. I thought the info of the PDL dims were showing what the problem is.Also, I tried reading the slatec module's doc, not comprehending much. Then I luckily found the Func module. I admit I have difficulties in reading signatures, although successfully working with PDL for years now.Anyways, your replies are greatly appreciated. PDL v. 2.4.10 Ingo PS: $ perl -V Summary of my perl5 (revision 5 version 12 subversion 4) configuration: Platform: osname=linux, osvers=3.1.5-gentoo, archname=x86_64-linux-thread-multiuname='linux riococo 3.1.5-gentoo #1 smp mon dec 12 16:41:23 cet 2011 x86_64 intel(r) xeon(r) cpu w3520 @ 2.67ghz genuineintel gnulinux ' config_args='-des -Duseshrplib -Darchname=x86_64-linux-thread -Dcc=x86_64-pc-linux-gnu-gcc -Doptimize=-O2 -pipe -march=core2 -fomit-frame-pointer -msse4 -msse4.1 -msse4.2 -mcx16 -msahf -Dldflags=-Wl,-O1 -Wl,--as-needed -Dprefix=/usr -Dinstallprefix=/usr -Dsiteprefix=/usr -Dvendorprefix=/usr -Dscriptdir=/usr/bin -Dprivlib=/usr/lib64/perl5/5.12.4 -Darchlib=/usr/lib64/perl5/5.12.4/x86_64-linux-thread-multi -Dsitelib=/usr/lib64/perl5/site_perl/5.12.4 -Dsitearch=/usr/lib64/perl5/site_perl/5.12.4/x86_64-linux-thread-multi -Dvendorlib=/usr/lib64/perl5/vendor_perl/5.12.4 -Dvendorarch=/usr/lib64/perl5/vendor_perl/5.12.4/x86_64-linux-thread-multi -Dman1dir=/usr/share/man/man1 -Dman3dir=/usr/share/man/man3 -Dsiteman1dir=/usr/share/man/man1 -Dsiteman3dir=/usr/share/man/man3 -Dvendorman1dir=/usr/share/man/man1 -Dvendorman3dir=/usr/share/man/man3 -Dman1ext=1 -Dman3ext=3pm -Dlibperl=libperl.so.5.12.4 -Dlocincpth=/usr/include -Dglibpth=/lib64 /usr/lib64 -Duselargefiles -Dd_semctl_semun -Dcf_by=Gentoo -Dmyhostname=localhost -Dperladmin=root@localhost -Dinstallusrbinperl=n -Ud_csh -Uusenm -Di_ndbm -Di_gdbm -Di_db -Dusethreads -DDEBUGGING=none -Dinc_version_list=5.12.3/x86_64-linux-thread-multi 5.12.3 5.12.2/x86_64-linux-thread-multi 5.12.2 5.12.1/x86_64-linux-thread-multi 5.12.1 5.12.0/x86_64-linux-thread-multi 5.12.0 -Dlibpth=/usr/local/lib64 /lib64 /usr/lib64 -Dnoextensions=ODBM_File'hint=recommended, useposix=true, d_sigaction=define useithreads=define, usemultiplicity=define useperlio=define, d_sfio=undef, uselargefiles=define, usesocks=undef use64bitint=define, use64bitall=define, uselongdouble=undef usemymalloc=n, bincompat5005=undef Compiler:cc='x86_64-pc-linux-gnu-gcc', ccflags ='-D_REENTRANT -D_GNU_SOURCE -fno-strict-aliasing -pipe -D_LARGEFILE_SOURCE -D_FILE_OFFSET_BITS=64', optimize='-O2 -pipe -march=core2 -fomit-frame-pointer -msse4 -msse4.1 -msse4.2 -mcx16 -msahf',cppflags='-D_REENTRANT -D_GNU_SOURCE -fno-strict-aliasing -pipe' ccversion='', gccversion='4.5.3', gccosandvers='' intsize=4, longsize=8, ptrsize=8, doublesize=8, byteorder=12345678 d_longlong=define, longlongsize=8, d_longdbl=define, longdblsize=16ivtype='long', ivsize=8, nvtype='double', nvsize=8, Off_t='off_t', lseeksize=8alignbytes=8, prototype=define Linker and Libraries: ld='x86_64-pc-linux-gnu-gcc', ldflags ='-Wl,-O1 -Wl,--as-needed' libpth=/usr/local/lib64 /lib64 /usr/lib64libs=-lnsl -lgdbm -ldb -ldl -lm -lcrypt -lutil -lpthread -lc -lgdbm_compatperllibs=-lnsl -ldl -lm -lcrypt -lutil -lpthread -lclibc=/lib/libc-2.14.1.so, so=so, useshrplib=true, libperl=libperl.so.5.12.4gnulibc_version='2.14.1' Dynamic Linking: dlsrc=dl_dlopen.xs, dlext=so, d_dlsymun=undef, ccdlflags='-Wl,-E'cccdlflags='-fPIC', lddlflags='-shared -O2 -pipe -march=core2 -fomit-frame-pointer -msse4 -msse4.1 -msse4.2 -mcx16 -msahf -Wl,-O1 -Wl,--as-needed'Characteristics of this binary (from libperl): Compile-time options: MULTIPLICITY PERL_DONT_CREATE_GVSVPERL_IMPLICIT_CONTEXT PERL_MALLOC_WRAP USE_64_BIT_ALLUSE_64_BIT_INT USE_ITHREADS USE_LARGE_FILES USE_PERLIO USE_PERL_ATOF USE_REENTRANT_API Locally applied patches: 0001-gentoo_MakeMaker-RUNPATH.diff 0002-gentoo_config_over.diff 0003-gentoo_cpan_definstalldirs.diff 0004-gentoo_cpanplus_definstalldirs.diff 0005-gentoo_create-libperl-soname.diff 0006-gentoo_MakeMaker-delete_packlist.diff 0007-fixes_8d66b3f9_h2hp_fix.diff 0008-fixes_f178b03b_h2ph_using_deprecated_goto.diff 0009-gentoo_mod-paths.diff 0010-gentoo_enc2xs.diff 0011-gentoo_IO-Compress_AutoLoader_dropped_from_Compress-Zlib.diff 0012-gentoo_drop-fstack-protector.diff Built under linux Compiled at Mar 5 2012 10:55:36 @INC: /etc/perl /usr/lib64/perl5/site_perl/5.12.4/x86_64-linux-thread-multi /usr/lib64/perl5/site_perl/5.12.4 /usr/lib64/perl5/vendor_perl/5.12.4/x86_64-linux-thread-multi /usr/lib64/perl5/vendor_perl/5.12.4 /usr/lib64/perl5/site_perl/5.12.1/x86_64-linux-thread-multi /usr/lib64/perl5/site_perl/5.12.1 /usr/lib64/perl5/site_perl /usr/lib64/perl5/vendor_perl/5.12.3/x86_64-linux-thread-multi /usr/lib64/perl5/vendor_perl/5.12.3 /usr/lib64/perl5/vendor_perl/5.12.2/x86_64-linux-thread-multi /usr/lib64/perl5/vendor_perl/5.12.2 /usr/lib64/perl5/vendor_perl/5.12.1/x86_64-linux-thread-multi /usr/lib64/perl5/vendor_perl/5.12.1 /usr/lib64/perl5/vendor_perl /usr/lib64/perl5/5.12.4/x86_64-linux-thread-multi /usr/lib64/perl5/5.12.4 /usr/local/lib/site_perl On 04/17/2012 04:27 PM, David Mertens wrote:Actually, seconding my second, a good working (failing) example can be the starting point for new failing tests, which are *always* useful in development.On Tue, Apr 17, 2012 at 9:26 AM, David Mertens <[email protected] <mailto:[email protected]>> wrote:On Tue, Apr 17, 2012 at 6:41 AM, Chris Marshall <[email protected] <mailto:[email protected]>> wrote: Hi Ingo, Do you have a simple test case to show the problem? What version of perl+PL... do you have? See BUGS in the PDL distribution for the information needed for problem reports. --Chris Ingo, I will second what Chris said. It's always helpful to give a complete working (i.e. failing) example so that we devs can just copy-and-paste it into an editor and start hacking. Also, the full error message might have been helpful here as well, as it likely would have given us this line number where the dimension mismatch occurred. To make a long story short, this is a bug with PDL::Func. Your best bet is to file a bug report on sourceforge.net <http://sourceforge.net> here: http://sourceforge.net/tracker/?func=browse&group_id=612&atid=100612&status=1 <http://sourceforge.net/tracker/?func=browse&group_id=612&atid=100612&status=1>. I suggest using a title such as "PDL::Func does thread properly" since the code works fine if you want the interpolation at a single x-value. Be sure to read the instructions when adding a bug (i.e. include the output of perldl -V, etc). Even though you did not provide a working example, a difference in dimension handling is usually pretty easy to pick out from the module's code itself, so I decided to simply look over PDL::Func's code. For devs and other interested parties, see PDL::Func (Lib/Func.pm in the repo) starting around line 723. Hermite interpolation uses the chfe routine provided by PDL::Slatec while linear interpolation uses the interpolate routine provided by PDL::Primitive. These two routines expect different kinds of inputs, as evidenced by their signatures (which are far from obvious and which I explain next): chfe (x(n);f(n);d(n);int check();xe(ne);[o]fe(ne);int [o]ierr()) interpolate (xi(); x(n); y(n); [o] yi(); int [o] err()) For the Slatec function chfe(), x(n) and f(n) are the known x and y values to be used in the interpolation. d(n) is... something (the derivatives at these points, perhaps?); the docs do not make it clear. The important point is that the value of x at which we want interpolation are given by xe(ne). For interpolate(), x(n) and y(n) are the known x and y values and the x at which we want interpolation is xi(). NOTICE THE DIFFERENCE IN THE DIMENSION between xi() and xe(ne). This difference was not properly accounted for by Doug and means that the hermite interpolation and linear interpolation interfaces are not interchangeable as written. That sucks, as the whole point of this module was to unify the interface between two distinct packages. Properly reconciling differences in dimension handling is doable but tricky, and writing cogent documentation to explain the dimension handling is even trickier. :-P David On 4/17/12, Ingo Schmid <[email protected] <mailto:[email protected]>> wrote: > Hi, > > I just noticed that when using interpolate, there is an inconsistency > when changing from 'Hermite' to 'Linear' and threading. This is my code. > When I try 'Linear' it barfs with thread dim 0 mismatch. > > my $seq=sequence(eval $p{kx})-$p{kx}/2; > say $r->info,cat($r,$i)->info,$resample->info,$seq->info; > >PDL: Double D [256,128,1,1,2,32]PDL: Double D [256,128,1,1,2,32,2]PDL: > Double D [256]PDL: Double D [256] > my $obj=PDL::Func->init(Interpolate => "Hermite", x > =>$resample, y=>cat($r,$i)); > my $kres=($obj->interpolate($seq))->mv(0,-1); > > > Btw., are there other schemes/modules to try? > Best > Ingo > > > _______________________________________________ Perldl mailing list [email protected] <mailto:[email protected]> http://mailman.jach.hawaii.edu/mailman/listinfo/perldl-- "Debugging is twice as hard as writing the code in the first place.Therefore, if you write the code as cleverly as possible, you are, by definition, not smart enough to debug it." -- Brian Kernighan -- "Debugging is twice as hard as writing the code in the first place. Therefore, if you write the code as cleverly as possible, you are, by definition, not smart enough to debug it." -- Brian Kernighan_______________________________________________ Perldl mailing list [email protected] http://mailman.jach.hawaii.edu/mailman/listinfo/perldl
test.pl
Description: Perl program
_______________________________________________ Perldl mailing list [email protected] http://mailman.jach.hawaii.edu/mailman/listinfo/perldl
