From bug-octave-request at bevo dot che dot wisc dot edu  Mon Feb  2 08:57:40 1998
Subject: Unidentified subject!
From: Christian Buchta <buchta at eternity dot wu-wien dot ac dot at>
To: bug-octave at bevo dot che dot wisc dot edu
Date: Mon, 2 Feb 1998 15:55:55 +0100

To: bug-octave at bevo dot che dot wisc dot edu
Subject: problem with built-in variable e

Bug report for Octave 2.0.9 configured for i386-suse-linux-gnu

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

  First of all, thanks to the people who have developped this
  software and made it available for everyone, it's really
  useful for my scientific work.

  I think the best way to see what the problem is, is that I
  copy the error-producing script and m-file into the next 
  section. The first should have name bug_script.m and the 
  second am_corr.m

  Since I am currently connected via ppp and I don' know
  what e-mail header this automatic bug-report will append,
  my adress is christian dot buchta at wu-wien dot ac dot at - there is no
  eternity.wu-wien.ac.at !!!

  Thanks

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

%
% bug_script
%
% e is a built-in variable that I can 'overload' with something
% I need to stand e for? It seems that top level 'overloading'
% does not survive function level 'overloading'.
%
% In statistics, for example, some vector of expectations is
% easily memorized by having name 'e'. So if this is not a 
% bug it is at least no feature  ;-)
%  

echo on

p=rand(5,5);
e=rand(5,1)

[p,v]=am_corr(p,e)

% show that e has changed to the value given in am_corr

e

%
% end of bug-script


% am_corr - calculate the correlations of an autologistic model
%
% [C,V]=am_corr(S,E)
%
% [C,V]=am_corr(C,E,V,A,D)
%
% C  on input an at least lower triangular without diagonal
%    matrix of marginal sums; usually PS from a generator like am_dfmg 
%    or am_rmg, or from am_sum. On output the sums are replaced
%    by correlations
% E  expectation vector; usally from the functions as mentioned above
% A  row/column index where a contiguous submodel block starts
% D  dimension of the submodel. The last row/column index of
%    the submodel is A+D-1
% V  vector of variances of the model variables
% TV the total variance of the autologistic model
%
% (C) c.buchta Oct. 1997 Vers.0.0
%     send bug-reports to Christian dot Buchta at wu-wien dot ac dot at
%

function [C,V,TV]=am_corr(C,E,V,A,D)
[mc,nc]=size(C);

if mc~=nc
   error('C must be square');
end

if nargin==2
   A=1;
   D=mc;
   e=mc;
else
   e=A+D-1;
end

if e>mc
   error('dimension of submodel exceeds dimension of P');
end

[me,ne]=size(E);

if ne~=1 
   error('E must be a column vector');
end

if me~=mc
   error('dimension of E must match that of P');
end

V(A:e,1)=(1-E(A:e)).*E(A:e);
TV=sum(V(A:e));


for i=A:e-1
   for j=i+1:e
      C(j,i)=(C(j,i)-E(i)*E(j))/sqrt(V(i)*V(j));
   end
end

%
% end of am_corr

Fix:
---

  * If possible, replace this item with a description of how to
    fix the problem (if you don't have a fix for the problem, don't
    include this section, but please do submit your report anyway).

  Sorry, no idea.


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

uname output:     Linux eternity 2.0.32 #12 Tue Nov 18 11:49:22 MET 1997 i686 unknown
configure opts:   --prefix=/usr --exec-prefix=/usr --enable-shared --enable-lite-kernel --host=i386-suse-linux
Fortran compiler: g77
FFLAGS:           -O
F2C:              
F2CFLAGS:         
FLIBS:            -lf2c -lm -L/usr/lib/gcc-lib/i486-linux/2.7.2.1 -L/usr/i486-linux/lib -lm -lgcc -lc
CPPFLAGS:         
INCFLAGS:         -I/usr/include -I/usr/include/octave-2.0.9
C compiler:       gcc, version 2.7.2.1
CFLAGS:           -DHAVE_CONFIG_H -mieee-fp -g -O2 -Wall
CPICFLAG:         -fPIC
C++ compiler:     c++, version 2.7.2.1
CXXFLAGS:         -DHAVE_CONFIG_H  -fno-implicit-templates -mieee-fp -g -O2 -Wall
CXXPICFLAG:       -fPIC
LDFLAGS:          -g
LIBFLAGS:         -L/usr/lib
RLD_FLAG:         -Xlinker -rpath -Xlinker /usr/lib
CXXLIBS:          -lstdc++ -lm -L/usr/lib/gcc-lib/i486-linux/2.7.2.1 -L/usr/i486-linux/lib -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 -DOCTAVE_LITE=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_LIBDL=1 -DHAVE_DLOPEN=1
  -DHAVE_DLSYM=1 -DHAVE_DLERROR=1 -DHAVE_DLCLOSE=1 -DWITH_DL=1
  -DWITH_DYNAMIC_LINKING=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

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

  EDITOR = "emacs"
  IMAGEPATH = ".:/usr/share/octave/2.0.9/imagelib//"
  INFO_FILE = "/usr/info/octave.info"
  LOADPATH = ".:/usr/libexec/octave/site/oct/i386-suse-linux-gnu//:/usr/share/octave/site/m//:/usr/libexec/octave/2.0.9/oct/i386-suse-linux-gnu//:/usr/share/octave/2.0.9/m//"
  OCTAVE_VERSION = "2.0.9"
  PAGER = "less"
  PS1 = "\\s:\\#> "
  PS2 = "> "
  automatic_replot = 0
  whitespace_in_literal_matrix = ""
  default_save_format = "ascii"
  do_fortran_indexing = 0
  empty_list_elements_ok = "warn"
  eps = 2.22045e-16
  gnuplot_binary = "gnuplot"
  ignore_function_time_stamp = "system"
  implicit_str_to_num_ok = 0
  ok_to_lose_imaginary_part = "warn"
  output_max_field_width = 10
  output_precision = 5
  page_screen_output = 1
  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