5. Examples

This section give you a tutorial of how to use the Ninf-G system for programming on the Grid. Simplicity of programming is the most beneficial aspect of the Ninf-G system, and we hope that users will be able to gridify his programs easily after reading this document. We hope to extend this example further to cover more advanced Ninf-G features. Examples are provided for GRPC API.


Gridifying a Numerical Library with Grid RPC API

We first cover the simple case where the library to be Gridifyied is defined as a linkable library function. Below is a sample code of a simple matrix multiply. The first scalar argument specifies the size of the matrix (n by n), parameters a and b are references to matrices to be multiplied, and c is the reference to the result matrix. Notice that, 1) the matrix (defined as arrays) do not itself embody size as type information, and 2) as a result there is a dependency between n and a, b, c. In fact, since array arguments are passed as a reference, one must assume the contents of the array are implicitly shared by the caller and the callee, with arbitrary choices as to using them as input, output, or temporary data structures.

void mmul(int n, double * a, double * b, double * c)
    double t;
    int i, j, k;
    for (i = 0; i < n; i++) {
        for (j = 0; j < n; j++) {
            t = 0;
            for (k = 0; k < n; k++){
                t += a[i * n + k] * b[k * n + j];
            c[i * n + j] = t;

The main routine which calls mmul() might be as follows:

    double A[N*N], B[N*N], C[N*N];

    initMatA(N, A); /* initialize */
    initMatB(N, B); /* initialize */

    mmul(N, A, B, C);

In order to "Gridify", or more precisely, allow mmul to be called remotely via GridRPC, we must describe the interface of the function so that information not embodied in the language type system becomes sufficiently available to the GridRPC system to make the remote call. Although future standardization is definitely conceivable, currently each GridRPC system has its own IDL (Interface Description Language); for example, Ninf has its own NinfIDL definition. Below we give the interface of mmul() defined by the NinfIDL syntax:

1: Module mmul;
3: Define mmul(IN int N, IN double A[N*N], 
4:             IN double B[N*N], OUT double C[N*N])
5: "matmul"
6: Required "mmul_lib.o"
7: Calls "C" mmul(N, A, B, C);

Line 1 declares the module name to be defined. There is a one-to-one correspondence between a module and an IDL file, and each module can have multiple entries to gridify multiple functions. Lines 3-7 are the definition for a particular entry mmul/mmul. Here, lines 3 and 4 declare the interface of the entry. The difference between a NinfIDL entry definition and the C prototype definition is that there are no return values (the return value of the Ninf call is used to return status info), argument input/output modes are specified, and array sizes are described in terms of the scalar arguments.

We note here that NinfIDL has special features to efficiently support gridifying of a library (similar features are found in Netsolve IDL). In contrast to standard procedure calls within a shared memory space, GridRPC needs to transfer data over the network. Transferring the entire contents of the array will be naturally very costly, especially for huge matrices appearing in real applications. Here, one will quickly observe that surprising number of numerical libraries take for granted the fact that address space of data structures, in particular arrays are shared, and (a) only use subarrays of the passed arrays, (b) write back results in the passed arrays, and (c) pass arrays as scratchpad temporaries. The portion of the arrays to be operated, etc., are determined by the semantics of the operation according to the input parameters passed to the function. For example in mmul, the whole arrays need to be passed, and their sizes are all N by N, where N is the first scalar parameter; A and B only need to be passed as input parameters and their contents do not change, while C is used as a return argument and thus need not be shipped to the server, but the result needs to be shipped back to the client. In general, in order to determine and minimize the size of transfer, NinfIDL allows flexible description of the array shape used by the remote library. One can specify leading dimensions, subarrays, and strides. In fact arbitrary arithmetic expressions involving constants and scalar arguments can be used in the array size expressions.

Line 5 is the comment describing the entry, while line 6 specifies the necessary object file when the executable for the particular file is to be linked. Line 7 gives the actual library function to be called, and the calling sequence; here "C" denotes C-style (row-major) array layout.

The user compiles this IDL file using the Ninf IDL compiler, and generates the stub code and its makefile. By executing this makefile a Ninf executable is generated. The user will subsequently register the executable to the server using the registry tool.

Now the client us ready to make the call of the network. In order to make a GridRPC call, the user modifies his original main program in the following manner. We notice that only the function call is modified---No need to change the program to adjust to the skeleton that the IDL generator generates as is with typical RPC systems such as CORBA. Moreover, we note that the IDL, the stub files and the executables are only resident on the server side, and the client only needs to link his program with a generic Ninf client library.

    double A[N*N], B[N*N], C[N*N];
    grpc_function_handle_t handle;
    initMatA(N, A); /* initialize */
    initMatB(N, B); /* initialize */
    grpc_function_handle_default(&handle, "mmul/mmul");
    if (grpc_call(&handle, N, A, B, C) != GRPC_NO_ERROR) {
        fprintf(stderr, "Error in grpc_call\n");

Gridifying Programs that use Files

The above example assumes that the numerical routine is supplied as a library with well-defined function API, or at least its source is available in a way such that it could easily converted into a library. In practice, many numerical routines are only available in a non-library executable and/or binary form, with input/output interfaces using files. In order to gridify such "canned" applications, GridRPC systems typically support remote files and their automatic management/transfer.

We take gnuplot as an example. Gnuplot in non-interactive mode inputs script from a specified file, and outputs the resulting graph to the standard output. Below is an example gnuplot script.

set terminal postscript
set xlabel "x"
set ylabel "y"
plot f(x) = sin(x*a), a = .2, f(x), a = .4, f(x)

If this script is saved under a filename "gplot":

> gnuplot gplot > graph.ps

will store the postscript representation of the graph to the file graph.ps. In order to execute gnuplot remotely, we must package it appropriately, and moreover must automatically transfer the input (gplot) and output (graph.ps) files between the client and the server.

Ninf-G IDL provides a type filename to specify that the particular argument is a file. Below is an example of using gnuplot via GridRPC.

Module plot;
Define plot(IN filename plotfile, OUT filename psfile )
"invoke gnuplot"
    char buffer[1000];
    sprintf(buffer, "gnuplot %s > %s", plotfile, psfile);

The IDL writes the string command sequence to invoke gnuplot into a variable buffer[], and invokes gnuplot as a system library. The file specified as an input is automatically transferred to the temporary directory of the server, and its temporary file name is passed to the stub function. As for the output file, only the temporary file name is created and passed to the stub function. After the stub program is executed, the files in output mode as specified in the IDL are automatically transferred to the client, and saved there under the name given in the argument.

Below is an example of how this function might be called via GridRPC.

#include <stdio.h>
#include "grpc.h"

main(int argc, char **argv)
    grpc_function_handle_t handle;


    grpc_function_handle_default(&handle, "plot/plot");

    if (grpc_call(&handle, argv[2], argv[3]) != GRPC_NO_ERROR) {
        fprintf(stderr, "Error in grpc_call\n");


We also note that, by combining this feature with the technique of using multiple servers simultaneously described in the next section, we can process large amount of data at once.

Using Multiple Servers for Parallel Programming on the Grid --- The Parameter Sweep Survey Example.

GridRPC can serve as a task-parallel programming abstraction, whose programs can scale from local workstations to the Grid. Here, we take an example of simple parameter sweep survey, and investigate how it can be easily programmed using GridRPC.

Calculating PI using a simple Monte Carlo Method

As an example, we compute the value of PI using a simple Monte Carlo Method. We generate a large number of random points within the square region that exactly encloses a unit circle (actually, 1/4 of a circle). We calculate the value of PI by inverse computing the area of the circle according to the probability that the points will fall within the circle. The program below shows the original sequential version.

long pi_trial(int seed, long times)
    long l, long counter = 0;

    for (l = 0; l < times; l++){
        double x = (double)random() / RAND_MAX;
        double y = (double)random() / RAND_MAX;
        if (x * x + y * y < 1.0)
    return counter;

main(int argc, char **argv)
    double pi;
    long times = atol(argv[1]);
    count = pi_trial(10, times);
    pi = 4.0 * (count / (double) times);
    printf("PI = %f\n", pi);

Gridifying the PI program.

First, we rewrite the program so that it does the appropriate GridRPC calls. The following steps are needed:

  1. Separate out the pi_trail() function into a separate file (say, trial_pi.c), and create its object file trial_pi.o using a standard C compiler.
  2. Describe the interface of pi_trial in an IDL file.
    Module pi;
    Define pi_trial(IN int seed, IN long times, OUT long * count)
    "monte carlo pi computation"
    Required "pi_trial.o"
        long counter;
        counter = pi_trial(seed, times);
        *count = counter;
  3. Rewrite the main program so that it makes a GridRPC call.
    main(int argc, char **argv)
        double pi;
        long times, count;
        grpc_function_handle_t handle;
        times = atol(argv[2]);
        grpc_function_handle_default(&handle, "pi/pi_trial");
        if (grpc_call(&handle, 10, times, &count) != GRPC_NO_ERROR) {
            fprintf(stderr, "Failed in grpc_call\n");
        pi = 4.0 * ( count / (double) times);
        printf("PI = %f\n", pi);

We now have made the body of the computation remote. The next phase is to parallelise it.

Employing Multiple Servers for Task Parallel Computation.

We next rewrite the main program so that parallel tasks are distributed to multiple servers. Although distribution of tasks are possible using metaserver scheduling with Ninf (and Agents with Netsolve), it is sometimes better to specify a host explicitly for performance reasons, for low overhead and explicit load balancing. Ninf-G allows explicit specification of servers by specifying the hostname in the initialization of the function handle.

The standard grpc_call() RPC is synchronous in that the client waits until the completion of the computation on the server side. For task-parallel execution, Ninf-G facilitates several asynchronous call APIs. For example, the most basic asynchronous call grpc_call_async is almost identical to grpc_call except that it returns immediately after all the arguments have been sent. The return value is the session ID of the asynchronous call; the ID is used for various synchronizations such as waiting for the return value of the call.

There are several calls for synchronization. The most basic is the grpc_wait(grpc_sessionid_t ID), where we wait for the result of the asynchronous call with the supplied session ID. grpc_wait_all() waits for all preceding asynchronous invocations made. Here, we employ grpc_wait_all() to parallelize the above PI client so that it uses multiple simultaneous remote server calls:

     1  #include "grpc.h"
     2  #define NUM_HOSTS 8
     3  char * hosts[] = {"node0.example.org", "node1.example.org", "node2.example.org", "node3.example.org",
     4                    "node4.example.org", "node5.example.org", "node6.example.org", "node7.example.org"};
     6  grpc_function_handle_t handles[NUM_HOSTS];
     7  grpc_sessionid_t ids[NUM_HOSTS];
     9  main(int argc, char **argv){
    10      double pi;
    11      long times, count[NUM_HOSTS], sum;
    12      char * config_file;
    13      int i;
    14      if (argc < 3) {
    15          fprintf(stderr, "USAGE: %s CONFIG_FILE TIMES \n", argv[0]);
    16          exit(2);
    17      }
    18      config_file = argv[1];
    19      times = atol(argv[2]) / NUM_HOSTS;
    21      /* Initialize GRPC runtimes. */
    22      if (grpc_initialize(config_file) != GRPC_NO_ERROR) {
    23          exit(2);
    24      }
    25      /* Initialize handles. */
    26      for (i = 0; i < NUM_HOSTS; i++)
    27          grpc_function_handle_init(&handles[i], hosts[i], "pi/pi_trial");
    29      for (i = 0; i < NUM_HOSTS; i++) {
    30          /* Parallel non-blocking remote function invocation. */
    31          if (grpc_call_async(&handles[i], &ids[i], i, times, &count[i]) != GRPC_NO_ERROR){
    32              grpc_perror_np("pi_trial");
    33              exit(2);
    34          }
    35      }
    36      /* Sync. */
    37      if (grpc_wait_all() != GRPC_NO_ERROR) {
    38          grpc_perror_np("wait_all");
    39          exit(2);
    40      }
    42      for (i = 0; i < NUM_HOSTS; i++)
    43          grpc_function_handle_destruct(&handles[i]);
    45      /* Compute and display pi. */
    46      for (i = 0, sum = 0; i < NUM_HOSTS; i++)
    47          sum += count[i];
    48      pi = 4.0 * ( sum / ((double) times * NUM_HOSTS));
    49      printf("PI = %f\n", pi);
    51      /* Finalize GRPC runtimes. */
    52      grpc_finalize();
    53 }

We specify the number of server hosts and their names in lines 2 and 3-4, respectively. Line 6 is the port number used, and line 19 divides the number of trials with the number of servers, determining the number of trials per server. The for loop in lines 29-35 invokes the servers asynchronously. Line 47 aggregates the results returned from all the servers.

In this manner, we can easily write a parallel parameter sweep survey program using the task parallel primitives of GridRPC. We next modify the program to perform dynamic load balancing.

last update : $Date: 2005/07/07 11:06:40 $