OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
VE_Mumps_symqamd.c
Go to the documentation of this file.
1/*
2 * -----------------------------------------------------------------------
3 * NEC Vector Host MUMPS Library Interface
4 * -----------------------------------------------------------------------
5 */
6
7#include <VE_Mumps.h>
8
9/*
10 * ---------------------------------------------------------------------------
11 * VE_Mumps_symqamd
12 * ---------------------------------------------------------------------------
13 */
14
15void VE_Mumps_symqamd( int THRESH , int * NDENSE , int N , int TOTEL , int64_t IWLEN ,
16 int64_t * PE , int64_t * PFREE , int * LEN , int * IW , int * NV ,
17 int * ELEN , int * LAST , int * NCMPA , int * DEGREE , int * HEAD ,
18 int * NEXT , int * W , int * PERM , int * LISTVAR_SCHUR , int SIZE_SCHUR ,
19 int AGG6 , int * PARENT )
20{
21 int ret=1,iarg;
22 int64_t symid;
23 vhcall_handle handle;
24 vhcall_args *ca;
25 uint64_t retval;
26
27 char * env_metis_lib_name;
28 char metis_lib_name[ NECVH_MAX_LIBNAME ];
29
30 /*
31 * -----------------------------------------------------------------------
32 * Set the library name ( check environment )
33 * -----------------------------------------------------------------------
34 */
35
36 if ( env_metis_lib_name = getenv( "VH_MUMPS_LIBRARY" ) ) {
37 size_t lenv = strlen( env_metis_lib_name );
38 if ( strlen( env_metis_lib_name ) > (size_t)NECVH_MAX_LIBNAME ) {
39 fprintf( stdout , "%s @ %d failed : metis library name too long (%s)\n" , basename(__FILE__) , __LINE__ , env_metis_lib_name );
40 fflush( stdout );
41 (void) mumps_abort();
42 }
43 else {
44 memcpy( ( void * ) metis_lib_name , ( void * ) env_metis_lib_name , lenv * sizeof( char) );
45 metis_lib_name[ lenv ] = '\0';
46 fprintf( stdout , "%s @ %d input library %s\n" , basename(__FILE__) , __LINE__ , metis_lib_name );
47 fflush( stdout );
48 }
49 }
50 else {
51 strcpy( metis_lib_name , "libvh.so" );
52 }
53
54 /*
55 * -----------------------------------------------------------------------
56 * Load VH C library
57 * -----------------------------------------------------------------------
58 */
59
60 handle = vhcall_install( metis_lib_name );
61 if (handle == (vhcall_handle)-1) {
62 perror("vhcall_install");
63 fprintf( stdout , "%s @ %d failed : can't install library %s\n" , basename(__FILE__) , __LINE__ , metis_lib_name );
64 fprintf( stdout , "Please check your LD_LIBRARY_PATH variable\n");
65 fflush( stdout );
66 (void) mumps_abort();
67 }
68
69 /*
70 * -----------------------------------------------------------------------
71 * Find VH C library function
72 * -----------------------------------------------------------------------
73 */
74
75 symid = vhcall_find(handle, "mumps_symqamd_");
76 if ( symid == -1 ) {
77 fprintf( stdout , "%s @ %d failed : can't find symbol MUMPS_SYMQAMD\n" , basename(__FILE__) , __LINE__ );
78 perror("vhcall_find");
79 fflush( stdout );
80 (void) mumps_abort();
81 }
82
83 /*
84 * -----------------------------------------------------------------------
85 * Allocate arguments object for VH library function
86 * -----------------------------------------------------------------------
87 */
88
89 ca = vhcall_args_alloc();
90 if ( ! ca ) {
91 perror("vhcall_args_alloc");
92 fprintf( stdout , "%s @ %d failed : unable to allocate function's argument list\n" , basename(__FILE__) , __LINE__ );
93 fflush( stdout );
94 (void) mumps_abort();
95 }
96
97 /*
98 * -----------------------------------------------------------------------
99 * Set arguments
100 * -----------------------------------------------------------------------
101 */
102 iarg = 0;
103 ret = vhcall_args_set_pointer( ca , VHCALL_INTENT_IN , iarg++ , &THRESH , sizeof( int ) );
104 if ( ret ) {
105 perror("vhcall_args_set_pointer");
106 fprintf( stdout , "%s @ %d failed : unable to set argument.\n" , basename(__FILE__) , __LINE__ );
107 fflush( stdout );
108 (void) mumps_abort();
109 }
110
111 ret = vhcall_args_set_pointer( ca , VHCALL_INTENT_OUT , iarg++ , NDENSE , N * sizeof( int ) );
112 if ( ret ) {
113 perror("vhcall_args_set_pointer");
114 fprintf( stdout , "%s @ %d failed : unable to set argument.\n" , basename(__FILE__) , __LINE__ );
115 fflush( stdout );
116 (void) mumps_abort();
117 }
118
119 ret = vhcall_args_set_pointer( ca , VHCALL_INTENT_IN , iarg++ , &N , sizeof( int ) );
120 if ( ret ) {
121 perror("vhcall_args_set_pointer");
122 fprintf( stdout , "%s @ %d failed : unable to set argument.\n" , basename(__FILE__) , __LINE__ );
123 fflush( stdout );
124 (void) mumps_abort();
125 }
126
127 ret = vhcall_args_set_pointer( ca , VHCALL_INTENT_IN , iarg++ , &TOTEL , sizeof( int ) );
128 if ( ret ) {
129 perror("vhcall_args_set_pointer");
130 fprintf( stdout , "%s @ %d failed : unable to set argument.\n" , basename(__FILE__) , __LINE__ );
131 fflush( stdout );
132 (void) mumps_abort();
133 }
134
135 ret = vhcall_args_set_pointer( ca , VHCALL_INTENT_IN , iarg++ , &IWLEN , sizeof( int64_t ) );
136 if ( ret ) {
137 perror("vhcall_args_set_pointer");
138 fprintf( stdout , "%s @ %d failed : unable to set argument.\n" , basename(__FILE__) , __LINE__ );
139 fflush( stdout );
140 (void) mumps_abort();
141 }
142
143 ret = vhcall_args_set_pointer( ca , VHCALL_INTENT_INOUT , iarg++ , PE , N * sizeof( int64_t ) );
144 if ( ret ) {
145 perror("vhcall_args_set_pointer");
146 fprintf( stdout , "%s @ %d failed : unable to set argument.\n" , basename(__FILE__) , __LINE__ );
147 fflush( stdout );
148 (void) mumps_abort();
149 }
150
151 ret = vhcall_args_set_pointer( ca , VHCALL_INTENT_INOUT , iarg++ , PFREE , sizeof( int64_t ) );
152 if ( ret ) {
153 perror("vhcall_args_set_pointer");
154 fprintf( stdout , "%s @ %d failed : unable to set argument.\n" , basename(__FILE__) , __LINE__ );
155 fflush( stdout );
156 (void) mumps_abort();
157 }
158
159 ret = vhcall_args_set_pointer( ca , VHCALL_INTENT_INOUT , iarg++ , LEN , N * sizeof( int ) );
160 if ( ret ) {
161 perror("vhcall_args_set_pointer");
162 fprintf( stdout , "%s @ %d failed : unable to set argument.\n" , basename(__FILE__) , __LINE__ );
163 fflush( stdout );
164 (void) mumps_abort();
165 }
166
167 ret = vhcall_args_set_pointer( ca , VHCALL_INTENT_INOUT , iarg++ , IW , IWLEN * sizeof( int ) );
168 if ( ret ) {
169 perror("vhcall_args_set_pointer");
170 fprintf( stdout , "%s @ %d failed : unable to set argument.\n" , basename(__FILE__) , __LINE__ );
171 fflush( stdout );
172 (void) mumps_abort();
173 }
174
175 ret = vhcall_args_set_pointer( ca , VHCALL_INTENT_INOUT , iarg++ , NV , N * sizeof( int ) );
176 if ( ret ) {
177 perror("vhcall_args_set_pointer");
178 fprintf( stdout , "%s @ %d failed : unable to set argument.\n" , basename(__FILE__) , __LINE__ );
179 fflush( stdout );
180 (void) mumps_abort();
181 }
182
183 ret = vhcall_args_set_pointer( ca , VHCALL_INTENT_OUT , iarg++ , ELEN , N * sizeof( int ) );
184 if ( ret ) {
185 perror("vhcall_args_set_pointer");
186 fprintf( stdout , "%s @ %d failed : unable to set argument.\n" , basename(__FILE__) , __LINE__ );
187 fflush( stdout );
188 (void) mumps_abort();
189 }
190
191 ret = vhcall_args_set_pointer( ca , VHCALL_INTENT_OUT , iarg++ , LAST , TOTEL * sizeof( int ) );
192 if ( ret ) {
193 perror("vhcall_args_set_pointer");
194 fprintf( stdout , "%s @ %d failed : unable to set argument.\n" , basename(__FILE__) , __LINE__ );
195 fflush( stdout );
196 (void) mumps_abort();
197 }
198
199 ret = vhcall_args_set_pointer( ca , VHCALL_INTENT_OUT , iarg++ , NCMPA , sizeof( int ) );
200 if ( ret ) {
201 perror("vhcall_args_set_pointer");
202 fprintf( stdout , "%s @ %d failed : unable to set argument.\n" , basename(__FILE__) , __LINE__ );
203 fflush( stdout );
204 (void) mumps_abort();
205 }
206
207 ret = vhcall_args_set_pointer( ca , VHCALL_INTENT_OUT , iarg++ , DEGREE , N * sizeof( int ) );
208 if ( ret ) {
209 perror("vhcall_args_set_pointer");
210 fprintf( stdout , "%s @ %d failed : unable to set argument.\n" , basename(__FILE__) , __LINE__ );
211 fflush( stdout );
212 (void) mumps_abort();
213 }
214
215 ret = vhcall_args_set_pointer( ca , VHCALL_INTENT_OUT , iarg++ , HEAD , TOTEL * sizeof( int ) );
216 if ( ret ) {
217 perror("vhcall_args_set_pointer");
218 fprintf( stdout , "%s @ %d failed : unable to set argument.\n" , basename(__FILE__) , __LINE__ );
219 fflush( stdout );
220 (void) mumps_abort();
221 }
222
223 ret = vhcall_args_set_pointer( ca , VHCALL_INTENT_OUT , iarg++ , NEXT , N * sizeof( int ) );
224 if ( ret ) {
225 perror("vhcall_args_set_pointer");
226 fprintf( stdout , "%s @ %d failed : unable to set argument.\n" , basename(__FILE__) , __LINE__ );
227 fflush( stdout );
228 (void) mumps_abort();
229 }
230
231 ret = vhcall_args_set_pointer( ca , VHCALL_INTENT_OUT , iarg++ , W , N * sizeof( int ) );
232 if ( ret ) {
233 perror("vhcall_args_set_pointer");
234 fprintf( stdout , "%s @ %d failed : unable to set argument.\n" , basename(__FILE__) , __LINE__ );
235 fflush( stdout );
236 (void) mumps_abort();
237 }
238
239 ret = vhcall_args_set_pointer( ca , VHCALL_INTENT_INOUT , iarg++ , PERM , N * sizeof( int ) );
240 if ( ret ) {
241 perror("vhcall_args_set_pointer");
242 fprintf( stdout , "%s @ %d failed : unable to set argument.\n" , basename(__FILE__) , __LINE__ );
243 fflush( stdout );
244 (void) mumps_abort();
245 }
246
247 ret = vhcall_args_set_pointer( ca , VHCALL_INTENT_OUT , iarg++ , LISTVAR_SCHUR , ( SIZE_SCHUR ? SIZE_SCHUR : 1 ) * sizeof( int ) );
248 if ( ret ) {
249 perror("vhcall_args_set_pointer");
250 fprintf( stdout , "%s @ %d failed : unable to set argument.\n" , basename(__FILE__) , __LINE__ );
251 fflush( stdout );
252 (void) mumps_abort();
253 }
254
255 ret = vhcall_args_set_pointer( ca , VHCALL_INTENT_INOUT , iarg++ , &SIZE_SCHUR , sizeof( int ) );
256 if ( ret ) {
257 perror("vhcall_args_set_pointer");
258 fprintf( stdout , "%s @ %d failed : unable to set argument.\n" , basename(__FILE__) , __LINE__ );
259 fflush( stdout );
260 (void) mumps_abort();
261 }
262
263 ret = vhcall_args_set_pointer( ca , VHCALL_INTENT_IN , iarg++ , &AGG6 , sizeof( int ) );
264 if ( ret ) {
265 perror("vhcall_args_set_pointer");
266 fprintf( stdout , "%s @ %d failed : unable to set argument.\n" , basename(__FILE__) , __LINE__ );
267 fflush( stdout );
268 (void) mumps_abort();
269 }
270
271 ret = vhcall_args_set_pointer( ca , VHCALL_INTENT_OUT , iarg++ , PARENT , N * sizeof( int ) );
272 if ( ret ) {
273 perror("vhcall_args_set_pointer");
274 fprintf( stdout , "%s @ %d failed : unable to set argument.\n" , basename(__FILE__) , __LINE__ );
275 fflush( stdout );
276 (void) mumps_abort();
277 }
278
279 /*
280 * -----------------------------------------------------------------------
281 * Invoke VH C library function
282 * -----------------------------------------------------------------------
283 */
284
285 ret = vhcall_invoke_with_args( symid , ca , &retval );
286 if ( ret ) {
287 perror("vhcall_invoke_with_args");
288 fprintf( stdout , "%s @ %d failed : unable to invoke the symbol from library\n" , basename(__FILE__) , __LINE__ );
289 fflush( stdout );
290 (void) mumps_abort();
291 }
292
293 /*
294 * -----------------------------------------------------------------------
295 * Free args
296 * -----------------------------------------------------------------------
297 */
298
299 vhcall_args_free( ca );
300
301 /*
302 * -----------------------------------------------------------------------
303 * Uninstall the library
304 * -----------------------------------------------------------------------
305 */
306
307 if ( vhcall_uninstall( handle ) ) {
308 perror("vhcall_uninstall");
309 fprintf( stdout , "cvhmetis_nodend failed : unable to uninstall the library\n" , basename(__FILE__) , __LINE__ );
310 fflush( stdout );
311 (void) mumps_abort();
312 }
313}
mumps_abort
#define mumps_abort
Definition VE_Metis.h:25
NECVH_MAX_LIBNAME
#define NECVH_MAX_LIBNAME
Definition VE_Metis.h:23
VE_Mumps.h
VE_Mumps_symqamd
void VE_Mumps_symqamd(int THRESH, int *NDENSE, int N, int TOTEL, int64_t IWLEN, int64_t *PE, int64_t *PFREE, int *LEN, int *IW, int *NV, int *ELEN, int *LAST, int *NCMPA, int *DEGREE, int *HEAD, int *NEXT, int *W, int *PERM, int *LISTVAR_SCHUR, int SIZE_SCHUR, int AGG6, int *PARENT)
Definition VE_Mumps_symqamd.c:15
handle
void * handle
Definition cpp_python_funct.cpp:67
N
#define N
Definition example_user.c:10
PERM
#define PERM
Definition rad2rad_c.c:91