From bug-request at octave dot org Sat Jan 15 15:14:42 2005
Subject: feval doesn't return all values/wrong values
From: Markus Jordan <mj604 at soton dot ac dot uk>
To: bug-octave at bevo dot che dot wisc dot edu
CC: mj604 at soton dot ac dot uk
Date: Sat, 15 Jan 2005 13:49:45 -0600

Bug report for Octave 2.1.49 configured for i686-pc-linux-gnu

Description:
-----------
 
  The Problem occurs while using the Netlab-software for matlab (  
http://www.ncrg.aston.ac.uk/netlab/index.php. In order to use Netlab the 
package octave-forge is also necessary.). Instead of returning the three 
variables, calculated in a function called via feval, it returns only 
the first one but three times. This is the according feval-call:
   [s{1:nargout}] = feval(errstr, net, x, t)
   The called function returns one value, and two matrices, but if you 
just print 's' onto the screen it only contains the first value three times.

The error occurs not only on this system, but also on a gentoo-based 
system with octave 2.1.50/57/64 and octave-forge/20040707.

Repeat-By:
---------

   %begin script for reproducing bug%
  mlp_type = 'linear';
  ALG='scg';
  nhidden = 10;
  ndata = 300;
  nin = 10;
  nout = 1;
  aw1 = 0.0001*ones(1,nin);
  ab1 = 0.0001;
  aw2 = 0.0001;
  ab2 = 0.0001;
  beta=200;
  [x y]=friedman_generate(ndata);
  prior = mlpprior(nin, nhidden, nout, aw1, ab1, aw2, ab2);
  net = mlp(nin, nhidden, nout, mlp_type, prior, beta);
  net = mlpinit(net, prior);

options = zeros(1,18);

for i = 1:1
   [net, aoptions] = netopt( net, options, x, y, ALG);
   [net, gamma ] = evidence(net, x, y, 30);

end
%end script for reproducing bug%

necessary function to generate the data  
function [X, y] = friedman_generate(ndata)
  X = randn(ndata,10);
  for i= 1:ndata
     y_wo(i) = 0.1*exp(4*X(i,1))+4/( 1 + exp( -20*( X(i,2)-1/2))) + 
3*X(i,3) + 2*X(i,4) + X(i,5);
     y(i)=y_wo(i)+ randn(1);
  end
  y = y';
%end of function friedman_generate

If the script is executed it produces an error that the index (2) is not 
available in file evidence.m line 88.
It seems that the datastructure eprior is a single value.
The value eprior is produced in line 55 from the function neterr.
This function neterr is calling the function mlperr via the feval-function.
The function mlperr generates the three return values in the last line.
If you print these values, you see the first one is a single value, and 
the 2nd and 3rd are matrices.
Back again in the neterr function with the feval-call it reassigns the 
returned values of the mlperr ( via feval ) for returning them to the 
calling function. But in this stage the returned value of the feval 
function contains only three times the first value of the three 
returning values of the mlperr-function.

This piece of code works in Matlab ( Without an error)

Fix:
---
 
 

Configuration (please do not edit this section):
-----------------------------------------------

uname output:     Linux corwin.ecs.soton.ac.uk 2.4.21-20.ELsmp #1 SMP 
Wed Aug 18 20:46:40 EDT 2004 i686 i686 i386 GNU/Linux
configure opts:   '--prefix=/usr/cslib/utils/octave-2.1.49'
Fortran compiler: g77
FFLAGS:           -O
F2C:              
F2CFLAGS:         
FLIBS:            -L/usr/lib/gcc-lib/i386-redhat-linux/2.96 
-L/usr/lib/gcc-lib/i386-redhat-linux/2.96/../../.. -lg2c -lm
CPPFLAGS:         
INCFLAGS:         -I. -I. -I./liboctave -I./src -I./libcruft/misc  
-I./glob -I./glob
C compiler:       gcc, version 2.96 20000731 (Red Hat Linux 7.3 2.96-113)
CFLAGS:           -g -O2
CPICFLAG:         -fPIC
C++ compiler:     g++, version 2.96
CXXFLAGS:         -g -O2 -Wall
CXXPICFLAG:       -fPIC
LD_CXX:           g++
LDFLAGS:          -g
LIBFLAGS:         -L.
RLD_FLAG:         
BLAS_LIBS:        -llapack -lblas
FFTW_LIBS:        
LIBS:             -lreadline  -lncurses -lm
LEXLIB:           
LIBPLPLOT:        
LIBDLFCN:         
LIBGLOB:          ./glob/glob.o ./glob/fnmatch.o
SED:              /bin/sed
DEFS:

  -DPACKAGE_NAME="" -DPACKAGE_TARNAME="" -DPACKAGE_VERSION=""
  -DPACKAGE_STRING="" -DPACKAGE_BUGREPORT="" -DOCTAVE_SOURCE=1
  -DSTDC_HEADERS=1 -DHAVE_SYS_TYPES_H=1 -DHAVE_SYS_STAT_H=1
  -DHAVE_STDLIB_H=1 -DHAVE_STRING_H=1 -DHAVE_MEMORY_H=1 -DHAVE_STRINGS_H=1
  -DHAVE_INTTYPES_H=1 -DHAVE_STDINT_H=1 -DHAVE_UNISTD_H=1 -DSEPCHAR=1
  -DSEPCHAR_STR=":" -D__NO_MATH_INLINES=1 -DCXX_NEW_FRIEND_TEMPLATE_DECL=1
  -DCXX_ABI=gnu_v2 -DHAVE_LIBM=1 -DF77_FUNC(name,NAME)=name ## _
  -DF77_FUNC_(name,NAME)=name ## __ -DHAVE_BLAS=1 -DHAVE_GETHOSTNAME=1
  -DHAVE_GETPWNAM=1 -DHAVE_DEV_T=1 -DHAVE_INO_T=1 -DHAVE_NLINK_T=1
  -DHAVE_NLINK_T=1 -DHAVE_LONG_LONG_INT=1 -DHAVE_UNSIGNED_LONG_LONG_INT=1
  -DHAVE_SIGSET_T=1 -DHAVE_SIG_ATOMIC_T=1 -DSIZEOF_SHORT=2 -DSIZEOF_INT=4
  -DSIZEOF_LONG=4 -DSIZEOF_LONG_LONG=8 -DHAVE_ALLOCA_H=1 -DHAVE_ALLOCA=1
  -DNPOS=std::string::npos -DHAVE_DYNAMIC_AUTO_ARRAYS=1 -DSTDC_HEADERS=1
  -DHAVE_DIRENT_H=1 -DTIME_WITH_SYS_TIME=1 -DHAVE_SYS_WAIT_H=1
  -DHAVE_ASSERT_H=1 -DHAVE_CURSES_H=1 -DHAVE_DLFCN_H=1 -DHAVE_FCNTL_H=1
  -DHAVE_FLOAT_H=1 -DHAVE_GRP_H=1 -DHAVE_LIMITS_H=1 -DHAVE_MEMORY_H=1
  -DHAVE_NCURSES_H=1 -DHAVE_POLL_H=1 -DHAVE_PWD_H=1 -DHAVE_STDLIB_H=1
  -DHAVE_STRING_H=1 -DHAVE_SYS_IOCTL_H=1 -DHAVE_SYS_PARAM_H=1
  -DHAVE_SYS_POLL_H=1 -DHAVE_SYS_RESOURCE_H=1 -DHAVE_SYS_SELECT_H=1
  -DHAVE_SYS_STAT_H=1 -DHAVE_SYS_TIME_H=1 -DHAVE_SYS_TIMES_H=1
  -DHAVE_SYS_TYPES_H=1 -DHAVE_SYS_UTSNAME_H=1 -DHAVE_TERMCAP_H=1
  -DHAVE_UNISTD_H=1 -DHAVE_VARARGS_H=1 -DHAVE_SSTREAM=1 -DHAVE_TERMIO_H=1
  -DHAVE_SGTTY_H=1 -DHAVE_GLOB_H=1 -DHAVE_FNMATCH_H=1 -DHAVE_ATEXIT=1
  -DHAVE_BASENAME=1 -DHAVE_BCOPY=1 -DHAVE_BZERO=1 -DHAVE_DUP2=1
  -DHAVE_ENDGRENT=1 -DHAVE_ENDPWENT=1 -DHAVE_EXECVP=1 -DHAVE_FCNTL=1
  -DHAVE_FORK=1 -DHAVE_GETCWD=1 -DHAVE_GETEGID=1 -DHAVE_GETEUID=1
  -DHAVE_GETGID=1 -DHAVE_GETGRENT=1 -DHAVE_GETGRGID=1 -DHAVE_GETGRNAM=1
  -DHAVE_GETPGRP=1 -DHAVE_GETPID=1 -DHAVE_GETPPID=1 -DHAVE_GETPWENT=1
  -DHAVE_GETPWUID=1 -DHAVE_GETTIMEOFDAY=1 -DHAVE_GETUID=1 -DHAVE_GETWD=1
  -DHAVE_KILL=1 -DHAVE_LINK=1 -DHAVE_LOCALTIME_R=1 -DHAVE_LSTAT=1
  -DHAVE_MEMMOVE=1 -DHAVE_MKDIR=1 -DHAVE_MKFIFO=1 -DHAVE_MKSTEMP=1
  -DHAVE_ON_EXIT=1 -DHAVE_PIPE=1 -DHAVE_POLL=1 -DHAVE_PUTENV=1
  -DHAVE_RAISE=1 -DHAVE_READLINK=1 -DHAVE_RENAME=1 -DHAVE_RINDEX=1
  -DHAVE_RMDIR=1 -DHAVE_SELECT=1 -DHAVE_SETGRENT=1 -DHAVE_SETPWENT=1
  -DHAVE_SETVBUF=1 -DHAVE_SIGACTION=1 -DHAVE_SIGLONGJMP=1
  -DHAVE_SIGPENDING=1 -DHAVE_SIGPROCMASK=1 -DHAVE_SIGSUSPEND=1
  -DHAVE_STAT=1 -DHAVE_STRCASECMP=1 -DHAVE_STRDUP=1 -DHAVE_STRERROR=1
  -DHAVE_STRFTIME=1 -DHAVE_STRNCASECMP=1 -DHAVE_STRPTIME=1 -DHAVE_SYMLINK=1
  -DHAVE_TEMPNAM=1 -DHAVE_UMASK=1 -DHAVE_UNLINK=1 -DHAVE_USLEEP=1
  -DHAVE_VFPRINTF=1 -DHAVE_VSPRINTF=1 -DHAVE_VSNPRINTF=1 -DHAVE_WAITPID=1
  -DHAVE_TIMEVAL=1 -DHAVE_FINITE=1 -DHAVE_ISNAN=1 -DHAVE_ISINF=1
  -DHAVE_COPYSIGN=1 -DHAVE_ACOSH=1 -DHAVE_ASINH=1 -DHAVE_ATANH=1
  -DHAVE_ERF=1 -DHAVE_ERFC=1 -DHAVE_STRUCT_STAT_ST_BLKSIZE=1
  -DHAVE_STRUCT_STAT_ST_BLOCKS=1 -DHAVE_STRUCT_STAT_ST_RDEV=1
  -DHAVE_STRUCT_TM_TM_ZONE=1 -DHAVE_TM_ZONE=1 -DUSE_READLINE=1
  -DEXCEPTION_IN_MATH=1 -DRETSIGTYPE=void -DSYS_SIGLIST_DECLARED=1
  -DHAVE_SYS_SIGLIST=1 -DHAVE_POSIX_SIGNALS=1 -DHAVE_GETRUSAGE=1
  -DHAVE_TIMES=1 -DYYTEXT_POINTER=1 -DGNUPLOT_BINARY="gnuplot"
  -DGNUPLOT_HAS_MULTIPLOT=1 -DGNUPLOT_HAS_FRAMES=1

User-preferences (please do not edit this section):
--------------------------------------------------

  EDITOR = "emacs"
  EXEC_PATH = 
":/usr/software/bin:/bin:/usr/cslib/bin:/usr/local/cslib:.:/usr/kerberos/bin:/usr/software/bin:/bin:/usr/bin:/usr/local/bin:/usr/X11R6/bin:/usr/X11R6/bin:/usr/openwin/bin:/usr/X11R6/bin:/usr/X11R6/bin:/usr/openwin/bin:/usr/X11R6/bin"
  IMAGEPATH = 
".:/usr/cslib/utils/octave-2.1.49/share/octave/2.1.49/imagelib//"
  INFO_FILE = "/usr/cslib/utils/octave-2.1.49/info/octave.info"
  INFO_PROGRAM = "info"
  LOADPATH = ":"
  PAGER = "less -e -P'-- less ?pB(%pB\\%):--. (f)orward, (b)ack, (q)uit$'"
  PS1 = "\\s:\\#> "
  PS2 = "> "
  PS4 = "+ "
  automatic_replot = 0
  beep_on_error = 0
  completion_append_char = " "
  default_eval_print_flag = 1
# default_global_variable_value = <no value or error in displaying it>
  default_return_value = []
  default_save_format = "ascii"
  define_all_return_values = 0
  do_fortran_indexing = 0
  echo_executing_commands = 0
  empty_list_elements_ok = 1
  fixed_point_format = 0
  gnuplot_binary = "gnuplot"
  gnuplot_command_end = "\n"
  gnuplot_command_plot = "pl"
  gnuplot_command_replot = "rep"
  gnuplot_command_splot = "sp"
  gnuplot_command_title = "t"
  gnuplot_command_using = "u"
  gnuplot_command_with = "w"
  gnuplot_has_frames = 1
  gnuplot_has_multiplot = 1
  history_file = "/home/mj604/.octave_hist"
  history_size = 1024
  ignore_function_time_stamp = "system"
  implicit_num_to_str_ok = 0
  implicit_str_to_num_ok = 0
  initialize_global_variables = 0
  max_recursion_depth = 256
  ok_to_lose_imaginary_part = "warn"
  output_max_field_width = 10
  output_precision = 5
  page_output_immediately = 0
  page_screen_output = 1
  prefer_column_vectors = 1
  print_answer_id_name = 1
  print_empty_dimensions = 1
  print_rhs_assign_val = 0
  propagate_empty_matrices = 1
  resize_on_range_error = 1
  return_last_computed_value = 0
  save_precision = 15
  saving_history = 1
  silent_functions = 0
  split_long_rows = 1
  string_fill_char = " "
  struct_levels_to_print = 2
  suppress_verbose_help_message = 0
  treat_neg_dim_as_zero = 0
  warn_assign_as_truth_value = 1
  warn_divide_by_zero = 1
  warn_function_name_clash = 1
  warn_future_time_stamp = 1
  warn_missing_semicolon = 0
  warn_variable_switch_label = 0
  whitespace_in_literal_matrix =



-------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.

Octave's home on the web:  http://www.octave.org
How to fund new projects:  http://www.octave.org/funding.html
Subscription information:  http://www.octave.org/archive.html
-------------------------------------------------------------