my_orders.c File Reference

#include <stdio.h>
#include <stdlib.h>

Go to the source code of this file.

Macros
#define	_FCALL

Functions
void	tri_direct (unsigned *data, unsigned *iwork, unsigned *index, int n, int irecl, unsigned *inds)
void	my_orders_ (int *mode, unsigned *iwork, unsigned *data, unsigned *index, int *n, int *irecl)
void _FCALL	MY_ORDERS (int *mode, int *iwork, int *data, int *index, int *n, int *irecl)
void	my_orders (int *mode, int *iwork, int *data, int *index, int *n, int *irecl)
void	my_orders__ (int *mode, int *iwork, int *data, int *index, int *n, int *irecl)

Macro Definition Documentation

_FCALL

#define _FCALL

Definition at line 26 of file my_orders.c.

Function Documentation

MY_ORDERS()

void _FCALL MY_ORDERS	(	int *	mode,
		int *	iwork,
		int *	data,
		int *	index,
		int *	n,
		int *	irecl )

Definition at line 79 of file my_orders.c.

{my_orders_(mode,(unsigned*)iwork,(unsigned*)data,(unsigned*)index,n,irecl);}

my_orders()

void my_orders	(	int *	mode,
		int *	iwork,
		int *	data,
		int *	index,
		int *	n,
		int *	irecl )

Definition at line 82 of file my_orders.c.

{my_orders_(mode,(unsigned*)iwork,(unsigned*)data,(unsigned*)index,n,irecl);}

my_orders_()

void my_orders_	(	int *	mode,
		unsigned *	iwork,
		unsigned *	data,
		unsigned *	index,
		int *	n,
		int *	irecl )

Definition at line 64 of file my_orders.c.

{
    int i;
    if(*mode == 0){
        for(i = 0 ;i < *n ; i++) index[i] = i ;
        tri_direct(data,iwork,index,*n,*irecl,index+*n);
        for(i = 0 ;i < *n ; i++) index[i] += 1 ;/* c => fortran */
    } else if(*mode == 10){
        for(i = 0 ;i < *n ; i++) index[i] -= 1 ;
        tri_direct(data,iwork,index,*n,*irecl,index+*n);
        for(i = 0 ;i < *n ; i++) index[i] += 1 ;/* c => fortran */
    } else {
        *mode = -1;
    }
}

my_orders__()

void my_orders__	(	int *	mode,
		int *	iwork,
		int *	data,
		int *	index,
		int *	n,
		int *	irecl )

Definition at line 85 of file my_orders.c.

{my_orders_(mode,(unsigned*)iwork,(unsigned*)data,(unsigned*)index,n,irecl);}

tri_direct()

void tri_direct	(	unsigned *	data,
		unsigned *	iwork,
		unsigned *	index,
		int	n,
		int	irecl,
		unsigned *	inds )

Definition at line 28 of file my_orders.c.

{
    /* tri a adressage direct                               */
    /*  0 <= x <= 2^16 - 1                                  */
 
    unsigned range=65535,mask=0xffff;
    int i,j,k;
 
    for(k=irecl-1;k >= 0;k--){
 
        /* poid faible (16bits)*/
        for(i = 0 ;i < range+2 ; i++) iwork[i] = 0 ;
        for(i = 0 ;i < n ; i++){
            iwork[(data[irecl*index[i]+k] & mask) + 1] += 1 ;
        }
        for(i = 0 ;i < range ; i++) iwork[i+1] += iwork[i];
        for(i = 0 ;i < n ; i++){
            j = data[irecl*index[i]+k] & mask;
            inds[iwork[j]] = index[i] ;
            iwork[j] += 1;
        }
 
        /* poid fort (16bits)*/
        for(i = 0 ;i < range+2 ; i++) iwork[i] = 0 ;
        for(i = 0 ;i < n ; i++) {
            iwork[((data[irecl*inds[i]+k] >> 16) & mask) + 1] += 1 ;
        }
        for(i = 0 ;i < range ; i++) iwork[i+1] += iwork[i];
        for(i = 0 ;i < n ; i++){
            j = (data[irecl*inds[i]+k] >> 16) & mask;
            index[iwork[j]] = inds[i] ;
            iwork[j] += 1;
        }
    }
 
}