<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;"><DIV>Thank you so much for clear things up Stephan.  Since I may need to pass large (GB) size matrices, how can I pass the parameters entirely by reference (for saving memory) to a BLAS function?   </DIV>
<DIV> </DIV>
<DIV>Thanks again for the big help,</DIV>
<DIV> </DIV>
<DIV>Salem </DIV>
<DIV> </DIV>
<DIV>P.S.  My understanding (and I could be wrong) is that SAC does not support call by reference for calling a SAC function.  Then in actual C translation, can SAC compiler clever enough to pass a read-only parameter by reference?</DIV>
<DIV><BR><BR>--- On <B>Fri, 8/6/10, <SPAN>Herhut</SPAN>, Stephan <I><s.a.herhut@herts.ac.uk></I></B> wrote:<BR></DIV>
<BLOCKQUOTE style="PADDING-LEFT: 5px; MARGIN-LEFT: 5px; BORDER-LEFT: rgb(16,16,255) 2px solid"><BR>From: Herhut, Stephan <s.a.herhut@herts.ac.uk><BR>Subject: Re: [sac-user] English Documentation for Calling C from SAC<BR>To: "Salem Reyen" <salemreyen@yahoo.com>, "sac-user@sac-home.org" <sac-user@sac-home.org><BR>Date: Friday, August 6, 2010, 10:23 AM<BR><BR>
<DIV id=yiv530853842><FONT face="Calibri, Verdana, Helvetica, Arial"><SPAN style="FONT-SIZE: 11pt">Dear Salem,<BR><BR>unfortunately no such documentation exists at the moment. There is a German master’s thesis but that is outdated, too. <BR><BR>We have not tried to use a BLAS C implementation in SAC. However, we have implemented some of the BLAS routines in SAC itself. It would be great to see how the two approaches compare. From our experience, if one compares single primitives, BLAS C implementations will almost always win as they have been tuned to the extreme. The power of an optimising compiler like sac2c only plays out as soon as multiple primitives are called in sequence. Only then additional optimisation potential arises. <BR><BR>You can find our BLAS implementation in the standard library (modules/numerical/blas).<BR><BR>Let my try to give you a brief intro into the C interface. This should suffice to integrate BLAS. As an example, consider
 the string module in the standard library (modules/structures/String.sac). In that file, the SAC function strtoi is defined as follows: <BR><BR>external int, string strtoi(string S, int BASE);<BR><BR>As the keyword external gives away, the function is not implemented in SAC but some other language. The compiler does not actually care what language that is, as long as it uses C calling conventions. To declare an external functions, one simply provides its signature prefixed with the keyword external. <BR>By just using the line above, the SAC compiler will make some assumptions on how to map the function to an actual C function. Firstly, it will map return values and arguments from the SAC side to arguments only on the C side. This is done to map the multiple return values of SAC to the C word. So, for the above example, the resulting C prototype would be<BR><BR>void strtoi( int *res1, char **res2, char *S, int BASE);<BR><BR>As you can see, the return
 values have been mapped into reference arguments. Often, this is not desired, like in the example of BLAS where the signatures are fixed. For those cases, we provide a means to remap SAC results/arguments to C results/arguments. For the above example, the following pragma is used in the standard library:<BR><BR>    #pragma linksign [0,1,2,3]<BR><BR>Pragmas in general directly follow an external function declaration. In this case, we use the linksign pragma. It maps each result/argument in SAC to the corresponding C result/arguments by means of a permutation vector. So, to achieve the desired C signature of strtoi, we want to have the first SAC result to be mapped to the C result and the arguments to be mapped in the same order. This is done using the permutation vector [0,1,2,3].  Position zero thereby represents the C function return value and numbers from one onwards are the function arguments. Using that vector, we
 get<BR><BR>int strtoi( char **res1, char *S, int BASE)<BR><BR>which is the standard C signature for that function. As you might guess, SAC by default uses the permutation vector [1,2,...].<BR><BR>There are some more tuneable settings:<BR><BR>    #pragma linkname " SACstrtoi"<BR><BR>The linkname pragma allows to use a different C function name. Using the above pragma, we would get<BR><BR>int SACstrtoi( char **res1, char *S, int BASE);<BR><BR>Further two pragmas handle linking:<BR><BR>    #pragma linkobj "src/String/strtoi.o"<BR><BR>would instruct the compiler to link the given file whenever the function is required in a program. Similarly<BR><BR>    #pragma linkobj "mylib"<BR><BR>would make the compiler link the final executable against the specified library. The library is looked for in the search path configured in the sac2crc resource file. The default paths for system wide libraries are
 preconfigured. By editing or creating a .sac2crc file in your home directory and adding the following lines<BR><BR>target default:<BR>EXTLIBPATH       := “path1:path2:pathn”<BR><BR>you can add additional paths to the search. A description of all valid field in that file can be found in its master copy in the setup subdirectory of the sac2c compiler. <BR><BR>Lastly, there is a the<BR><BR>#pragma refcounting [0,1,2]<BR><BR>pragma. This configures whether parameters are passed by value or reference. Per default, all parameters are passed by value and thus are immutable. If you want to modify a parameter, you have to use a slightly more complex calling convention. I will cover that subject as demand arises. <BR><BR>I hope this suffices to get you started. If you have any further questions, please feel free to reply to this list. Please use the list instead of contacting me directly as this will increase the audience and
 ensure a timely reply<BR><BR>Regards<BR>  Stephan<BR><BR>On 05/08/2010 04:24, "Salem Reyen" <<A rel=nofollow>salemreyen@yahoo.com</A>> wrote:<BR><BR></SPAN></FONT>
<BLOCKQUOTE><FONT face="Calibri, Verdana, Helvetica, Arial"><SPAN style="FONT-SIZE: 11pt">Hi,<BR> <BR>  Is there any documentation (in English) about how to call C functions from a SAC Module?  More importantly, has anyone successfully tried to call BLAS/LAPACK (eg, Intel MKL or AMD AML) from SAC?  Any help will be greatly appreciated.<BR> <BR>Salem<BR><BR></SPAN></FONT></BLOCKQUOTE><FONT face="Calibri, Verdana, Helvetica, Arial"><SPAN style="FONT-SIZE: 11pt"><BR>-- <BR>Stephan Herhut<BR>Centre for Computer Science and Informatics Research<BR>Science and Technology Research Institute<BR>University of Hertfordshire<BR><BR></SPAN></FONT></DIV></BLOCKQUOTE></td></tr></table><br>