--- gforth/prim 2000/07/14 08:55:15 1.47 +++ gforth/prim 2000/08/11 19:49:39 1.53 @@ -36,8 +36,9 @@ \ be separated by at least one empty line \ \ Both pronounciation and stack items (in the stack effect) must -\ conform to the C name syntax or the C compiler will complain. -\ +\ conform to the C identifier syntax or the C compiler will complain. +\ If you don't have a pronounciation field, the Forth name is used, +\ and has to conform to the C identifier syntax. \ \ These specifications are automatically translated into C-code for the \ interpreter and into some other files. I hope that your C compiler has @@ -444,8 +445,9 @@ n = rp[4]; \ digit is high-level: 0/0% move ( c_from c_to ucount -- ) core -""Copy the contents of @i{ucount} address units at @i{c-from} to +""Copy the contents of @i{ucount} aus at @i{c-from} to @i{c-to}. @code{move} works correctly even if the two areas overlap."" +/* !! note that the standard specifies addr, not c-addr */ memmove(c_to,c_from,ucount); /* make an Ifdef for bsd and others? */ : @@ -474,8 +476,7 @@ while (u-- > 0) DO 1- dup c@ I c! -1 +LOOP drop ; fill ( c_addr u c -- ) core -"" If @i{u}>0, store character @i{c} in each of @i{u} consecutive -@code{char} addresses in memory, starting at address @i{c-addr}."" +""Store @i{c} in @i{u} chars starting at @i{c-addr}."" memset(c_addr,c,u); : -rot bounds @@ -510,7 +511,7 @@ else if (n>0) : swap bounds ?DO dup c@ I c@ = WHILE 1+ LOOP drop 0 - ELSE c@ I c@ - unloop THEN -text-flag ; + ELSE c@ I c@ - unloop THEN sgn ; : sgn ( n -- -1/0/1 ) dup 0= IF EXIT THEN 0< 2* 1+ ; @@ -532,7 +533,7 @@ else if (n>0) ?DO dup c@ I c@ <> IF dup c@ toupper I c@ toupper = ELSE true THEN WHILE 1+ LOOP drop 0 - ELSE c@ toupper I c@ toupper - unloop THEN -text-flag ; + ELSE c@ toupper I c@ toupper - unloop THEN sgn ; -trailing ( c_addr u1 -- c_addr u2 ) string dash_trailing ""Adjust the string specified by @i{c-addr, u1} to remove all trailing @@ -599,11 +600,11 @@ else : 2dup > IF swap THEN drop ; -abs ( n1 -- n2 ) core -if (n1<0) - n2 = -n1; +abs ( n -- u ) core +if (n<0) + u = -n; else - n2 = n1; + u = n; : dup 0< IF negate THEN ; @@ -629,12 +630,14 @@ n3 = n1%n2; /* !! is this correct? look >r s>d r> fm/mod ; 2* ( n1 -- n2 ) core two_star +""Shift left by 1; also works on unsigned numbers"" n2 = 2*n1; : dup + ; 2/ ( n1 -- n2 ) core two_slash -/* !! is this still correct? */ +""Arithmetic shift right by 1. For signed numbers this is a floored +division by 2 (note that @code{/} not necessarily floors)."" n2 = n1>>1; : dup MINI and IF 1 ELSE 0 THEN @@ -776,6 +779,7 @@ d2 = -d1; invert swap negate tuck 0= - ; d2* ( d1 -- d2 ) double d_two_star +""Shift left by 1; also works on unsigned numbers"" #ifdef BUGGY_LONG_LONG d2.lo = d1.lo<<1; d2.hi = (d1.hi<<1) | (d1.lo>>(CELL_BITS-1)); @@ -786,6 +790,8 @@ d2 = 2*d1; 2dup d+ ; d2/ ( d1 -- d2 ) double d_two_slash +""Arithmetic shift right by 1. For signed numbers this is a floored +division by 2."" #ifdef BUGGY_LONG_LONG d2.hi = d1.hi>>1; d2.lo= (d1.lo>>1) | (d1.hi<<(CELL_BITS-1)); @@ -813,6 +819,7 @@ w2 = ~w1; MAXU xor ; rshift ( u1 n -- u2 ) core r_shift +""Logical shift right by @i{n} bits."" u2 = u1>>n; : 0 ?DO 2/ MAXI and LOOP ; @@ -938,7 +945,7 @@ about this word is to consider the numbe around from @code{max-u} to 0 for unsigned, and from @code{max-n} to min-n for signed numbers); now consider the range from u2 towards increasing numbers up to and excluding u3 (giving an empty range if -u2=u3; if u1 is in this range, @code{within} returns true."" +u2=u3); if u1 is in this range, @code{within} returns true."" f = FLAG(u1-u2 < u3-u2); : over - >r - r> u< ; @@ -971,41 +978,48 @@ fp = f_addr; SET_IP((Xt *)(*rp++)); >r ( w -- ) core to_r +""@code{( R: -- w )}"" *--rp = w; : (>r) ; : (>r) rp@ cell+ @ rp@ ! rp@ cell+ ! ; r> ( -- w ) core r_from +""@code{( R: w -- )}"" w = *rp++; : rp@ cell+ @ rp@ @ rp@ cell+ ! (rdrop) rp@ ! ; Create (rdrop) ' ;s A, rdrop ( -- ) gforth +""@code{( R: w -- )}"" rp++; : r> r> drop >r ; 2>r ( w1 w2 -- ) core-ext two_to_r +""@code{( R: -- w1 w2 )}"" *--rp = w1; *--rp = w2; : swap r> swap >r swap >r >r ; 2r> ( -- w1 w2 ) core-ext two_r_from +""@code{( R: w1 w2 -- )}"" w2 = *rp++; w1 = *rp++; : r> r> swap r> swap >r swap ; 2r@ ( -- w1 w2 ) core-ext two_r_fetch +""@code{( R: w1 w2 -- w1 w2 )}"" w2 = rp[0]; w1 = rp[1]; : i' j ; 2rdrop ( -- ) gforth two_r_drop +""@code{( R: w1 w2 -- )}"" rp+=2; : r> r> drop r> drop >r ; @@ -1049,6 +1063,8 @@ tuck ( w1 w2 -- w2 w1 w2 ) core-ext swap over ; ?dup ( w -- w ) core question_dupe +""Actually the stack effect is: @code{( w -- 0 | w w )}. It performs a +@code{dup} if w is nonzero."" if (w!=0) { IF_TOS(*sp-- = w;) #ifndef USE_TOS @@ -1059,6 +1075,7 @@ if (w!=0) { dup IF dup THEN ; pick ( u -- w ) core-ext +""Actually the stack effect is @code{ x0 ... xu u -- x0 ... xu x0 }."" w = sp[u+1]; : 1+ cells sp@ + @ ; @@ -1094,21 +1111,21 @@ w = sp[u+1]; \ toggle is high-level: 0.11/0.42% @ ( a_addr -- w ) core fetch -"" Read from the cell at address @i{a-addr}, and return its contents, @i{w}."" +""@i{w} is the cell stored at @i{a_addr}."" w = *a_addr; ! ( w a_addr -- ) core store -"" Write the value @i{w} to the cell at address @i{a-addr}."" +""Store @i{w} into the cell at @i{a-addr}."" *a_addr = w; +! ( n a_addr -- ) core plus_store -"" Add @i{n} to the value stored in the cell at address @i{a-addr}."" +""Add @i{n} to the cell at @i{a-addr}."" *a_addr += n; : tuck @ + swap ! ; c@ ( c_addr -- c ) core c_fetch -"" Read from the char at address @i{c-addr}, and return its contents, @i{c}."" +""@i{c} is the char stored at @i{c_addr}."" c = *c_addr; : [ bigendian [IF] ] @@ -1134,7 +1151,7 @@ c = *c_addr; : 8>> 2/ 2/ 2/ 2/ 2/ 2/ 2/ 2/ ; c! ( c c_addr -- ) core c_store -"" Write the value @i{c} to the char at address @i{c-addr}."" +""Store @i{c} into the char at @i{c-addr}."" *c_addr = c; : [ bigendian [IF] ] @@ -1164,28 +1181,28 @@ c! ( c c_addr -- ) core c_store : 8<< 2* 2* 2* 2* 2* 2* 2* 2* ; 2! ( w1 w2 a_addr -- ) core two_store -"" Write the value @i{w1, w2} to the double at address @i{a-addr}."" +""Store @i{w2} into the cell at @i{c-addr} and @i{w1} into the next cell."" a_addr[0] = w2; a_addr[1] = w1; : tuck ! cell+ ! ; 2@ ( a_addr -- w1 w2 ) core two_fetch -"" Read from the double at address @i{a-addr}, and return its contents, @i{w1, w2}."" +""@i{w2} is the content of the cell stored at @i{a-addr}, @i{w1} is +the content of the next cell."" w2 = a_addr[0]; w1 = a_addr[1]; : dup cell+ @ swap @ ; cell+ ( a_addr1 -- a_addr2 ) core cell_plus -"" Increment @i{a-addr1} by the number of address units corresponding to the size of -one cell, to give @i{a-addr2}."" +""@code{1 cells +}"" a_addr2 = a_addr1+1; : cell + ; cells ( n1 -- n2 ) core -"" @i{n2} is the number of address units corresponding to @i{n1} cells."" +"" @i{n2} is the number of address units of @i{n1} cells."" n2 = n1 * sizeof(Cell); : [ cell @@ -1196,8 +1213,7 @@ n2 = n1 * sizeof(Cell); drop ] ; char+ ( c_addr1 -- c_addr2 ) core char_plus -"" Increment @i{c-addr1} by the number of address units corresponding to the size of -one char, to give @i{c-addr2}."" +""@code{1 chars +}."" c_addr2 = c_addr1 + 1; : 1+ ; @@ -1503,6 +1519,8 @@ struct timeval time1; struct timezone zone1; struct tm *ltime; gettimeofday(&time1,&zone1); +/* !! Single Unix specification: + If tzp is not a null pointer, the behaviour is unspecified. */ ltime=localtime((time_t *)&time1.tv_sec); nyear =ltime->tm_year+1900; nmonth=ltime->tm_mon+1; @@ -1523,16 +1541,15 @@ allocate ( u -- a_addr wior ) memory contents of the data space is undefined. If the allocation is successful, @i{a-addr} is the start address of the allocated region and @i{wior} is 0. If the allocation fails, @i{a-addr} is undefined and @i{wior} -is an implementation-defined I/O result code."" +is a non-zero I/O result code."" a_addr = (Cell *)malloc(u?u:1); wior = IOR(a_addr==NULL); free ( a_addr -- wior ) memory ""Return the region of data space starting at @i{a-addr} to the system. -The regon must originally have been obtained using @code{allocate} or +The region must originally have been obtained using @code{allocate} or @code{resize}. If the operational is successful, @i{wior} is 0. -If the operation fails, @i{wior} is an implementation-defined -I/O result code."" +If the operation fails, @i{wior} is a non-zero I/O result code."" free(a_addr); wior = 0; @@ -1540,8 +1557,8 @@ resize ( a_addr1 u -- a_addr2 wior ) mem ""Change the size of the allocated area at @i{a-addr1} to @i{u} address units, possibly moving the contents to a different area. @i{a-addr2} is the address of the resulting area. -If the operational is successful, @i{wior} is 0. -If the operation fails, @i{wior} is an implementation-defined +If the operation is successful, @i{wior} is 0. +If the operation fails, @i{wior} is a non-zero I/O result code. If @i{a-addr1} is 0, Gforth's (but not the Standard) @code{resize} @code{allocate}s @i{u} address units."" /* the following check is not necessary on most OSs, but it is needed @@ -1584,10 +1601,6 @@ wior = IOR(fclose((FILE *)wfileid)==EOF) open-file ( c_addr u ntype -- wfileid wior ) file open_file wfileid = (Cell)fopen(tilde_cstr(c_addr, u, 1), fileattr[ntype]); -#if defined(GO32) && defined(MSDOS) -if(wfileid && !(ntype & 1)) - setbuf((FILE*)wfileid, NULL); -#endif wior = IOR(wfileid == 0); create-file ( c_addr u ntype -- wfileid wior ) file create_file @@ -1595,10 +1608,6 @@ Cell fd; fd = open(tilde_cstr(c_addr, u, 1), O_CREAT|O_TRUNC|ufileattr[ntype], 0666); if (fd != -1) { wfileid = (Cell)fdopen(fd, fileattr[ntype]); -#if defined(GO32) && defined(MSDOS) - if(wfileid && !(ntype & 1)) - setbuf((FILE*)wfileid, NULL); -#endif wior = IOR(wfileid == 0); } else { wfileid = 0; @@ -1748,15 +1757,15 @@ d = r; #endif f! ( r f_addr -- ) float f_store -"" Store the floating-point value @i{r} to address @i{f-addr}."" +""Store @i{r} into the float at address @i{f-addr}."" *f_addr = r; f@ ( f_addr -- r ) float f_fetch -"" Fetch floating-point value @i{r} from address @i{f-addr}."" +""@i{r} is the float at address @i{f-addr}."" r = *f_addr; df@ ( df_addr -- r ) float-ext d_f_fetch -"" Fetch the double-precision IEEE floating-point value @i{r} from the address @i{df-addr}."" +""Fetch the double-precision IEEE floating-point value @i{r} from the address @i{df-addr}."" #ifdef IEEE_FP r = *df_addr; #else @@ -1764,7 +1773,8 @@ r = *df_addr; #endif df! ( r df_addr -- ) float-ext d_f_store -"" Store the double-precision IEEE floating-point value @i{r} to the address @i{df-addr}."" +""Store @i{r} as double-precision IEEE floating-point value to the +address @i{df-addr}."" #ifdef IEEE_FP *df_addr = r; #else @@ -1772,7 +1782,7 @@ df! ( r df_addr -- ) float-ext d_f_store #endif sf@ ( sf_addr -- r ) float-ext s_f_fetch -"" Fetch the single-precision IEEE floating-point value @i{r} from the address @i{sf-addr}."" +""Fetch the single-precision IEEE floating-point value @i{r} from the address @i{sf-addr}."" #ifdef IEEE_FP r = *sf_addr; #else @@ -1780,7 +1790,8 @@ r = *sf_addr; #endif sf! ( r sf_addr -- ) float-ext s_f_store -"" Store the single-precision IEEE floating-point value @i{r} to the address @i{sf-addr}."" +""Store @i{r} as single-precision IEEE floating-point value to the +address @i{sf-addr}."" #ifdef IEEE_FP *sf_addr = r; #else @@ -1821,12 +1832,11 @@ fnip ( r1 r2 -- r2 ) gforth f_nip ftuck ( r1 r2 -- r2 r1 r2 ) gforth f_tuck float+ ( f_addr1 -- f_addr2 ) float float_plus -"" Increment @i{f-addr1} by the number of address units corresponding to the size of -one floating-point number, to give @i{f-addr2}."" +""@code{1 floats +}."" f_addr2 = f_addr1+1; floats ( n1 -- n2 ) float -""@i{n2} is the number of address units corresponding to @i{n1} floating-point numbers."" +""@i{n2} is the number of address units of @i{n1} floats."" n2 = n1*sizeof(Float); floor ( r1 -- r2 ) float @@ -2021,24 +2031,24 @@ r2 = atanh(r1); r> IF fnegate THEN ; sfloats ( n1 -- n2 ) float-ext s_floats -""@i{n2} is the number of address units corresponding to @i{n1} +""@i{n2} is the number of address units of @i{n1} single-precision IEEE floating-point numbers."" n2 = n1*sizeof(SFloat); dfloats ( n1 -- n2 ) float-ext d_floats -""@i{n2} is the number of address units corresponding to @i{n1} +""@i{n2} is the number of address units of @i{n1} double-precision IEEE floating-point numbers."" n2 = n1*sizeof(DFloat); sfaligned ( c_addr -- sf_addr ) float-ext s_f_aligned -"" @i{sf-addr} is the first single-float-aligned address greater +""@i{sf-addr} is the first single-float-aligned address greater than or equal to @i{c-addr}."" sf_addr = (SFloat *)((((Cell)c_addr)+(sizeof(SFloat)-1))&(-sizeof(SFloat))); : [ 1 sfloats 1- ] Literal + [ -1 sfloats ] Literal and ; dfaligned ( c_addr -- df_addr ) float-ext d_f_aligned -"" @i{df-addr} is the first double-float-aligned address greater +""@i{df-addr} is the first double-float-aligned address greater than or equal to @i{c-addr}."" df_addr = (DFloat *)((((Cell)c_addr)+(sizeof(DFloat)-1))&(-sizeof(DFloat))); : @@ -2118,6 +2128,7 @@ lp -= sizeof(Float); *(Float *)lp = r; fpick ( u -- r ) gforth +""Actually the stack effect is @code{ r0 ... ru u -- r0 ... ru r0 }."" r = fp[u+1]; /* +1, because update of fp happens before this fragment */ : floats fp@ + f@ ; @@ -2248,3 +2259,48 @@ char newline[] = { }; c_addr=newline; u=sizeof(newline); +: + "newline count ; +Create "newline 1 c, $0A c, + +utime ( -- dtime ) gforth +""Report the current time in microseconds since some epoch."" +struct timeval time1; +gettimeofday(&time1,NULL); +dtime = timeval2us(&time1); + +cputime ( -- duser dsystem ) gforth +""duser and dsystem are the respective user- and system-level CPU +times used since the start of the Forth system (excluding child +processes), in microseconds (the granularity may be much larger, +however). On platforms without the getrusage call, it reports elapsed +time (since some epoch) for duser and 0 for dsystem."" +#ifdef HAVE_GETRUSAGE +struct rusage usage; +getrusage(RUSAGE_SELF, &usage); +duser = timeval2us(&usage.ru_utime); +dsystem = timeval2us(&usage.ru_stime); +#else +struct timeval time1; +gettimeofday(&time1,NULL); +duser = timeval2us(&time1); +dsystem = (DCell)0; +#endif + +v* ( f_addr1 nstride1 f_addr2 nstride2 ucount -- r ) gforth v_star +""dot-product: r=v1*v2. The first element of v1 is at f_addr1, the +next at f_addr1+nstride1 and so on (similar for v2). Both vectors have +ucount elements."" +for (r=0.; ucount>0; ucount--) { + r += *f_addr1 * *f_addr2; + f_addr1 = (Float *)(((Address)f_addr1)+nstride1); + f_addr2 = (Float *)(((Address)f_addr2)+nstride2); +} + +faxpy ( ra f_x nstridex f_y nstridey ucount -- ) gforth +""vy=ra*vx+vy"" +for (; ucount>0; ucount--) { + *f_y += ra * *f_x; + f_x = (Float *)(((Address)f_x)+nstridex); + f_y = (Float *)(((Address)f_y)+nstridey); +}