1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
|
/* Copyright (c) 1997 by Inria Lorraine. All Rights Reserved */
/***
NAME
sci_tools
PURPOSE
NOTES
HISTORY
fleury - Dec 17, 1997: Created.
$Log: sci_tools.c,v $
Revision 1.5 1998/03/27 12:20:22 fleury
Version pvm OK.
TODO: faire des tests de compil sur plateforme separee (POPC0
TODO: commenter source (-;
TODO: faire un peu de netoyage
Revision 1.4 1998/03/17 11:49:33 fleury
Broadcast OK.
TODO: mettre les listes.
TODO: faire qcq tests
Revision 1.3 1998/03/13 13:57:07 fleury
Version send/recv avec pack. A tester.
TODO: ajouter les listes + BROADCAST
TODO: faire un clean du dir et des fichiers...
Revision 1.2 1998/01/06 13:23:48 fleury
Use memcopy instead of fori
Revision 1.1 1997/12/18 18:35:57 fleury
Premier commit
Ajout de fct permettant de tester le type d une variable scilab
conversion de complex format scilab to complex format f77
TODO:use memcpy
TODO:use imatrix
***/
#include "../machine.h"
#ifndef __ABSC__
#include <malloc.h>
#else
#include <stdlib.h>
#endif
#include <stdio.h>
#include <string.h>
#include "sci_tools.h"
#ifdef __STDC__
void
C2F(ccomplexf)(int *n, double **ip, double *op)
#else
void
C2F(ccomplexf)(n, ip, op)
int *n;
double **ip;
double *op;
#endif
{
memcpy(op, *ip, *n * sizeof(double));
/* int i */
/* for (i = *n; --i >= 0; ) { */
/* op[i] = (*ip)[i];*/ /* TODO: replace by memcpy */
/* } */
SET_TYPE_COMPLEX(op); /* type is complex */
SET_NB_ROW(op, NB_ROW(op) / 2); /* nb row is halfed */
free((char*) (*ip));
} /* ccomplexf */
#ifdef __STDC__
void
(SciToF77)(double *ptr, int size, int lda)
#else
void
SciToF77(ptr, size, lda)
double *ptr;
int size;
int lda;
#endif
{
int i;
double *tab;
if ((tab = (double *) malloc(size * sizeof(double))) == NULL) {
(void) fprintf(stderr, "SciToF77: Error malloc\n");
return;
}
/* for (i = size; --i >= 0; ) { */
/* tab[i] = ptr[i]; */
/* } */
memcpy(tab, ptr, size * sizeof(double));
for (i = 0; i < size; ++i) {
ptr[2*i] = tab[i];
ptr[2*i+1] = ptr[lda+i];
}
free(tab);
} /* SciToF77 */
#ifdef __STDC__
void
(F77ToSci)(double *ptr, int size, int lda)
#else
void
F77ToSci(ptr, size, lda)
double *ptr;
int size;
int lda;
#endif
{
int i;
double *tab;
if ((tab = (double *) malloc(size * sizeof(double))) == NULL) {
(void) fprintf(stderr, "F77ToSci: Error malloc\n");
return;
}
for (i = 0; i < size; ++i) {
tab[i] = ptr[2*i+1];
ptr[i] = ptr[2*i];
}
memcpy(ptr + lda, tab, size * sizeof(double));
/* for (i = size; --i >= 0; ) { */
/* ptr[lda+i] = tab[i]; */
/* } */
free(tab);
} /* F77ToSci */
|
