Subversion Repositories forwardcom

/*************************  trigonometric.as  *********************************

* Author:        Agner Fog

* date created:  2018-02-29

* last modified: 2020-04-24

* Version:       1.09

* Project:       ForwardCom example, assembly code

* Description:   Makes a table of sine, cosine, and tangent

*                Uses a simple loop without vectorization

*

* Link with libraries libc.li and math.li

*

* Copyright 2018-2020 GNU General Public License http://www.gnu.org/licenses

*****************************************************************************/

const section read ip                            // read-only data section

intro: int8 "\n   x        sin(x)        cos(x)        tan(x)",0  // table heading

form: int8 "\n%6.3f  %12.8f  %12.8f  %12.8f",0   // format string for printf

newline: int8 "\n",0                             // newline

const end

bss section datap uninitialized                  // uninitialized read/write data section

int64 parlist[5]                                 // parameter list for printf

bss end

code section execute align = 4                   // code section

extern _sincosf: function, reguse = 0, 0x7BF      // library function: sine and cosine in radians

extern _printf: function                         // library function: formatted output to stdout

extern _puts:   function                         // library function: print string to stdout

_main function public                            // program begins here

int64 r0 = address([intro])

call _puts

for (double v12 = -1; v12 < 4; v12 += 0.1) {     // loop

   double v0 = compress(v12, 0)

   float v0 = v0

   call _sincosf                                  // sin(x) in v0, cos(x) in v1

   double v0 = expand(v0, 0)

   double v1 = expand(v1, 0)

   double v2 = v0 / v1                           // tan(x) = sin(x) / cos(x)

   int64 r0 = address([form])                    // format string

   int64 r1 = address([parlist])                 // parameter list for _printf

   double [parlist,    scalar] = v12             // x

   double [parlist+8,  scalar] = v0              // sin(x)

   double [parlist+16, scalar] = v1              // cos(x)

   double [parlist+24, scalar] = v2              // tan(x)

   call _printf                                  // print formatted results

}

int64 r0 = address([newline])                    // end with newline

call _puts

int64 r0 = 0                                     // program return value

return                                           // return from main

_main end

code end