p_Procs_Dynamic.cc
Go to the documentation of this file.
1 /****************************************
2 * Computer Algebra System SINGULAR *
3 ****************************************/
4 /***************************************************************
5  * File: p_Procs_Dynamic.cc
6  * Purpose: source for dynamically loaded version of p_Procs
7  * Author: obachman (Olaf Bachmann)
8  * Created: 12/00
9  *******************************************************************/
10 #include <misc/auxiliary.h>
11 
12 #include <factory/factory.h>
13 
14 #include <reporter/reporter.h>
15 
16 // #include <polys/structs.h>
17 #include <polys/monomials/ring.h>
19 
20 #include <polys/kbuckets.h>
21 
23 
28 
29 
30 #ifdef HAVE_DL
32 
33 #define WARN_MSG "Singular will work properly, but much slower."
34 
35 // need external linkage, so that dynl_sym works
36 #undef LINKAGE
37 #define LINKAGE extern "C"
38 #define p_Procs_Kernel
39 
40 #include "templates/p_Procs.inc"
41 
43 // include general p_Proc stuff
44 #include "templates/p_Procs_Impl.h"
45 
46 #include "mod_raw.h"
47 
48 // define to bound for length of p_Proc name
49 #define MAX_PROCNAME_LEN 200
50 
53 static void* p_procs_handle_FieldQ = NULL;
55 
56 static void* p_ProcInitHandle(void** handle, const char* module)
57 {
58  if (*handle == NULL)
59  {
60  char name[25];
61  sprintf(name, "p_Procs_%s", module);
62  *handle = dynl_open_binary_warn(name, WARN_MSG);
63  }
64  return *handle;
65 }
66 
67 static inline void* p_ProcGetHandle(p_Proc proc, p_Field field)
68 {
69  const char* module = p_ProcField_2_Module(proc, field);
70 
71  if (strcmp(module, "FieldIndep") == 0)
73  else if (strcmp(module, "FieldZp") == 0)
75  else if (strcmp(module, "FieldQ") == 0)
76  return p_ProcInitHandle(&p_procs_handle_FieldQ, module);
77  else if (strcmp(module, "FieldGeneral") == 0)
79  else
80  {
81  assume(0);
82  return NULL;
83  }
84 }
85 
86 
87 static void* GetGeneralProc(p_Proc proc)
88 {
89  switch(proc)
90  {
91  case p_Copy_Proc:
92  return cast_A_to_vptr(p_Copy__FieldGeneral_LengthGeneral_OrdGeneral);
93  case p_Delete_Proc:
94  return cast_A_to_vptr(p_Delete__FieldGeneral_LengthGeneral_OrdGeneral);
96  return cast_A_to_vptr(p_ShallowCopyDelete__FieldGeneral_LengthGeneral_OrdGeneral);
97  case p_Add_q_Proc:
98  return cast_A_to_vptr(p_Add_q__FieldGeneral_LengthGeneral_OrdGeneral);
99  case p_Neg_Proc:
100  return cast_A_to_vptr(p_Neg__FieldGeneral_LengthGeneral_OrdGeneral);
101  case p_Merge_q_Proc:
102  return cast_A_to_vptr(p_Merge_q__FieldGeneral_LengthGeneral_OrdGeneral);
103  case p_kBucketSetLm_Proc:
104  return cast_A_to_vptr(p_kBucketSetLm__FieldGeneral_LengthGeneral_OrdGeneral);
105 #ifdef HAVE_RINGS
106  case p_Mult_nn_Proc:
107  return cast_A_to_vptr(p_Mult_nn__RingGeneral_LengthGeneral_OrdGeneral);
108  case pp_Mult_nn_Proc:
109  return cast_A_to_vptr(pp_Mult_nn__RingGeneral_LengthGeneral_OrdGeneral);
110  case pp_Mult_mm_Proc:
111  return cast_A_to_vptr(pp_Mult_mm__RingGeneral_LengthGeneral_OrdGeneral);
113  return cast_A_to_vptr(pp_Mult_mm_Noether__RingGeneral_LengthGeneral_OrdGeneral);
114  case p_Mult_mm_Proc:
115  return cast_A_to_vptr(p_Mult_mm__RingGeneral_LengthGeneral_OrdGeneral);
117  return cast_A_to_vptr(p_Minus_mm_Mult_qq__RingGeneral_LengthGeneral_OrdGeneral);
119  return cast_A_to_vptr(pp_Mult_Coeff_mm_DivSelect__RingGeneral_LengthGeneral_OrdGeneral);
121  return cast_A_to_vptr(pp_Mult_Coeff_mm_DivSelectMult__RingGeneral_LengthGeneral_OrdGeneral);
122 #else
123  case p_Mult_nn_Proc:
124  return cast_A_to_vptr(p_Mult_nn__FieldGeneral_LengthGeneral_OrdGeneral);
125  case pp_Mult_nn_Proc:
126  return cast_A_to_vptr(pp_Mult_nn__FieldGeneral_LengthGeneral_OrdGeneral);
127  case pp_Mult_mm_Proc:
128  return cast_A_to_vptr(pp_Mult_mm__FieldGeneral_LengthGeneral_OrdGeneral);
130  return cast_A_to_vptr(pp_Mult_mm_Noether__FieldGeneral_LengthGeneral_OrdGeneral);
131  case p_Mult_mm_Proc:
132  return cast_A_to_vptr(p_Mult_mm__FieldGeneral_LengthGeneral_OrdGeneral);
134  return cast_A_to_vptr(p_Minus_mm_Mult_qq__FieldGeneral_LengthGeneral_OrdGeneral);
136  return cast_A_to_vptr(pp_Mult_Coeff_mm_DivSelect__FieldGeneral_LengthGeneral_OrdGeneral);
138  return cast_A_to_vptr(pp_Mult_Coeff_mm_DivSelectMult__FieldGeneral_LengthGeneral_OrdGeneral);
139 #endif
140  case p_Unknown_Proc:
141  break;
142  }
143  dReportBug("p_Unknown_Proc");
144  return NULL;
145 }
146 
147 #ifdef RDEBUG
148 #include <omalloc/omalloc.h>
149 
150 static const char* GetGeneralProcName(p_Proc proc)
151 {
152  switch(proc)
153  {
154  case p_Copy_Proc:
155  return "p_Copy__FieldGeneral_LengthGeneral_OrdGeneral";
156  case p_Delete_Proc:
157  return "p_Delete__FieldGeneral_LengthGeneral_OrdGeneral";
159  return "p_ShallowCopyDelete__FieldGeneral_LengthGeneral_OrdGeneral";
160  case p_Add_q_Proc:
161  return "p_Add_q__FieldGeneral_LengthGeneral_OrdGeneral";
162  case p_Neg_Proc:
163  return "p_Neg__FieldGeneral_LengthGeneral_OrdGeneral";
164  case p_Merge_q_Proc:
165  return "p_Merge_q__FieldGeneral_LengthGeneral_OrdGeneral";
166  case p_kBucketSetLm_Proc:
167  return "p_kBucketSetLm__FieldGeneral_LengthGeneral_OrdGeneral";
168  case p_Unknown_Proc:
169  break;
170 #ifdef HAVE_RINGS
171  case p_Mult_nn_Proc:
172  return "p_Mult_nn__RingGeneral_LengthGeneral_OrdGeneral";
173  case pp_Mult_nn_Proc:
174  return "pp_Mult_nn__RingGeneral_LengthGeneral_OrdGeneral";
175  case pp_Mult_mm_Proc:
176  return "pp_Mult_mm__RingGeneral_LengthGeneral_OrdGeneral";
178  return "pp_Mult_mm_Noether__RingGeneral_LengthGeneral_OrdGeneral";
179  case p_Mult_mm_Proc:
180  return "p_Mult_mm__RingGeneral_LengthGeneral_OrdGeneral";
182  return "p_Minus_mm_Mult_qq__RingGeneral_LengthGeneral_OrdGeneral";
184  return "pp_Mult_Coeff_mm_DivSelect__RingGeneral_LengthGeneral_OrdGeneral";
186  return "pp_Mult_Coeff_mm_DivSelectMult__RingGeneral_LengthGeneral_OrdGeneral";
187 #else
188  case p_Mult_nn_Proc:
189  return "p_Mult_nn__FieldGeneral_LengthGeneral_OrdGeneral";
190  case pp_Mult_nn_Proc:
191  return "pp_Mult_nn__FieldGeneral_LengthGeneral_OrdGeneral";
192  case pp_Mult_mm_Proc:
193  return "pp_Mult_mm__FieldGeneral_LengthGeneral_OrdGeneral";
195  return "pp_Mult_mm_Noether__FieldGeneral_LengthGeneral_OrdGeneral";
196  case p_Mult_mm_Proc:
197  return "p_Mult_mm__FieldGeneral_LengthGeneral_OrdGeneral";
199  return "p_Minus_mm_Mult_qq__FieldGeneral_LengthGeneral_OrdGeneral";
201  return "pp_Mult_Coeff_mm_DivSelect__FieldGeneral_LengthGeneral_OrdGeneral";
203  return "pp_Mult_Coeff_mm_DivSelectMult__FieldGeneral_LengthGeneral_OrdGeneral";
204 #endif
205  }
206  return "p_Unknown_Proc";
207 }
208 #endif
209 
210 
211 static void* GetDynamicProc(const char* proc_s, p_Proc proc,
212  p_Field field, p_Length length, p_Ord ord
213 #ifdef RDEBUG
214  , int get_name = 0
215 #endif
216  )
217 {
218  void* proc_ptr = NULL;
219  char proc_name[MAX_PROCNAME_LEN];
220  sprintf(proc_name, "%s__%s_%s_%s", proc_s,
221  p_FieldEnum_2_String(field),
222  p_LengthEnum_2_String(length),
223  p_OrdEnum_2_String(ord));
224  // printf("set %s\n",proc_name);
225  // first, try to get the proc from the kernel
226  proc_ptr = dynl_sym(DYNL_KERNEL_HANDLE, proc_name);
227  if (proc_ptr == NULL)
228  {
229  proc_ptr = dynl_sym_warn(p_ProcGetHandle(proc, field), proc_name, WARN_MSG);
230  // last but not least use general proc
231  if (proc_ptr == NULL)
232  {
233  proc_ptr = GetGeneralProc(proc);
234 #ifdef RDEBUG
235  sprintf(proc_name,"%s", GetGeneralProcName(proc));
236 #endif
237  }
238  }
239 #ifdef RDEBUG
240  if (get_name)
241  {
242  char* name = omStrDup(proc_name);
243 #if (!defined(SING_NDEBUG)) && (!defined(OM_NDEBUG))
244  omMarkAsStaticAddr(name);
245 #endif
246  return (void*) name;
247  }
248 #endif
249  return proc_ptr;
250 }
251 
252 
253 #define DoReallySetProc(what, field, length, ord) \
254  _p_procs->what = cast_vptr_to_A<what##_Proc_Ptr>( \
255  GetDynamicProc(#what, what##_Proc, field, length, ord))
256 
257 #ifdef RDEBUG
258 #define DoSetProc(what, field, length, ord) \
259 do \
260 { \
261  if (set_names) \
262  _p_procs->what = cast_vptr_to_A<what##_Proc_Ptr>( \
263  GetDynamicProc(#what, what##_Proc, field, length, ord, 1)); \
264  else \
265  DoReallySetProc(what, field, length, ord); \
266 } \
267 while(0)
268 #else
269 #define DoSetProc DoReallySetProc
270 #endif
271 
272 #include "templates/p_Procs_Set.h"
273 
274 #endif
275 
static void * p_ProcInitHandle(void **handle, const char *module)
unsigned char * proc[NUM_PROC]
Definition: checklibs.c:14
static void * p_procs_handle_FieldIndep
void * dynl_open_binary_warn(const char *binary_name, const char *msg)
Definition: mod_raw.cc:56
static void * p_procs_handle_FieldGeneral
static const char * p_OrdEnum_2_String(p_Ord ord)
Definition: p_Procs_Impl.h:237
#define TRUE
Definition: auxiliary.h:99
static void * GetDynamicProc(const char *proc_s, p_Proc proc, p_Field field, p_Length length, p_Ord ord, int get_name=0)
void * dynl_sym(void *handle, const char *symbol)
Definition: mod_raw.cc:167
void * dynl_sym_warn(void *handle, const char *proc, const char *msg)
Definition: mod_raw.cc:100
p_Length
Definition: p_Procs_Impl.h:122
p_Proc
Definition: p_Procs_Impl.h:178
#define WARN_MSG
#define RDEBUG
Definition: mod2.h:323
static void * p_ProcGetHandle(p_Proc proc, p_Field field)
static void * GetGeneralProc(p_Proc proc)
#define assume(x)
Definition: mod2.h:403
p_Field
Definition: p_Procs_Impl.h:106
static void * p_procs_handle_FieldZp
All the auxiliary stuff.
#define dReportBug(s)
Definition: reporter.h:112
char name(const Variable &v)
Definition: factory.h:178
static const char * p_LengthEnum_2_String(p_Length length)
Definition: p_Procs_Impl.h:219
p_Ord
Definition: p_Procs_Impl.h:135
static void * p_procs_handle_FieldQ
#define NULL
Definition: omList.c:10
void * cast_A_to_vptr(A a)
Definition: auxiliary.h:385
const BOOLEAN p_procs_dynamic
static const char * p_ProcField_2_Module(p_Proc proc, p_Field field)
Definition: p_Procs_Impl.h:331
static const char * GetGeneralProcName(p_Proc proc)
#define DYNL_KERNEL_HANDLE
Definition: mod_raw.h:32
void omMarkAsStaticAddr(void *addr)
#define MAX_PROCNAME_LEN
int BOOLEAN
Definition: auxiliary.h:86
static const char * p_FieldEnum_2_String(p_Field field)
Definition: p_Procs_Impl.h:198
#define omStrDup(s)
Definition: omAllocDecl.h:263