From bug-octave-request at bevo dot che dot wisc dot edu  Mon Dec 29 06:25:09 1997
Subject: incorrect computation of interpolated values in polyfit.m
From: Mario Storti <mstorti at minerva dot unl dot edu dot ar>
To: bug-octave at bevo dot che dot wisc dot edu
cc: mstorti at minerva dot unl dot edu dot ar
Date: Mon, 29 Dec 1997 09:26:35 -0300

Bug report for Octave 2.0.9 configured for i586-pc-linux-gnu

Description:
-----------

  * In polyfit.m, when computing the interpolated values `yf', it
    seems that the coefficients of the polynomial should be reversed
    (flipud) again before the matrix multiplication. 


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

  * 

Fix:
---

  * Replace (in polyfit.m):

>  if (nargout == 2)
>    yf = X * p;
>  endif

    by:

>  if (nargout == 2)
>    yf = X * flipud(p);
>  endif



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

uname output:     Linux minerva 2.0.18 #1 Tue Sep 10 10:15:48 EDT 1996 i586
configure opts:   --disable-shared --disable-dl --disable-shl
Fortran compiler: g77
FFLAGS:           -O
F2C:              
F2CFLAGS:         
FLIBS:            -lf2c -lm -L/usr/lib/gcc-lib/i486-linux/2.7.2 -lm -lgcc -lc
CPPFLAGS:         
INCFLAGS:         -I/usr/local/include -I/usr/local/include/octave-2.0.9
C compiler:       gcc, version 2.7.2
CFLAGS:           -DHAVE_CONFIG_H -mieee-fp -O
CPICFLAG:         -fPIC
C++ compiler:     c++, version 2.7.2
CXXFLAGS:         -DHAVE_CONFIG_H  -fno-implicit-templates -mieee-fp -O
CXXPICFLAG:       -fPIC
LDFLAGS:          -static
LIBFLAGS:         -L/usr/local/lib
RLD_FLAG:         -Xlinker -rpath -Xlinker /usr/local/lib
CXXLIBS:          -lstdc++ -lm -L/usr/lib/gcc-lib/i486-linux/2.7.2 -lstdc++ -lm -lgcc -lc -lgcc
TERMLIBS:         -lncurses
LIBS:             
LEXLIB:           
LIBPLPLOT:        
LIBDLFCN:         
DEFS:

  -DOCTAVE_SOURCE=1 -DSEPCHAR=':' -DSEPCHAR_STR=":" -DUSE_GNU_INFO=1
  -DUSE_READLINE=1 -DF77_APPEND_UNDERSCORE=1 -DSIZEOF_SHORT=2 -DSIZEOF_INT=4
  -DSIZEOF_LONG=4 -DHAVE_ALLOCA_H=1 -DHAVE_ALLOCA=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_FCNTL_H=1 -DHAVE_FLOAT_H=1
  -DHAVE_GRP_H=1 -DHAVE_LIMITS_H=1 -DHAVE_MEMORY_H=1 -DHAVE_NAN_H=1
  -DHAVE_NCURSES_H=1 -DHAVE_PWD_H=1 -DHAVE_STDLIB_H=1 -DHAVE_STRING_H=1
  -DHAVE_SYS_PARAM_H=1 -DHAVE_SYS_RESOURCE_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_TERMIO_H=1
  -DHAVE_TERMIOS_H=1 -DHAVE_UNISTD_H=1 -DHAVE_VARARGS_H=1
  -DHAVE_ATEXIT=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_GETHOSTNAME=1 -DHAVE_GETPGRP=1 -DHAVE_GETPID=1 -DHAVE_GETPPID=1
  -DHAVE_GETPWENT=1 -DHAVE_GETPWNAM=1 -DHAVE_GETPWUID=1 -DHAVE_GETUID=1
  -DHAVE_GETWD=1 -DHAVE_LSTAT=1 -DHAVE_MEMMOVE=1 -DHAVE_MKDIR=1
  -DHAVE_MKFIFO=1 -DHAVE_ON_EXIT=1 -DHAVE_PIPE=1 -DHAVE_PUTENV=1
  -DHAVE_RENAME=1 -DHAVE_RINDEX=1 -DHAVE_RMDIR=1 -DHAVE_SETGRENT=1
  -DHAVE_SETPWENT=1 -DHAVE_SETVBUF=1 -DHAVE_SIGACTION=1 -DHAVE_SIGPENDING=1
  -DHAVE_SIGPROCMASK=1 -DHAVE_SIGSUSPEND=1 -DHAVE_STAT=1 -DHAVE_STRCASECMP=1
  -DHAVE_STRDUP=1 -DHAVE_STRERROR=1 -DHAVE_STRNCASECMP=1 -DHAVE_TEMPNAM=1
  -DHAVE_UMASK=1 -DHAVE_UNLINK=1 -DHAVE_USLEEP=1 -DHAVE_VFPRINTF=1
  -DHAVE_VSPRINTF=1 -DHAVE_WAITPID=1 -DHAVE_PROGRAM_INVOCATION_NAME=1
  -DHAVE_LIBM=1 -DHAVE_FINITE=1 -DHAVE_ISNAN=1 -DHAVE_ISINF=1
  -DHAVE_ACOSH=1 -DHAVE_ASINH=1 -DHAVE_ATANH=1 -DHAVE_ERF=1 -DHAVE_ERFC=1
  -DHAVE_ST_BLKSIZE=1 -DHAVE_ST_BLOCKS=1 -DHAVE_ST_RDEV=1 -DHAVE_TZNAME=1
  -DHAVE_GR_PASSWD=1 -DRETSIGTYPE=void -DSYS_SIGLIST_DECLARED=1
  -DHAVE_SYS_SIGLIST=1 -DHAVE_POSIX_SIGNALS=1 -DHAVE_GETRUSAGE=1
  -DHAVE_TIMES=1 -DGNUPLOT_HAS_MULTIPLOT=1 -DGNUPLOT_HAS_FRAMES=1

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

  EDITOR = "emacsclient"
  IMAGEPATH = ".:/usr/local/share/octave/2.0.9/imagelib//"
  INFO_FILE = "/usr/local/info/octave.info"
  LOADPATH = ".:~/OCTAVE/TOOLS//:~/SHIP/INTEF/RUN/TOOLS/:/usr/local/share/octave/2.0.9/m//:/usr/local/share/octave/site//"
  OCTAVE_VERSION = "2.0.9"
  PAGER = "less -e -P'-- less ?pB(%pB\\%):--. (f)orward, (b)ack, (q)uit$'"
  PS1 = "\\s> "
  PS2 = "> "
  automatic_replot = 0
  whitespace_in_literal_matrix = ""
  default_save_format = "ascii"
  do_fortran_indexing = 0
  empty_list_elements_ok = 1
  eps = 2.22045e-16
  gnuplot_binary = "gnuplot -geometry 300x400+1+1"
  ignore_function_time_stamp = "system"
  implicit_str_to_num_ok = 1
  ok_to_lose_imaginary_part = "warn"
  output_max_field_width = 10
  output_precision = 5
  page_screen_output = 0
  prefer_column_vectors = 1
  prefer_zero_one_indexing = 0
  print_answer_id_name = 1
  print_empty_dimensions = 1
  propagate_empty_matrices = 1
  resize_on_range_error = 1
  return_last_computed_value = 0
  save_precision = 15
  silent_functions = 0
  split_long_rows = 1
  treat_neg_dim_as_zero = 0
  warn_assign_as_truth_value = 1
  warn_comma_in_global_decl = 1
  warn_divide_by_zero = 1