Simple Example

The following C source code calls a Fortran function to multiply the input vector by 10.

/*
////////////////////////////////////////////////////////////////////
// File      : simple.c
// Date      : 7/14/98
// Copyright : Altair Software Development, Inc. 1998
//
// Calculate nothing and return the same vector that you were given.
// freeman
//////////////////////////////////////////////////////////////////////
//                                                                  //
// This source code intended to serve as an example of              //
// the external function capability of MotionView/MotionView.       //
// Neither Altair nor the authors assume any responsibility         //
// for the validity, accuracy, or applicability of any              //
// results obtained from this program. The user must verify         //
// the results on his/her own.                                      //
//                                                                  //
//////////////////////////////////////////////////////////////////////
// arg1 = dummy vector
// arg2 = vector to return
////////////////////////////////////////////////////////////////////
*/
#include <stdio.h>
#include <stdlib.h>
#include <signal.h>
#include <sys/types.h>
/*
#include <sys/socket.h>
#include <unistd.h>
*/
#if defined(SGI_PLATFORM)
#include <strings.h>
#else
#include <string.h>
#endif

#include "../base/altair_ipc.h"

#if defined(_WINDOWS)
extern void __stdcall SIMPLE_S_(); 
#endif

int main(argc, argv)
        int argc;
        char *argv[];
{
    AltairIpcVectorList  vec_list;
    char                 msg_buff[1024];
    double              *return_vec;
    double              *in_vec1;
    double              *in_vec2;
    int                  nfp=1;
    int                  nlp;

    /*
    * You must do this before anything else...
    */
    if (AltairIpcParseCommandLine(argc, argv) != ALTAIR_IPC_SUCCESS)
    {
        AltairIpcPutError("Could not get data.");
    }
    if (AltairIpcGetVectorList(&vec_list) != ALTAIR_IPC_SUCCESS)
    {
        AltairIpcPutError("Could not get data.");
    }
    /*
    * Do some program/routine specific checking...
    */
    AltairIpcPutMessage("Checking for validity of data.");
    if (vec_list.num_vectors < 2)
    {
        sprintf(msg_buff, "Expecting two arguments, got %d.", vec_list.num_vectors);
        AltairIpcPutError(msg_buff);
    }
    else
    {
        AltairIpcPutMessage("Passed validity check.");
        in_vec1 = vec_list.vec_val[0];
        in_vec2 = vec_list.vec_val[1];
        nlp = vec_list.vec_cnt[0];
        sprintf(msg_buff, "delta for first vector = %f", ((in_vec1[(nlp-1)]-
        n_vec1[0])/(nlp-1)));
        AltairIpcPutMessage(msg_buff);
		/* get some memory for the return vector */
        return_vec = malloc(nlp*sizeof(double));
        if ( !return_vec )
        {
            AltairIpcPutError("'malloc(...)' failed");
        }
        /* call the FORTRAN subroutine simple_s from this C program */
#if defined(HP_PLATFORM) || defined(IBM_PLATFORM)
        simple_s( in_vec2, return_vec, &nfp, &nlp);
#elif defined(_WINDOWS)
        SIMPLE_S_( in_vec2, return_vec, &nfp, &nlp);
#else
		simple_s_( in_vec2, return_vec, &nfp, &nlp);
#endif
		AltairIpcPutVector(nlp, return_vec);
        free(return_vec);
    }
    return 0;
}

	SUBROUTINE SIMPLE_S( INVEC, OUTVEC, NFP, NLP)

	! Parameters ___________________________________________________________

	INTEGER NFP ! Index of first data point

	INTEGER NLP ! Index of last data point

	DOUBLEPRECISION INVEC(NFP:NLP)

	DOUBLEPRECISION OUTVEC(NFP:NLP) 
	

	! Variables 

	INTEGER I

	DO I=NFP,NLP,1
		OUTVEC(I) = INVEC(I)*10
	END DO

	RETURN
	END