SuperLU Distributed 8.2.1
Distributed memory sparse direct solver
superlu_ddefs.h
Go to the documentation of this file.
1
28#ifndef __SUPERLU_DDEFS /* allow multiple inclusions */
29#define __SUPERLU_DDEFS
30
31/*
32 * File name: superlu_ddefs.h
33 * Purpose: Distributed SuperLU data types and function prototypes
34 * History:
35 */
36
37#include "superlu_defs.h"
38
39/*
40 *-- The structure used to store matrix A of the linear system and
41 * several vectors describing the transformations done to matrix A.
42 *
43 * A (SuperMatrix*)
44 * Matrix A in A*X=B, of dimension (A->nrow, A->ncol).
45 * The number of linear equations is A->nrow. The type of A can be:
46 * Stype = SLU_NC; Dtype = SLU_D; Mtype = SLU_GE.
47 *
48 * DiagScale (DiagScale_t)
49 * Specifies the form of equilibration that was done.
50 * = NOEQUIL: No equilibration.
51 * = ROW: Row equilibration, i.e., A was premultiplied by diag(R).
52 * = COL: Column equilibration, i.e., A was postmultiplied by diag(C).
53 * = BOTH: Both row and column equilibration, i.e., A was replaced
54 * by diag(R)*A*diag(C).
55 *
56 * R double*, dimension (A->nrow)
57 * The row scale factors for A.
58 * If DiagScale = ROW or BOTH, A is multiplied on the left by diag(R).
59 * If DiagScale = NOEQUIL or COL, R is not defined.
60 *
61 * C double*, dimension (A->ncol)
62 * The column scale factors for A.
63 * If DiagScale = COL or BOTH, A is multiplied on the right by diag(C).
64 * If DiagScale = NOEQUIL or ROW, C is not defined.
65 *
66 * perm_r (int*) dimension (A->nrow)
67 * Row permutation vector which defines the permutation matrix Pr,
68 * perm_r[i] = j means row i of A is in position j in Pr*A.
69 *
70 * perm_c (int*) dimension (A->ncol)
71 * Column permutation vector, which defines the
72 * permutation matrix Pc; perm_c[i] = j means column i of A is
73 * in position j in A*Pc.
74 *
75 */
76typedef struct {
77 DiagScale_t DiagScale;
78 double *R;
79 double *C;
80 int_t *perm_r;
81 int_t *perm_c;
82} dScalePermstruct_t;
83
84#if 0 // Sherry: move to superlu_defs.h
85/*-- Auxiliary data type used in PxGSTRS/PxGSTRS1. */
86typedef struct {
87 int_t lbnum; /* Row block number (local). */
88 int_t indpos; /* Starting position in Uindex[]. */
89} Ucb_indptr_t;
90#endif
91
92/*
93 * On each processor, the blocks in L are stored in compressed block
94 * column format, the blocks in U are stored in compressed block row format.
95 */
96#define MAX_LOOKAHEADS 50
97typedef struct {
98 int_t **Lrowind_bc_ptr; /* size ceil(NSUPERS/Pc);
99 free'd in distribution routinies */
100 int_t *Lrowind_bc_dat; /* size sum of sizes of Lrowind_bc_ptr[lk]) */
101 long int *Lrowind_bc_offset; /* size ceil(NSUPERS/Pc) */
102 long int Lrowind_bc_cnt;
103
104 double **Lnzval_bc_ptr; /* size ceil(NSUPERS/Pc);
105 free'd in distribution routinies */
106 double *Lnzval_bc_dat; /* size sum of sizes of Lnzval_bc_ptr[lk]) */
107 long int *Lnzval_bc_offset; /* size ceil(NSUPERS/Pc) */
108 long int Lnzval_bc_cnt;
109
110 double **Linv_bc_ptr; /* size ceil(NSUPERS/Pc);
111 free'd in distribution routinies */
112 double *Linv_bc_dat; /* size sum of sizes of Linv_bc_ptr[lk]) */
113 long int *Linv_bc_offset; /* size ceil(NSUPERS/Pc) */
114 long int Linv_bc_cnt;
115
116 int_t **Lindval_loc_bc_ptr; /* size ceil(NSUPERS/Pc);
117 pointers to locations in Lrowind_bc_ptr and Lnzval_bc_ptr;
118 free'd in distribution routinies */
119
120 int_t *Lindval_loc_bc_dat; /* size: sum of sizes of Lindval_loc_bc_ptr[lk]) */
121 long int *Lindval_loc_bc_offset; /* size ceil(NSUPERS/Pc) */
122 long int Lindval_loc_bc_cnt;
123
124 /* for new U format -> */
125 int_t **Ucolind_bc_ptr; /* size ceil(NSUPERS/Pc) */
126 int_t *Ucolind_bc_dat; /* size: sum of sizes of Ucolind_bc_ptr[lk]) */
127 int64_t *Ucolind_bc_offset; /* size ceil(NSUPERS/Pc) */
128 int64_t Ucolind_bc_cnt;
129
130 double **Unzval_bc_ptr; /* size ceil(NSUPERS/Pc) */
131 double *Unzval_bc_dat; /* size: sum of sizes of Unzval_bc_ptr[lk]) */
132 int64_t *Unzval_bc_offset; /* size ceil(NSUPERS/Pc) */
133 int64_t Unzval_bc_cnt;
134
135 int_t **Uindval_loc_bc_ptr; /* size ceil(NSUPERS/Pc) pointers to locations in Ucolind_bc_ptr and Unzval_bc_ptr */
136 int_t *Uindval_loc_bc_dat; /* size: sum of sizes of Uindval_loc_bc_ptr[lk]) */
137 int64_t *Uindval_loc_bc_offset; /* size ceil(NSUPERS/Pc) */
138 int64_t Uindval_loc_bc_cnt;
139 /* end for new U format <- */
140
141 int_t *Unnz; /* number of nonzeros per block column in U*/
142 int_t **Lrowind_bc_2_lsum; /* size ceil(NSUPERS/Pc) map indices of Lrowind_bc_ptr to indices of lsum */
143 double **Uinv_bc_ptr; /* size ceil(NSUPERS/Pc) */
144 double *Uinv_bc_dat; /* size sum of sizes of Linv_bc_ptr[lk]) */
145 long int *Uinv_bc_offset; /* size ceil(NSUPERS/Pc) */
146 long int Uinv_bc_cnt;
147
148 int_t **Ufstnz_br_ptr; /* size ceil(NSUPERS/Pr) */
149 int_t *Ufstnz_br_dat; /* size sum of sizes of Ufstnz_br_ptr[lk]) */
150 long int *Ufstnz_br_offset; /* size ceil(NSUPERS/Pr) */
151 long int Ufstnz_br_cnt;
152
153 double **Unzval_br_ptr; /* size ceil(NSUPERS/Pr) */
154 double *Unzval_br_dat; /* size sum of sizes of Unzval_br_ptr[lk]) */
155 long int *Unzval_br_offset; /* size ceil(NSUPERS/Pr) */
156 long int Unzval_br_cnt;
157
158 /*-- Data structures used for broadcast and reduction trees. --*/
159 C_Tree *LBtree_ptr; /* size ceil(NSUPERS/Pc) */
160 C_Tree *LRtree_ptr; /* size ceil(NSUPERS/Pr) */
161 C_Tree *UBtree_ptr; /* size ceil(NSUPERS/Pc) */
162 C_Tree *URtree_ptr; /* size ceil(NSUPERS/Pr) */
163#if 0
164 int_t *Lsub_buf; /* Buffer for the remote subscripts of L */
165 double *Lval_buf; /* Buffer for the remote nonzeros of L */
166 int_t *Usub_buf; /* Buffer for the remote subscripts of U */
167 double *Uval_buf; /* Buffer for the remote nonzeros of U */
168#endif
169 int_t *Lsub_buf_2[MAX_LOOKAHEADS]; /* Buffers for the remote subscripts of L*/
170 double *Lval_buf_2[MAX_LOOKAHEADS]; /* Buffers for the remote nonzeros of L */
171 int_t *Usub_buf_2[MAX_LOOKAHEADS]; /* Buffer for the remote subscripts of U */
172 double *Uval_buf_2[MAX_LOOKAHEADS]; /* Buffer for the remote nonzeros of U */
173 double *ujrow; /* used in panel factorization. */
174 int_t bufmax[NBUFFERS]; /* Maximum buffer size across all MPI ranks:
175 * 0 : maximum size of Lsub_buf[]
176 * 1 : maximum size of Lval_buf[]
177 * 2 : maximum size of Usub_buf[]
178 * 3 : maximum size of Uval_buf[]
179 * 4 : maximum size of tempv[LDA]
180 */
181
182 /*-- Record communication schedule for factorization. --*/
183 int *ToRecv; /* Recv from no one (0), left (1), and up (2).*/
184 int *ToSendD; /* Whether need to send down block row. */
185 int **ToSendR; /* List of processes to send right block col. */
186
187 /*-- Record communication schedule for forward/back solves. --*/
188 /* 1/15/22 Sherry: changed int_t to int type */
189 int *fmod; /* Modification count for L-solve */
190 int **fsendx_plist; /* Column process list to send down Xk */
191 int *frecv; /* Modifications to be recv'd in proc row */
192 int nfrecvx; /* Number of Xk I will receive in L-solve */
193 int nfsendx; /* Number of Xk I will send in L-solve */
194 int *bmod; /* Modification count for U-solve */
195 int **bsendx_plist; /* Column process list to send down Xk */
196 int *brecv; /* Modifications to be recv'd in proc row */
197 int nbrecvx; /* Number of Xk I will receive in U-solve */
198 int nbsendx; /* Number of Xk I will send in U-solve */
199 int *mod_bit; /* Flag contribution from each row blocks */
200
201 /*-- Auxiliary arrays used for forward/back solves. --*/
202 int_t *ilsum; /* Starting position of each supernode in lsum
203 (local) */
204 int_t ldalsum; /* LDA of lsum (local) */
205 int_t SolveMsgSent; /* Number of actual messages sent in LU-solve */
206 int_t SolveMsgVol; /* Volume of messages sent in the solve phase */
207
208
209 /*********************/
210 /* The following variables are used in the hybrid solver */
211
212 /*-- Counts to be used in U^{-T} triangular solve. -- */
213 int_t UT_SOLVE;
214 int_t L_SOLVE;
215 int_t FRECV;
216 int_t ut_ldalsum; /* LDA of lsum (local) */
217 int_t *ut_ilsum; /* ilsum in column-wise */
218 int_t *utmod; /* Modification count for Ut-solve. */
219 int_t **ut_sendx_plist; /* Row process list to send down Xk */
220 int_t *utrecv; /* Modifications to be recev'd in proc column. */
221 int_t n_utsendx; /* Number of Xk I will receive */
222 int_t n_utrecvx; /* Number of Xk I will send */
223 int_t n_utrecvmod;
224 int_t nroot;
225 int_t *ut_modbit;
226 int_t *Urbs;
227 Ucb_indptr_t **Ucb_indptr;/* Vertical linked list pointing to Uindex[] */
228 Ucb_indptr_t *Ucb_inddat;
229 long int *Ucb_indoffset;
230 long int Ucb_indcnt;
231
232 int_t **Ucb_valptr; /* Vertical linked list pointing to Unzval[] */
233 int_t *Ucb_valdat;
234 long int *Ucb_valoffset;
235 long int Ucb_valcnt;
236
237 /* some additional counters for L solve */
238 int_t n;
239 int_t nleaf;
240 int_t nfrecvmod;
241 int_t inv; /* whether the diagonal block is inverted*/
242
243#ifdef GPU_ACC
244 /* The following variables are used in GPU trisolve */
245
246 int_t *d_Lrowind_bc_dat;
247 long int *d_Lrowind_bc_offset;
248 double *d_Lnzval_bc_dat;
249 long int *d_Lnzval_bc_offset;
250 int_t *d_Ucolind_bc_dat;
251 int64_t *d_Ucolind_bc_offset;
252 double *d_Unzval_bc_dat;
253 long int *d_Unzval_bc_offset;
254
255 double *d_Linv_bc_dat ;
256 double *d_Uinv_bc_dat ;
257 long int *d_Linv_bc_offset ;
258 long int *d_Uinv_bc_offset ;
259 int_t *d_Lindval_loc_bc_dat ;
260 int64_t *d_Lindval_loc_bc_offset ;
261 int_t *d_Uindval_loc_bc_dat ;
262 int64_t *d_Uindval_loc_bc_offset ;
263
264 // long int *d_Lindval_loc_bc_offset ;
265 // int_t *d_Urbs;
266 // int_t *d_Ufstnz_br_dat;
267 // long int *d_Ufstnz_br_offset;
268 // double *d_Unzval_br_dat;
269 // long int *d_Unzval_br_offset;
270 // int_t *d_Ucb_valdat;
271 // long int *d_Ucb_valoffset;
272 // Ucb_indptr_t *d_Ucb_inddat;
273 // long int *d_Ucb_indoffset;
274
275 int_t *d_ilsum ;
276 int_t *d_xsup ;
277 C_Tree *d_LBtree_ptr ;
278 C_Tree *d_LRtree_ptr ;
279 C_Tree *d_UBtree_ptr ;
280 C_Tree *d_URtree_ptr ;
281#endif
282
283} dLocalLU_t;
284
285
286typedef struct {
287 int_t *etree;
288 Glu_persist_t *Glu_persist;
289 dLocalLU_t *Llu;
290 char dt;
291} dLUstruct_t;
292
293
294/*-- Data structure for communication during matrix-vector multiplication. */
295typedef struct {
296 int_t *extern_start;
297 int_t *ind_tosend; /* X indeices to be sent to other processes */
298 int_t *ind_torecv; /* X indeices to be received from other processes */
299 int_t *ptr_ind_tosend;/* Printers to ind_tosend[] (Size procs)
300 (also point to val_torecv) */
301 int_t *ptr_ind_torecv;/* Printers to ind_torecv[] (Size procs)
302 (also point to val_tosend) */
303 int *SendCounts; /* Numbers of X indices to be sent
304 (also numbers of X values to be received) */
305 int *RecvCounts; /* Numbers of X indices to be received
306 (also numbers of X values to be sent) */
307 void *val_tosend; /* X values to be sent to other processes */
308 void *val_torecv; /* X values to be received from other processes */
309 int_t TotalIndSend; /* Total number of indices to be sent
310 (also total number of values to be received) */
311 int_t TotalValSend; /* Total number of values to be sent.
312 (also total number of indices to be received) */
313} pdgsmv_comm_t;
314
315/*-- Data structure holding the information for the solution phase --*/
316typedef struct {
317 int_t *row_to_proc;
318 int_t *inv_perm_c;
319 int_t num_diag_procs, *diag_procs, *diag_len;
320 pdgsmv_comm_t *gsmv_comm; /* communication metadata for SpMV,
321 required by IterRefine. */
322 pxgstrs_comm_t *gstrs_comm; /* communication metadata for SpTRSV. */
323 int_t *A_colind_gsmv; /* After pdgsmv_init(), the global column
324 indices of A are translated into the relative
325 positions in the gathered x-vector.
326 This is re-used in repeated calls to pdgsmv() */
327 int_t *xrow_to_proc; /* used by PDSLin */
328 NRformat_loc3d* A3d; /* Point to 3D {A, B} gathered on 2D layer 0.
329 This needs to be peresistent between
330 3D factorization and solve. */
331} dSOLVEstruct_t;
332
333
334
335/*==== For 3D code ====*/
336
337// new structures for pdgstrf_4_8
338
339#if 0 // Sherry: moved to superlu_defs.h
340typedef struct
341{
342 int_t nub;
343 int_t klst;
344 int_t ldu;
345 int_t* usub;
346 double* uval;
347} uPanelInfo_t;
348
349typedef struct
350{
351 int_t *lsub;
352 double *lusup;
353 int_t luptr0;
354 int_t nlb; //number of l blocks
355 int_t nsupr;
356} lPanelInfo_t;
357
358/* HyP_t is the data structure to assist HALO offload of Schur-complement. */
359typedef struct
360{
361 Remain_info_t *lookAhead_info, *Remain_info;
362 Ublock_info_t *Ublock_info, *Ublock_info_Phi;
363
364 int_t first_l_block_acc , first_u_block_acc;
365 int_t last_offload ;
366 int_t *Lblock_dirty_bit, * Ublock_dirty_bit;
367 double *lookAhead_L_buff, *Remain_L_buff;
368 int_t lookAheadBlk; /* number of blocks in look-ahead window */
369 int_t RemainBlk ; /* number of blocks outside look-ahead window */
370 int_t num_look_aheads, nsupers;
371 int_t ldu, ldu_Phi;
372 int_t num_u_blks, num_u_blks_Phi;
373
374 int_t jj_cpu;
375 double *bigU_Phi;
376 double *bigU_host;
377 int_t Lnbrow;
378 int_t Rnbrow;
379
380 int_t buffer_size;
381 int_t bigu_size;
382 int offloadCondition;
383 int superlu_acc_offload;
384 int nGPUStreams;
385} HyP_t;
386
387#endif // Above are moved to superlu_defs.h
388
389
390typedef struct
391{
392 int_t * Lsub_buf ;
393 double * Lval_buf ;
394 int_t * Usub_buf ;
395 double * Uval_buf ;
396} dLUValSubBuf_t;
397
398int_t scuStatUpdate(
399 int_t knsupc,
400 HyP_t* HyP,
401 SCT_t* SCT,
402 SuperLUStat_t *stat
403 );
404
405typedef struct
406{
407 gEtreeInfo_t gEtreeInfo;
408 int_t* iperm_c_supno;
409 int_t* myNodeCount;
410 int_t* myTreeIdxs;
411 int_t* myZeroTrIdxs;
412 int_t** treePerm;
413 sForest_t** sForests;
414 int_t* supernode2treeMap;
415 dLUValSubBuf_t *LUvsb;
416} dtrf3Dpartition_t;
417
418typedef struct
419{
420 double *bigU;
421 double *bigV;
422} dscuBufs_t;
423
424typedef struct
425{
426 double* BlockLFactor;
427 double* BlockUFactor;
428} ddiagFactBufs_t;
429
430/*=====================*/
431
432/***********************************************************************
433 * Function prototypes
434 ***********************************************************************/
435
436#ifdef __cplusplus
437extern "C" {
438#endif
439
440
441/* Supernodal LU factor related */
442extern void
443dCreate_CompCol_Matrix_dist(SuperMatrix *, int_t, int_t, int_t, double *,
444 int_t *, int_t *, Stype_t, Dtype_t, Mtype_t);
445extern void
446dCreate_CompRowLoc_Matrix_dist(SuperMatrix *, int_t, int_t, int_t, int_t,
447 int_t, double *, int_t *, int_t *,
448 Stype_t, Dtype_t, Mtype_t);
449extern void
450dCompRow_to_CompCol_dist(int_t, int_t, int_t, double *, int_t *, int_t *,
451 double **, int_t **, int_t **);
452extern int
453pdCompRow_loc_to_CompCol_global(int_t, SuperMatrix *, gridinfo_t *,
454 SuperMatrix *);
455extern void
456dCopy_CompCol_Matrix_dist(SuperMatrix *, SuperMatrix *);
457extern void
458dCreate_Dense_Matrix_dist(SuperMatrix *, int_t, int_t, double *, int_t,
459 Stype_t, Dtype_t, Mtype_t);
460extern void
461dCreate_SuperNode_Matrix_dist(SuperMatrix *, int_t, int_t, int_t, double *,
462 int_t *, int_t *, int_t *, int_t *, int_t *,
463 Stype_t, Dtype_t, Mtype_t);
464extern void
465dCopy_Dense_Matrix_dist(int_t, int_t, double *, int_t,
466 double *, int_t);
467
468extern void dallocateA_dist (int_t, int_t, double **, int_t **, int_t **);
469extern void dGenXtrue_dist (int_t, int_t, double *, int_t);
470extern void dFillRHS_dist (char *, int_t, double *, int_t,
471 SuperMatrix *, double *, int_t);
472extern int dcreate_matrix(SuperMatrix *, int, double **, int *,
473 double **, int *, FILE *, gridinfo_t *);
474extern int dcreate_matrix_rb(SuperMatrix *, int, double **, int *,
475 double **, int *, FILE *, gridinfo_t *);
476extern int dcreate_matrix_dat(SuperMatrix *, int, double **, int *,
477 double **, int *, FILE *, gridinfo_t *);
478extern int dcreate_matrix_postfix(SuperMatrix *, int, double **, int *,
479 double **, int *, FILE *, char *, gridinfo_t *);
480
481extern void dScalePermstructInit(const int_t, const int_t,
482 dScalePermstruct_t *);
483extern void dScalePermstructFree(dScalePermstruct_t *);
484
485/* Driver related */
486extern void dgsequ_dist (SuperMatrix *, double *, double *, double *,
487 double *, double *, int_t *);
488extern double dlangs_dist (char *, SuperMatrix *);
489extern void dlaqgs_dist (SuperMatrix *, double *, double *, double,
490 double, double, char *);
491extern void pdgsequ (SuperMatrix *, double *, double *, double *,
492 double *, double *, int_t *, gridinfo_t *);
493extern double pdlangs (char *, SuperMatrix *, gridinfo_t *);
494extern void pdlaqgs (SuperMatrix *, double *, double *, double,
495 double, double, char *);
496extern int pdPermute_Dense_Matrix(int_t, int_t, int_t [], int_t[],
497 double [], int, double [], int, int,
498 gridinfo_t *);
499
500extern int sp_dtrsv_dist (char *, char *, char *, SuperMatrix *,
501 SuperMatrix *, double *, int *);
502extern int sp_dgemv_dist (char *, double, SuperMatrix *, double *,
503 int, double, double *, int);
504extern int sp_dgemm_dist (char *, int, double, SuperMatrix *,
505 double *, int, double, double *, int);
506
507extern float ddistribute(superlu_dist_options_t *,
508 int_t, SuperMatrix *, Glu_freeable_t *,
509 dLUstruct_t *, gridinfo_t *);
510extern void pdgssvx_ABglobal(superlu_dist_options_t *, SuperMatrix *,
511 dScalePermstruct_t *, double *,
512 int, int, gridinfo_t *, dLUstruct_t *, double *,
513 SuperLUStat_t *, int *);
514extern float pddistribute(superlu_dist_options_t *, int_t, SuperMatrix *,
515 dScalePermstruct_t *, Glu_freeable_t *,
516 dLUstruct_t *, gridinfo_t *);
517extern void pdgssvx(superlu_dist_options_t *, SuperMatrix *,
518 dScalePermstruct_t *, double *,
519 int, int, gridinfo_t *, dLUstruct_t *,
520 dSOLVEstruct_t *, double *, SuperLUStat_t *, int *);
521extern void pdCompute_Diag_Inv(int_t, dLUstruct_t *,gridinfo_t *, SuperLUStat_t *, int *);
522extern int dSolveInit(superlu_dist_options_t *, SuperMatrix *, int_t [], int_t [],
523 int_t, dLUstruct_t *, gridinfo_t *, dSOLVEstruct_t *);
524extern void dSolveFinalize(superlu_dist_options_t *, dSOLVEstruct_t *);
525extern void dDestroy_A3d_gathered_on_2d(dSOLVEstruct_t *, gridinfo3d_t *);
526extern int_t pdgstrs_init(int_t, int_t, int_t, int_t,
527 int_t [], int_t [], gridinfo_t *grid,
528 Glu_persist_t *, dSOLVEstruct_t *);
529extern void pxgstrs_finalize(pxgstrs_comm_t *);
530extern int dldperm_dist(int, int, int_t, int_t [], int_t [],
531 double [], int_t *, double [], double []);
532extern int dstatic_schedule(superlu_dist_options_t *, int, int,
533 dLUstruct_t *, gridinfo_t *, SuperLUStat_t *,
534 int_t *, int_t *, int *);
535extern void dLUstructInit(const int_t, dLUstruct_t *);
536extern void dLUstructFree(dLUstruct_t *);
537extern void dDestroy_LU(int_t, gridinfo_t *, dLUstruct_t *);
538extern void dDestroy_Tree(int_t, gridinfo_t *, dLUstruct_t *);
539extern void dscatter_l (int ib, int ljb, int nsupc, int_t iukp, int_t* xsup,
540 int klst, int nbrow, int_t lptr, int temp_nbrow,
541 int_t* usub, int_t* lsub, double *tempv,
542 int* indirect_thread, int* indirect2,
543 int_t ** Lrowind_bc_ptr, double **Lnzval_bc_ptr,
544 gridinfo_t * grid);
545extern void dscatter_u (int ib, int jb, int nsupc, int_t iukp, int_t * xsup,
546 int klst, int nbrow, int_t lptr, int temp_nbrow,
547 int_t* lsub, int_t* usub, double* tempv,
548 int_t ** Ufstnz_br_ptr, double **Unzval_br_ptr,
549 gridinfo_t * grid);
550extern int_t pdgstrf(superlu_dist_options_t *, int, int, double anorm,
551 dLUstruct_t*, gridinfo_t*, SuperLUStat_t*, int*);
552
553/* #define GPU_PROF
554#define IPM_PROF */
555
556/* Solve related */
557extern void pdgstrs_Bglobal(superlu_dist_options_t *,
558 int_t, dLUstruct_t *, gridinfo_t *,
559 double *, int_t, int, SuperLUStat_t *, int *);
560extern void pdgstrs(superlu_dist_options_t *, int_t,
561 dLUstruct_t *, dScalePermstruct_t *, gridinfo_t *,
562 double *, int_t, int_t, int_t, int, dSOLVEstruct_t *,
563 SuperLUStat_t *, int *);
564extern void pdgstrf2_trsm(superlu_dist_options_t * options, int_t k0, int_t k,
565 double thresh, Glu_persist_t *, gridinfo_t *,
566 dLocalLU_t *, MPI_Request *, int tag_ub,
567 SuperLUStat_t *, int *info);
568extern void pdgstrs2_omp(int_t k0, int_t k, Glu_persist_t *, gridinfo_t *,
569 dLocalLU_t *, Ublock_info_t *, SuperLUStat_t *);
570extern int_t pdReDistribute_B_to_X(double *B, int_t m_loc, int nrhs, int_t ldb,
571 int_t fst_row, int_t *ilsum, double *x,
572 dScalePermstruct_t *, Glu_persist_t *,
573 gridinfo_t *, dSOLVEstruct_t *);
574extern void dlsum_fmod(double *, double *, double *, double *,
575 int, int, int_t , int *fmod, int_t, int_t, int_t,
576 int_t *, gridinfo_t *, dLocalLU_t *,
577 MPI_Request [], SuperLUStat_t *);
578extern void dlsum_bmod(double *, double *, double *,
579 int, int_t, int *bmod, int_t *, Ucb_indptr_t **,
580 int_t **, int_t *, gridinfo_t *, dLocalLU_t *,
581 MPI_Request [], SuperLUStat_t *);
582
583extern void dlsum_fmod_inv(double *, double *, double *, double *,
584 int, int_t , int *fmod,
585 int_t *, gridinfo_t *, dLocalLU_t *,
586 SuperLUStat_t **, int_t *, int_t *, int_t, int_t, int_t, int_t, int, int);
587extern void dlsum_fmod_inv_master(double *, double *, double *, double *,
588 int, int, int_t , int *fmod, int_t,
589 int_t *, gridinfo_t *, dLocalLU_t *,
590 SuperLUStat_t **, int_t, int_t, int_t, int_t, int, int);
591extern void dlsum_bmod_inv(double *, double *, double *, double *,
592 int, int_t, int *bmod, int_t *, Ucb_indptr_t **,
593 int_t **, int_t *, gridinfo_t *, dLocalLU_t *,
594 SuperLUStat_t **, int_t *, int_t *, int_t, int_t, int, int);
595extern void dlsum_bmod_inv_master(double *, double *, double *, double *,
596 int, int_t, int *bmod, int_t *, Ucb_indptr_t **,
597 int_t **, int_t *, gridinfo_t *, dLocalLU_t *,
598 SuperLUStat_t **, int_t, int_t, int, int);
599
600extern void dComputeLevelsets(int , int_t , gridinfo_t *,
601 Glu_persist_t *, dLocalLU_t *, int_t *);
602
603#ifdef GPU_ACC
604extern void pdconvertU(superlu_dist_options_t *, gridinfo_t *, dLUstruct_t *, SuperLUStat_t *, int);
605extern void dlsum_fmod_inv_gpu_wrap(int_t, int_t, int_t, int_t, double *, double *, int, int, int_t , int *fmod, C_Tree *, C_Tree *, int_t *, int_t *, int64_t *, double *, int64_t *, double *, int64_t *, int_t *, int64_t *, int_t *, gridinfo_t *, double * , double * , int_t );
606extern void dlsum_bmod_inv_gpu_wrap(superlu_dist_options_t *, int_t, int_t, int_t, int_t, double *, double *, int, int, int_t , int *bmod, C_Tree *, C_Tree *, int_t *, int_t *, int64_t *, double *, int64_t *, double *, int64_t *, int_t *, int64_t *,int_t *,gridinfo_t *);
607#endif
608
609extern void pdgsrfs(superlu_dist_options_t *, int_t,
610 SuperMatrix *, double, dLUstruct_t *,
611 dScalePermstruct_t *, gridinfo_t *,
612 double [], int_t, double [], int_t, int,
613 dSOLVEstruct_t *, double *, SuperLUStat_t *, int *);
614extern void pdgsrfs_ABXglobal(superlu_dist_options_t *, int_t,
615 SuperMatrix *, double, dLUstruct_t *,
616 gridinfo_t *, double *, int_t, double *, int_t,
617 int, double *, SuperLUStat_t *, int *);
618extern int pdgsmv_AXglobal_setup(SuperMatrix *, Glu_persist_t *,
619 gridinfo_t *, int_t *, int_t *[],
620 double *[], int_t *[], int_t []);
621extern int pdgsmv_AXglobal(int_t, int_t [], double [], int_t [],
622 double [], double []);
623extern int pdgsmv_AXglobal_abs(int_t, int_t [], double [], int_t [],
624 double [], double []);
625extern void pdgsmv_init(SuperMatrix *, int_t *, gridinfo_t *,
626 pdgsmv_comm_t *);
627extern void pdgsmv(int_t, SuperMatrix *, gridinfo_t *, pdgsmv_comm_t *,
628 double x[], double ax[]);
629extern void pdgsmv_finalize(pdgsmv_comm_t *);
630
631/* Memory-related */
632extern double *doubleMalloc_dist(int_t);
633extern double *doubleCalloc_dist(int_t);
634extern void *duser_malloc_dist (int_t, int_t);
635extern void duser_free_dist (int_t, int_t);
636extern int_t dQuerySpace_dist(int_t, dLUstruct_t *, gridinfo_t *,
637 SuperLUStat_t *, superlu_dist_mem_usage_t *);
638
639/* Auxiliary routines */
640
641extern void dClone_CompRowLoc_Matrix_dist(SuperMatrix *, SuperMatrix *);
642extern void dCopy_CompRowLoc_Matrix_dist(SuperMatrix *, SuperMatrix *);
643extern void dZero_CompRowLoc_Matrix_dist(SuperMatrix *);
644extern void dScaleAddId_CompRowLoc_Matrix_dist(SuperMatrix *, double);
645extern void dScaleAdd_CompRowLoc_Matrix_dist(SuperMatrix *, SuperMatrix *, double);
646extern void dZeroLblocks(int, int, gridinfo_t *, dLUstruct_t *);
647extern void dZeroUblocks(int iam, int n, gridinfo_t *, dLUstruct_t *);
648extern void dfill_dist (double *, int_t, double);
649extern void dinf_norm_error_dist (int_t, int_t, double*, int_t,
650 double*, int_t, gridinfo_t*);
651extern void pdinf_norm_error(int, int_t, int_t, double [], int_t,
652 double [], int_t , MPI_Comm);
653extern void dreadhb_dist (int, FILE *, int_t *, int_t *, int_t *,
654 double **, int_t **, int_t **);
655extern void dreadtriple_dist(FILE *, int_t *, int_t *, int_t *,
656 double **, int_t **, int_t **);
657extern void dreadtriple_noheader(FILE *, int_t *, int_t *, int_t *,
658 double **, int_t **, int_t **);
659extern void dreadrb_dist(int, FILE *, int_t *, int_t *, int_t *,
660 double **, int_t **, int_t **);
661extern void dreadMM_dist(FILE *, int_t *, int_t *, int_t *,
662 double **, int_t **, int_t **);
663extern int dread_binary(FILE *, int_t *, int_t *, int_t *,
664 double **, int_t **, int_t **);
665
666/* Distribute the data for numerical factorization */
667extern float ddist_psymbtonum(superlu_dist_options_t *, int_t, SuperMatrix *,
668 dScalePermstruct_t *, Pslu_freeable_t *,
669 dLUstruct_t *, gridinfo_t *);
670extern void pdGetDiagU(int_t, dLUstruct_t *, gridinfo_t *, double *);
671
672extern int d_c2cpp_GetHWPM(SuperMatrix *, gridinfo_t *, dScalePermstruct_t *);
673
674/* Routines for debugging */
675extern void dPrintLblocks(int, int_t, gridinfo_t *, Glu_persist_t *,
676 dLocalLU_t *);
677extern void dPrintUblocks(int, int_t, gridinfo_t *, Glu_persist_t *,
678 dLocalLU_t *);
679extern void dPrint_CompCol_Matrix_dist(SuperMatrix *);
680extern void dPrint_Dense_Matrix_dist(SuperMatrix *);
681extern int dPrint_CompRowLoc_Matrix_dist(SuperMatrix *);
682extern int file_dPrint_CompRowLoc_Matrix_dist(FILE *fp, SuperMatrix *A);
683extern void Printdouble5(char *, int_t, double *);
684extern int file_Printdouble5(FILE *, char *, int_t, double *);
685
686extern void dGenCOOLblocks(int, int_t, gridinfo_t*,
687 Glu_persist_t*, dLocalLU_t *, int_t** , int_t** , double ** , int_t* , int_t* );
688extern void dGenCSCLblocks(int, int_t, gridinfo_t*,
689 Glu_persist_t*, dLocalLU_t *, double **, int_t **, int_t **, int_t*, int_t*);
690extern void dGenCSRLblocks(int, int_t, gridinfo_t*,
691 Glu_persist_t*, dLocalLU_t *, double **, int_t **, int_t **, int_t*, int_t*);
692
693
694/* BLAS */
695
696#ifdef USE_VENDOR_BLAS
697extern void dgemm_(const char*, const char*, const int*, const int*, const int*,
698 const double*, const double*, const int*, const double*,
699 const int*, const double*, double*, const int*, int, int);
700extern void dtrsv_(char*, char*, char*, int*, double*, int*,
701 double*, int*, int, int, int);
702extern void dtrsm_(const char*, const char*, const char*, const char*,
703 const int*, const int*, const double*, const double*, const int*,
704 double*, const int*, int, int, int, int);
705extern void dgemv_(const char *, const int *, const int *, const double *,
706 const double *a, const int *, const double *, const int *,
707 const double *, double *, const int *, int);
708
709#else
710extern int dgemm_(const char*, const char*, const int*, const int*, const int*,
711 const double*, const double*, const int*, const double*,
712 const int*, const double*, double*, const int*);
713extern int dtrsv_(char*, char*, char*, int*, double*, int*,
714 double*, int*);
715extern int dtrsm_(const char*, const char*, const char*, const char*,
716 const int*, const int*, const double*, const double*, const int*,
717 double*, const int*);
718extern void dgemv_(const char *, const int *, const int *, const double *,
719 const double *a, const int *, const double *, const int *,
720 const double *, double *, const int *);
721#endif
722
723extern void dger_(const int*, const int*, const double*,
724 const double*, const int*, const double*, const int*,
725 double*, const int*);
726
727extern int dscal_(const int *n, const double *alpha, double *dx, const int *incx);
728extern int daxpy_(const int *n, const double *alpha, const double *x,
729 const int *incx, double *y, const int *incy);
730
731/* SuperLU BLAS interface: dsuperlu_blas.c */
732extern int superlu_dgemm(const char *transa, const char *transb,
733 int m, int n, int k, double alpha, double *a,
734 int lda, double *b, int ldb, double beta, double *c, int ldc);
735extern int superlu_dtrsm(const char *sideRL, const char *uplo,
736 const char *transa, const char *diag, const int m, const int n,
737 const double alpha, const double *a,
738 const int lda, double *b, const int ldb);
739extern int superlu_dger(const int m, const int n, const double alpha,
740 const double *x, const int incx, const double *y,
741 const int incy, double *a, const int lda);
742extern int superlu_dscal(const int n, const double alpha, double *x, const int incx);
743extern int superlu_daxpy(const int n, const double alpha,
744 const double *x, const int incx, double *y, const int incy);
745extern int superlu_dgemv(const char *trans, const int m,
746 const int n, const double alpha, const double *a,
747 const int lda, const double *x, const int incx,
748 const double beta, double *y, const int incy);
749extern int superlu_dtrsv(char *uplo, char *trans, char *diag,
750 int n, double *a, int lda, double *x, int incx);
751
752#ifdef SLU_HAVE_LAPACK
753extern void dtrtri_(char*, char*, int*, double*, int*, int*);
754#endif
755
756/*==== For 3D code ====*/
757extern int dcreate_matrix3d(SuperMatrix *A, int nrhs, double **rhs,
758 int *ldb, double **x, int *ldx,
759 FILE *fp, gridinfo3d_t *grid3d);
760extern int dcreate_matrix_postfix3d(SuperMatrix *A, int nrhs, double **rhs,
761 int *ldb, double **x, int *ldx,
762 FILE *fp, char * postfix, gridinfo3d_t *grid3d);
763
764/* Matrix distributed in NRformat_loc in 3D process grid. It converts
765 it to a NRformat_loc distributed in 2D grid in grid-0 */
766extern void dGatherNRformat_loc3d(fact_t Fact, NRformat_loc *A, double *B,
767 int ldb, int nrhs, gridinfo3d_t *grid3d,
768 NRformat_loc3d **);
769extern int dScatter_B3d(NRformat_loc3d *A3d, gridinfo3d_t *grid3d);
770
771extern void pdgssvx3d (superlu_dist_options_t *, SuperMatrix *,
772 dScalePermstruct_t *, double B[], int ldb, int nrhs,
773 gridinfo3d_t *, dLUstruct_t *, dSOLVEstruct_t *,
774 double *berr, SuperLUStat_t *, int *info);
775extern int_t pdgstrf3d(superlu_dist_options_t *, int m, int n, double anorm,
776 dtrf3Dpartition_t*, SCT_t *, dLUstruct_t *,
777 gridinfo3d_t *, SuperLUStat_t *, int *);
778extern void dInit_HyP(HyP_t* HyP, dLocalLU_t *Llu, int_t mcb, int_t mrb );
779extern void Free_HyP(HyP_t* HyP);
780extern int updateDirtyBit(int_t k0, HyP_t* HyP, gridinfo_t* grid);
781
782 /* from scatter.h */
783extern void
784dblock_gemm_scatter( int_t lb, int_t j, Ublock_info_t *Ublock_info,
785 Remain_info_t *Remain_info, double *L_mat, int ldl,
786 double *U_mat, int ldu, double *bigV,
787 // int_t jj0,
788 int_t knsupc, int_t klst,
789 int_t *lsub, int_t *usub, int_t ldt,
790 int_t thread_id,
791 int *indirect, int *indirect2,
792 int_t **Lrowind_bc_ptr, double **Lnzval_bc_ptr,
793 int_t **Ufstnz_br_ptr, double **Unzval_br_ptr,
794 int_t *xsup, gridinfo_t *, SuperLUStat_t *
795#ifdef SCATTER_PROFILE
796 , double *Host_TheadScatterMOP, double *Host_TheadScatterTimer
797#endif
798 );
799
800#ifdef _OPENMP
801/*this version uses a lock to prevent multiple thread updating the same block*/
802extern void
803dblock_gemm_scatter_lock( int_t lb, int_t j, omp_lock_t* lock,
804 Ublock_info_t *Ublock_info, Remain_info_t *Remain_info,
805 double *L_mat, int_t ldl, double *U_mat, int_t ldu,
806 double *bigV,
807 // int_t jj0,
808 int_t knsupc, int_t klst,
809 int_t *lsub, int_t *usub, int_t ldt,
810 int_t thread_id,
811 int *indirect, int *indirect2,
812 int_t **Lrowind_bc_ptr, double **Lnzval_bc_ptr,
813 int_t **Ufstnz_br_ptr, double **Unzval_br_ptr,
814 int_t *xsup, gridinfo_t *
815#ifdef SCATTER_PROFILE
816 , double *Host_TheadScatterMOP, double *Host_TheadScatterTimer
817#endif
818 );
819#endif
820
821extern int_t
822dblock_gemm_scatterTopLeft( int_t lb, int_t j, double* bigV,
823 int_t knsupc, int_t klst, int_t* lsub,
824 int_t * usub, int_t ldt,
825 int* indirect, int* indirect2,
826 HyP_t* HyP, dLUstruct_t *, gridinfo_t*,
827 SCT_t*SCT, SuperLUStat_t *
828 );
829extern int_t
830dblock_gemm_scatterTopRight( int_t lb, int_t j, double* bigV,
831 int_t knsupc, int_t klst, int_t* lsub,
832 int_t * usub, int_t ldt,
833 int* indirect, int* indirect2,
834 HyP_t* HyP, dLUstruct_t *, gridinfo_t*,
835 SCT_t*SCT, SuperLUStat_t * );
836extern int_t
837dblock_gemm_scatterBottomLeft( int_t lb, int_t j, double* bigV,
838 int_t knsupc, int_t klst, int_t* lsub,
839 int_t * usub, int_t ldt,
840 int* indirect, int* indirect2,
841 HyP_t* HyP, dLUstruct_t *, gridinfo_t*,
842 SCT_t*SCT, SuperLUStat_t * );
843extern int_t
844dblock_gemm_scatterBottomRight( int_t lb, int_t j, double* bigV,
845 int_t knsupc, int_t klst, int_t* lsub,
846 int_t * usub, int_t ldt,
847 int* indirect, int* indirect2,
848 HyP_t* HyP, dLUstruct_t *, gridinfo_t*,
849 SCT_t*SCT, SuperLUStat_t * );
850
851 /* from gather.h */
852extern void dgather_u(int_t num_u_blks,
853 Ublock_info_t *Ublock_info, int_t * usub,
854 double *uval, double *bigU, int_t ldu,
855 int_t *xsup, int_t klst /* for SuperSize */
856 );
857
858extern void dgather_l( int_t num_LBlk, int_t knsupc,
859 Remain_info_t *L_info,
860 double * lval, int_t LD_lval,
861 double * L_buff );
862
863extern void dRgather_L(int_t k, int_t *lsub, double *lusup, gEtreeInfo_t*,
864 Glu_persist_t *, gridinfo_t *, HyP_t *,
865 int_t *myIperm, int_t *iperm_c_supno );
866extern void dRgather_U(int_t k, int_t jj0, int_t *usub, double *uval,
867 double *bigU, gEtreeInfo_t*, Glu_persist_t *,
868 gridinfo_t *, HyP_t *, int_t *myIperm,
869 int_t *iperm_c_supno, int_t *perm_u);
870
871 /* from xtrf3Dpartition.h */
872extern dtrf3Dpartition_t* dinitTrf3Dpartition(int_t nsupers,
873 superlu_dist_options_t *options,
874 dLUstruct_t *LUstruct, gridinfo3d_t * grid3d);
875extern void dDestroy_trf3Dpartition(dtrf3Dpartition_t *trf3Dpartition, gridinfo3d_t *grid3d);
876
877extern void d3D_printMemUse(dtrf3Dpartition_t* trf3Dpartition,
878 dLUstruct_t *LUstruct, gridinfo3d_t * grid3d);
879
880//extern int* getLastDep(gridinfo_t *grid, SuperLUStat_t *stat,
881// superlu_dist_options_t *options, dLocalLU_t *Llu,
882// int_t* xsup, int_t num_look_aheads, int_t nsupers,
883// int_t * iperm_c_supno);
884
885extern void dinit3DLUstructForest( int_t* myTreeIdxs, int_t* myZeroTrIdxs,
886 sForest_t** sForests, dLUstruct_t* LUstruct,
887 gridinfo3d_t* grid3d);
888
889extern int_t dgatherAllFactoredLUFr(int_t* myZeroTrIdxs, sForest_t* sForests,
890 dLUstruct_t* LUstruct, gridinfo3d_t* grid3d,
891 SCT_t* SCT );
892
893 /* The following are from pdgstrf2.h */
894extern int_t dLpanelUpdate(int_t off0, int_t nsupc, double* ublk_ptr,
895 int_t ld_ujrow, double* lusup, int_t nsupr, SCT_t*);
896extern void Local_Dgstrf2(superlu_dist_options_t *options, int_t k,
897 double thresh, double *BlockUFactor, Glu_persist_t *,
898 gridinfo_t *, dLocalLU_t *,
899 SuperLUStat_t *, int *info, SCT_t*);
900extern int_t dTrs2_GatherU(int_t iukp, int_t rukp, int_t klst,
901 int_t nsupc, int_t ldu, int_t *usub,
902 double* uval, double *tempv);
903extern int_t dTrs2_ScatterU(int_t iukp, int_t rukp, int_t klst,
904 int_t nsupc, int_t ldu, int_t *usub,
905 double* uval, double *tempv);
906extern int_t dTrs2_GatherTrsmScatter(int_t klst, int_t iukp, int_t rukp,
907 int_t *usub, double* uval, double *tempv,
908 int_t knsupc, int nsupr, double* lusup,
909 Glu_persist_t *Glu_persist) ;
910extern void pdgstrs2
911#ifdef _CRAY
912(
913 int_t m, int_t k0, int_t k, Glu_persist_t *Glu_persist, gridinfo_t *grid,
914 dLocalLU_t *Llu, SuperLUStat_t *stat, _fcd ftcs1, _fcd ftcs2, _fcd ftcs3
915);
916#else
917(
918 int_t m, int_t k0, int_t k, Glu_persist_t *Glu_persist, gridinfo_t *grid,
919 dLocalLU_t *Llu, SuperLUStat_t *stat
920);
921#endif
922
923extern void pdgstrf2(superlu_dist_options_t *, int_t nsupers, int_t k0,
924 int_t k, double thresh, Glu_persist_t *, gridinfo_t *,
925 dLocalLU_t *, MPI_Request *, int, SuperLUStat_t *, int *);
926
927 /* from p3dcomm.h */
928extern int_t dAllocLlu_3d(int_t nsupers, dLUstruct_t * LUstruct, gridinfo3d_t* grid3d);
929extern int_t dp3dScatter(int_t n, dLUstruct_t * LUstruct, gridinfo3d_t* grid3d);
930extern int_t dscatter3dLPanels(int_t nsupers,
931 dLUstruct_t * LUstruct, gridinfo3d_t* grid3d);
932extern int_t dscatter3dUPanels(int_t nsupers,
933 dLUstruct_t * LUstruct, gridinfo3d_t* grid3d);
934extern int_t dcollect3dLpanels(int_t layer, int_t nsupers, dLUstruct_t * LUstruct, gridinfo3d_t* grid3d);
935extern int_t dcollect3dUpanels(int_t layer, int_t nsupers, dLUstruct_t * LUstruct, gridinfo3d_t* grid3d);
936extern int_t dp3dCollect(int_t layer, int_t n, dLUstruct_t * LUstruct, gridinfo3d_t* grid3d);
937/*zero out LU non zero entries*/
938extern int_t dzeroSetLU(int_t nnodes, int_t* nodeList , dLUstruct_t *, gridinfo3d_t*);
939extern int dAllocGlu_3d(int_t n, int_t nsupers, dLUstruct_t *);
940extern int dDeAllocLlu_3d(int_t n, dLUstruct_t *, gridinfo3d_t*);
941extern int dDeAllocGlu_3d(dLUstruct_t *);
942
943/* Reduces L and U panels of nodes in the List nodeList (size=nnnodes)
944receiver[L(nodelist)] =sender[L(nodelist)] +receiver[L(nodelist)]
945receiver[U(nodelist)] =sender[U(nodelist)] +receiver[U(nodelist)]
946*/
947int_t dreduceAncestors3d(int_t sender, int_t receiver,
948 int_t nnodes, int_t* nodeList,
949 double* Lval_buf, double* Uval_buf,
950 dLUstruct_t* LUstruct, gridinfo3d_t* grid3d, SCT_t* SCT);
951/*reduces all nodelists required in a level*/
952extern int dreduceAllAncestors3d(int_t ilvl, int_t* myNodeCount,
953 int_t** treePerm,
954 dLUValSubBuf_t* LUvsb,
955 dLUstruct_t* LUstruct,
956 gridinfo3d_t* grid3d,
957 SCT_t* SCT );
958/*
959 Copies factored L and U panels from sender grid to receiver grid
960 receiver[L(nodelist)] <-- sender[L(nodelist)];
961 receiver[U(nodelist)] <-- sender[U(nodelist)];
962*/
963int_t dgatherFactoredLU(int_t sender, int_t receiver,
964 int_t nnodes, int_t *nodeList, dLUValSubBuf_t* LUvsb,
965 dLUstruct_t* LUstruct, gridinfo3d_t* grid3d,SCT_t* SCT );
966
967/*Gathers all the L and U factors to grid 0 for solve stage
968 By repeatidly calling above function*/
969int_t dgatherAllFactoredLU(dtrf3Dpartition_t* trf3Dpartition, dLUstruct_t* LUstruct,
970 gridinfo3d_t* grid3d, SCT_t* SCT );
971
972/*Distributes data in each layer and initilizes ancestors
973 as zero in required nodes*/
974int_t dinit3DLUstruct( int_t* myTreeIdxs, int_t* myZeroTrIdxs,
975 int_t* nodeCount, int_t** nodeList,
976 dLUstruct_t* LUstruct, gridinfo3d_t* grid3d);
977
978int_t dzSendLPanel(int_t k, int_t receiver,
979 dLUstruct_t* LUstruct, gridinfo3d_t* grid3d, SCT_t* SCT);
980int_t dzRecvLPanel(int_t k, int_t sender, double alpha,
981 double beta, double* Lval_buf,
982 dLUstruct_t* LUstruct, gridinfo3d_t* grid3d, SCT_t* SCT);
983int_t dzSendUPanel(int_t k, int_t receiver,
984 dLUstruct_t* LUstruct, gridinfo3d_t* grid3d, SCT_t* SCT);
985int_t dzRecvUPanel(int_t k, int_t sender, double alpha,
986 double beta, double* Uval_buf,
987 dLUstruct_t* LUstruct, gridinfo3d_t* grid3d, SCT_t* SCT);
988
989 /* from communication_aux.h */
990extern int_t dIBcast_LPanel (int_t k, int_t k0, int_t* lsub, double* lusup,
991 gridinfo_t *, int* msgcnt, MPI_Request *,
992 int **ToSendR, int_t *xsup, int );
993extern int_t dBcast_LPanel(int_t k, int_t k0, int_t* lsub, double* lusup,
994 gridinfo_t *, int* msgcnt, int **ToSendR,
995 int_t *xsup , SCT_t*, int);
996extern int_t dIBcast_UPanel(int_t k, int_t k0, int_t* usub, double* uval,
997 gridinfo_t *, int* msgcnt, MPI_Request *,
998 int *ToSendD, int );
999extern int_t dBcast_UPanel(int_t k, int_t k0, int_t* usub, double* uval,
1000 gridinfo_t *, int* msgcnt, int *ToSendD, SCT_t*, int);
1001extern int_t dIrecv_LPanel (int_t k, int_t k0, int_t* Lsub_buf,
1002 double* Lval_buf, gridinfo_t *,
1003 MPI_Request *, dLocalLU_t *, int);
1004extern int_t dIrecv_UPanel(int_t k, int_t k0, int_t* Usub_buf, double*,
1005 dLocalLU_t *, gridinfo_t*, MPI_Request *, int);
1006extern int_t dWait_URecv(MPI_Request *, int* msgcnt, SCT_t *);
1007extern int_t dWait_LRecv(MPI_Request*, int* msgcnt, int* msgcntsU,
1008 gridinfo_t *, SCT_t*);
1009extern int_t dISend_UDiagBlock(int_t k0, double *ublk_ptr, int_t size,
1010 MPI_Request *, gridinfo_t *, int);
1011extern int_t dRecv_UDiagBlock(int_t k0, double *ublk_ptr, int_t size,
1012 int_t src, gridinfo_t *, SCT_t*, int);
1013extern int_t dPackLBlock(int_t k, double* Dest, Glu_persist_t *,
1014 gridinfo_t *, dLocalLU_t *);
1015extern int_t dISend_LDiagBlock(int_t k0, double *lblk_ptr, int_t size,
1016 MPI_Request *, gridinfo_t *, int);
1017extern int_t dIRecv_UDiagBlock(int_t k0, double *ublk_ptr, int_t size,
1018 int_t src, MPI_Request *, gridinfo_t *,
1019 SCT_t*, int);
1020extern int_t dIRecv_LDiagBlock(int_t k0, double *L_blk_ptr, int_t size,
1021 int_t src, MPI_Request *, gridinfo_t*, SCT_t*, int);
1022extern int_t dUDiagBlockRecvWait( int_t k, int_t* IrecvPlcd_D, int_t* factored_L,
1023 MPI_Request *, gridinfo_t *, dLUstruct_t *, SCT_t *);
1024extern int_t LDiagBlockRecvWait( int_t k, int_t* factored_U, MPI_Request *, gridinfo_t *);
1025
1026#if (MPI_VERSION>2)
1027extern int_t dIBcast_UDiagBlock(int_t k, double *ublk_ptr, int_t size,
1028 MPI_Request *, gridinfo_t *);
1029extern int_t dIBcast_LDiagBlock(int_t k, double *lblk_ptr, int_t size,
1030 MPI_Request *, gridinfo_t *);
1031#endif
1032
1033 /* from trfCommWrapper.h */
1034extern int_t dDiagFactIBCast(int_t k, int_t k0,
1035 double *BlockUFactor, double *BlockLFactor,
1036 int_t* IrecvPlcd_D, MPI_Request *, MPI_Request *,
1037 MPI_Request *, MPI_Request *, gridinfo_t *,
1038 superlu_dist_options_t *, double thresh,
1039 dLUstruct_t *LUstruct, SuperLUStat_t *, int *info,
1040 SCT_t *, int tag_ub);
1041extern int_t dUPanelTrSolve( int_t k, double* BlockLFactor, double* bigV,
1042 int_t ldt, Ublock_info_t*, gridinfo_t *,
1043 dLUstruct_t *, SuperLUStat_t *, SCT_t *);
1044extern int_t dLPanelUpdate(int_t k, int_t* IrecvPlcd_D, int_t* factored_L,
1045 MPI_Request *, double* BlockUFactor, gridinfo_t *,
1046 dLUstruct_t *, SCT_t *);
1047extern int_t dUPanelUpdate(int_t k, int_t* factored_U, MPI_Request *,
1048 double* BlockLFactor, double* bigV,
1049 int_t ldt, Ublock_info_t*, gridinfo_t *,
1050 dLUstruct_t *, SuperLUStat_t *, SCT_t *);
1051extern int_t dIBcastRecvLPanel(int_t k, int_t k0, int* msgcnt,
1052 MPI_Request *, MPI_Request *,
1053 int_t* Lsub_buf, double* Lval_buf,
1054 int * factored, gridinfo_t *, dLUstruct_t *,
1055 SCT_t *, int tag_ub);
1056extern int_t dIBcastRecvUPanel(int_t k, int_t k0, int* msgcnt, MPI_Request *,
1057 MPI_Request *, int_t* Usub_buf, double* Uval_buf,
1058 gridinfo_t *, dLUstruct_t *, SCT_t *, int tag_ub);
1059extern int_t dWaitL(int_t k, int* msgcnt, int* msgcntU, MPI_Request *,
1060 MPI_Request *, gridinfo_t *, dLUstruct_t *, SCT_t *);
1061extern int_t dWaitU(int_t k, int* msgcnt, MPI_Request *, MPI_Request *,
1062 gridinfo_t *, dLUstruct_t *, SCT_t *);
1063extern int_t dLPanelTrSolve(int_t k, int_t* factored_L, double* BlockUFactor,
1064 gridinfo_t *, dLUstruct_t *);
1065
1066 /* from trfAux.h */
1067extern int getNsupers(int, Glu_persist_t *);
1068extern int_t initPackLUInfo(int_t nsupers, packLUInfo_t* packLUInfo);
1069extern int freePackLUInfo(packLUInfo_t* packLUInfo);
1070extern int_t dSchurComplementSetup(int_t k, int *msgcnt, Ublock_info_t*,
1071 Remain_info_t*, uPanelInfo_t *,
1072 lPanelInfo_t *, int_t*, int_t *, int_t *,
1073 double *bigU, int_t* Lsub_buf,
1074 double* Lval_buf, int_t* Usub_buf,
1075 double* Uval_buf, gridinfo_t *, dLUstruct_t *);
1076extern int_t dSchurComplementSetupGPU(int_t k, msgs_t* msgs, packLUInfo_t*,
1077 int_t*, int_t*, int_t*, gEtreeInfo_t*,
1078 factNodelists_t*, dscuBufs_t*,
1079 dLUValSubBuf_t* LUvsb, gridinfo_t *,
1080 dLUstruct_t *, HyP_t*);
1081extern double* dgetBigV(int_t, int_t);
1082extern double* dgetBigU(superlu_dist_options_t *,
1083 int_t, gridinfo_t *, dLUstruct_t *);
1084// permutation from superLU default
1085
1086 /* from treeFactorization.h */
1087extern int_t dLluBufInit(dLUValSubBuf_t*, dLUstruct_t *);
1088extern int_t dinitScuBufs(superlu_dist_options_t *,
1089 int_t ldt, int_t num_threads, int_t nsupers,
1090 dscuBufs_t*, dLUstruct_t*, gridinfo_t *);
1091extern int dfreeScuBufs(dscuBufs_t* scuBufs);
1092
1093#if 0 // NOT CALLED
1094// the generic tree factoring code
1095extern int_t treeFactor(
1096 int_t nnnodes, // number of nodes in the tree
1097 int_t *perm_c_supno, // list of nodes in the order of factorization
1098 commRequests_t *comReqs, // lists of communication requests
1099 dscuBufs_t *scuBufs, // contains buffers for schur complement update
1100 packLUInfo_t*packLUInfo,
1101 msgs_t*msgs,
1102 dLUValSubBuf_t* LUvsb,
1103 ddiagFactBufs_t *dFBuf,
1104 factStat_t *factStat,
1105 factNodelists_t *fNlists,
1106 superlu_dist_options_t *options,
1107 int_t * gIperm_c_supno,
1108 int_t ldt,
1109 dLUstruct_t *LUstruct, gridinfo3d_t * grid3d, SuperLUStat_t *stat,
1110 double thresh, SCT_t *SCT,
1111 int *info
1112);
1113#endif
1114
1115extern int_t dsparseTreeFactor(
1116 int_t nnodes, // number of nodes in the tree
1117 int_t *perm_c_supno, // list of nodes in the order of factorization
1118 treeTopoInfo_t* treeTopoInfo,
1119 commRequests_t *comReqs, // lists of communication requests
1120 dscuBufs_t *scuBufs, // contains buffers for schur complement update
1121 packLUInfo_t*packLUInfo,
1122 msgs_t*msgs,
1123 dLUValSubBuf_t* LUvsb,
1124 ddiagFactBufs_t *dFBuf,
1125 factStat_t *factStat,
1126 factNodelists_t *fNlists,
1127 superlu_dist_options_t *options,
1128 int_t * gIperm_c_supno,
1129 int_t ldt,
1130 dLUstruct_t *LUstruct, gridinfo3d_t * grid3d, SuperLUStat_t *stat,
1131 double thresh, SCT_t *SCT,
1132 int *info
1133);
1134
1135extern int_t ddenseTreeFactor(
1136 int_t nnnodes, // number of nodes in the tree
1137 int_t *perm_c_supno, // list of nodes in the order of factorization
1138 commRequests_t *comReqs, // lists of communication requests
1139 dscuBufs_t *scuBufs, // contains buffers for schur complement update
1140 packLUInfo_t*packLUInfo,
1141 msgs_t*msgs,
1142 dLUValSubBuf_t* LUvsb,
1143 ddiagFactBufs_t *dFBuf,
1144 factStat_t *factStat,
1145 factNodelists_t *fNlists,
1146 superlu_dist_options_t *options,
1147 int_t * gIperm_c_supno,
1148 int_t ldt,
1149 dLUstruct_t *LUstruct, gridinfo3d_t * grid3d, SuperLUStat_t *stat,
1150 double thresh, SCT_t *SCT, int tag_ub,
1151 int *info
1152);
1153
1154extern int_t dsparseTreeFactor_ASYNC(
1155 sForest_t* sforest,
1156 commRequests_t **comReqss, // lists of communication requests // size maxEtree level
1157 dscuBufs_t *scuBufs, // contains buffers for schur complement update
1158 packLUInfo_t*packLUInfo,
1159 msgs_t**msgss, // size=num Look ahead
1160 dLUValSubBuf_t** LUvsbs, // size=num Look ahead
1161 ddiagFactBufs_t **dFBufs, // size maxEtree level
1162 factStat_t *factStat,
1163 factNodelists_t *fNlists,
1164 gEtreeInfo_t* gEtreeInfo, // global etree info
1165 superlu_dist_options_t *options,
1166 int_t * gIperm_c_supno,
1167 int_t ldt,
1168 HyP_t* HyP,
1169 dLUstruct_t *LUstruct, gridinfo3d_t * grid3d, SuperLUStat_t *stat,
1170 double thresh, SCT_t *SCT, int tag_ub,
1171 int *info
1172);
1173extern dLUValSubBuf_t** dLluBufInitArr(int_t numLA, dLUstruct_t *LUstruct);
1174extern int dLluBufFreeArr(int_t numLA, dLUValSubBuf_t **LUvsbs);
1175extern ddiagFactBufs_t** dinitDiagFactBufsArr(int_t mxLeafNode, int_t ldt, gridinfo_t* grid);
1176extern int dfreeDiagFactBufsArr(int_t mxLeafNode, ddiagFactBufs_t** dFBufs);
1177extern int_t dinitDiagFactBufs(int_t ldt, ddiagFactBufs_t* dFBuf);
1178extern int_t checkRecvUDiag(int_t k, commRequests_t *comReqs,
1179 gridinfo_t *grid, SCT_t *SCT);
1180extern int_t checkRecvLDiag(int_t k, commRequests_t *comReqs, gridinfo_t *, SCT_t *);
1181
1182#if 0 // NOT CALLED
1183/* from ancFactorization.h (not called) */
1184extern int_t ancestorFactor(
1185 int_t ilvl, // level of factorization
1186 sForest_t* sforest,
1187 commRequests_t **comReqss, // lists of communication requests // size maxEtree level
1188 dscuBufs_t *scuBufs, // contains buffers for schur complement update
1189 packLUInfo_t*packLUInfo,
1190 msgs_t**msgss, // size=num Look ahead
1191 dLUValSubBuf_t** LUvsbs, // size=num Look ahead
1192 ddiagFactBufs_t **dFBufs, // size maxEtree level
1193 factStat_t *factStat,
1194 factNodelists_t *fNlists,
1195 gEtreeInfo_t* gEtreeInfo, // global etree info
1196 superlu_dist_options_t *options,
1197 int_t * gIperm_c_supno,
1198 int_t ldt,
1199 HyP_t* HyP,
1200 dLUstruct_t *LUstruct, gridinfo3d_t * grid3d, SuperLUStat_t *stat,
1201 double thresh, SCT_t *SCT, int tag_ub, int *info
1202);
1203#endif
1204
1205/*== end 3D prototypes ===================*/
1206
1207
1208#ifdef __cplusplus
1209 }
1210#endif
1211
1212#endif /* __SUPERLU_dDEFS */
1213
j
int j
Definition: dutil_dist.c:248
dtrtri_
#define dtrtri_
Definition: superlu_FCnames.h:161
superlu_dger
int superlu_dger(const int m, const int n, const double alpha, const double *x, const int incx, const double *y, const int incy, double *a, const int lda)
dLluBufInit
int_t dLluBufInit(dLUValSubBuf_t *, dLUstruct_t *)
dscatter3dLPanels
int_t dscatter3dLPanels(int_t nsupers, dLUstruct_t *LUstruct, gridinfo3d_t *grid3d)
dstatic_schedule
int dstatic_schedule(superlu_dist_options_t *, int, int, dLUstruct_t *, gridinfo_t *, SuperLUStat_t *, int_t *, int_t *, int *)
Definition: dstatic_schedule.c:46
dSolveFinalize
void dSolveFinalize(superlu_dist_options_t *, dSOLVEstruct_t *)
Release the resources used for the solution phase.
Definition: pdutil.c:874
dScaleAddId_CompRowLoc_Matrix_dist
void dScaleAddId_CompRowLoc_Matrix_dist(SuperMatrix *, double)
Scale and add I: scales a matrix and adds an identity. A_{i,j} = c * A_{i,j} + \delta_{i,...
Definition: dutil_dist.c:356
initPackLUInfo
int_t initPackLUInfo(int_t nsupers, packLUInfo_t *packLUInfo)
Definition: treeFactorization.c:367
dScaleAdd_CompRowLoc_Matrix_dist
void dScaleAdd_CompRowLoc_Matrix_dist(SuperMatrix *, SuperMatrix *, double)
Scale and add: adds a scalar multiple of one matrix to another. A_{i,j} = c * A_{i,...
Definition: dutil_dist.c:381
dscatter_l
void dscatter_l(int ib, int ljb, int nsupc, int_t iukp, int_t *xsup, int klst, int nbrow, int_t lptr, int temp_nbrow, int_t *usub, int_t *lsub, double *tempv, int *indirect_thread, int *indirect2, int_t **Lrowind_bc_ptr, double **Lnzval_bc_ptr, gridinfo_t *grid)
dDeAllocLlu_3d
int dDeAllocLlu_3d(int_t n, dLUstruct_t *, gridinfo3d_t *)
Definition: dutil_dist.c:450
pdReDistribute_B_to_X
int_t pdReDistribute_B_to_X(double *B, int_t m_loc, int nrhs, int_t ldb, int_t fst_row, int_t *ilsum, double *x, dScalePermstruct_t *, Glu_persist_t *, gridinfo_t *, dSOLVEstruct_t *)
Definition: pdgstrs.c:347
dPackLBlock
int_t dPackLBlock(int_t k, double *Dest, Glu_persist_t *, gridinfo_t *, dLocalLU_t *)
updateDirtyBit
int updateDirtyBit(int_t k0, HyP_t *HyP, gridinfo_t *grid)
Definition: sec_structs.c:618
pdGetDiagU
void pdGetDiagU(int_t, dLUstruct_t *, gridinfo_t *, double *)
Definition: pdGetDiagU.c:66
superlu_daxpy
int superlu_daxpy(const int n, const double alpha, const double *x, const int incx, double *y, const int incy)
dcreate_matrix
int dcreate_matrix(SuperMatrix *, int, double **, int *, double **, int *, FILE *, gridinfo_t *)
Definition: dcreate_matrix.c:349
pdgsrfs
void pdgsrfs(int_t, SuperMatrix *, double, dLUstruct_t *, dScalePermstruct_t *, gridinfo_t *, double[], int_t, double[], int_t, int, dSOLVEstruct_t *, double *, SuperLUStat_t *, int *)
dQuerySpace_dist
int_t dQuerySpace_dist(int_t, dLUstruct_t *, gridinfo_t *, SuperLUStat_t *, superlu_dist_mem_usage_t *)
Definition: dmemory_dist.c:73
dcreate_matrix3d
int dcreate_matrix3d(SuperMatrix *A, int nrhs, double **rhs, int *ldb, double **x, int *ldx, FILE *fp, gridinfo3d_t *grid3d)
dGenCSCLblocks
void dGenCSCLblocks(int, int_t, gridinfo_t *, Glu_persist_t *, dLocalLU_t *, double **, int_t **, int_t **, int_t *, int_t *)
Definition: dutil_dist.c:901
dZeroUblocks
void dZeroUblocks(int iam, int n, gridinfo_t *, dLUstruct_t *)
Sets all entries of matrix U to zero.
Definition: dutil_dist.c:1249
pdgsmv_init
void pdgsmv_init(SuperMatrix *, int_t *, gridinfo_t *, pdgsmv_comm_t *)
Definition: pdgsmv.c:27
dISend_UDiagBlock
int_t dISend_UDiagBlock(int_t k0, double *ublk_ptr, int_t size, MPI_Request *, gridinfo_t *, int)
dfreeScuBufs
int dfreeScuBufs(dscuBufs_t *scuBufs)
dUDiagBlockRecvWait
int_t dUDiagBlockRecvWait(int_t k, int_t *IrecvPlcd_D, int_t *factored_L, MPI_Request *, gridinfo_t *, dLUstruct_t *, SCT_t *)
pdgsmv
void pdgsmv(int_t, SuperMatrix *, gridinfo_t *, pdgsmv_comm_t *, double x[], double ax[])
Definition: pdgsmv.c:235
dreduceAncestors3d
int_t dreduceAncestors3d(int_t sender, int_t receiver, int_t nnodes, int_t *nodeList, double *Lval_buf, double *Uval_buf, dLUstruct_t *LUstruct, gridinfo3d_t *grid3d, SCT_t *SCT)
dGenCSRLblocks
void dGenCSRLblocks(int, int_t, gridinfo_t *, Glu_persist_t *, dLocalLU_t *, double **, int_t **, int_t **, int_t *, int_t *)
Definition: dutil_dist.c:1054
dTrs2_GatherTrsmScatter
int_t dTrs2_GatherTrsmScatter(int_t klst, int_t iukp, int_t rukp, int_t *usub, double *uval, double *tempv, int_t knsupc, int nsupr, double *lusup, Glu_persist_t *Glu_persist)
Definition: pdgstrf2.c:716
dcollect3dUpanels
int_t dcollect3dUpanels(int_t layer, int_t nsupers, dLUstruct_t *LUstruct, gridinfo3d_t *grid3d)
dIRecv_UDiagBlock
int_t dIRecv_UDiagBlock(int_t k0, double *ublk_ptr, int_t size, int_t src, MPI_Request *, gridinfo_t *, SCT_t *, int)
dinitScuBufs
int_t dinitScuBufs(int_t ldt, int_t num_threads, int_t nsupers, dscuBufs_t *, dLUstruct_t *, gridinfo_t *)
ddenseTreeFactor
int_t ddenseTreeFactor(int_t nnnodes, int_t *perm_c_supno, commRequests_t *comReqs, dscuBufs_t *scuBufs, packLUInfo_t *packLUInfo, msgs_t *msgs, dLUValSubBuf_t *LUvsb, ddiagFactBufs_t *dFBuf, factStat_t *factStat, factNodelists_t *fNlists, superlu_dist_options_t *options, int_t *gIperm_c_supno, int_t ldt, dLUstruct_t *LUstruct, gridinfo3d_t *grid3d, SuperLUStat_t *stat, double thresh, SCT_t *SCT, int tag_ub, int *info)
dinit3DLUstructForest
void dinit3DLUstructForest(int_t *myTreeIdxs, int_t *myZeroTrIdxs, sForest_t **sForests, dLUstruct_t *LUstruct, gridinfo3d_t *grid3d)
scuStatUpdate
int_t scuStatUpdate(int_t knsupc, HyP_t *HyP, SCT_t *SCT, SuperLUStat_t *stat)
Definition: sec_structs.c:635
pdgsmv_AXglobal_abs
int pdgsmv_AXglobal_abs(int_t, int_t[], double[], int_t[], double[], double[])
Definition: pdgsmv_AXglobal.c:285
dgatherAllFactoredLU
int_t dgatherAllFactoredLU(trf3Dpartition_t *trf3Dpartition, dLUstruct_t *LUstruct, gridinfo3d_t *grid3d, SCT_t *SCT)
dreadrb_dist
void dreadrb_dist(int, FILE *, int_t *, int_t *, int_t *, double **, int_t **, int_t **)
Definition: dreadrb.c:275
pdCompute_Diag_Inv
void pdCompute_Diag_Inv(int_t, dLUstruct_t *, gridinfo_t *, SuperLUStat_t *, int *)
Definition: pdgstrs.c:839
pdgsequ
void pdgsequ(SuperMatrix *, double *, double *, double *, double *, double *, int_t *, gridinfo_t *)
Definition: pdgsequ.c:86
dscal_
int dscal_(const int *n, const double *alpha, double *dx, const int *incx)
dWaitL
int_t dWaitL(int_t k, int *msgcnt, int *msgcntU, MPI_Request *, MPI_Request *, gridinfo_t *, dLUstruct_t *, SCT_t *)
dtrsm_
int dtrsm_(const char *, const char *, const char *, const char *, const int *, const int *, const double *, const double *, const int *, double *, const int *)
dGenXtrue_dist
void dGenXtrue_dist(int_t, int_t, double *, int_t)
Definition: dutil_dist.c:486
d3D_printMemUse
void d3D_printMemUse(trf3Dpartition_t *trf3Dpartition, dLUstruct_t *LUstruct, gridinfo3d_t *grid3d)
Definition: dmemory_dist.c:243
dgather_l
void dgather_l(int_t num_LBlk, int_t knsupc, Remain_info_t *L_info, double *lval, int_t LD_lval, double *L_buff)
sp_dgemv_dist
int sp_dgemv_dist(char *, double, SuperMatrix *, double *, int, double, double *, int)
SpGEMV.
Definition: dsp_blas2_dist.c:391
dzRecvLPanel
int_t dzRecvLPanel(int_t k, int_t sender, double alpha, double beta, double *Lval_buf, dLUstruct_t *LUstruct, gridinfo3d_t *grid3d, SCT_t *SCT)
dIBcast_UPanel
int_t dIBcast_UPanel(int_t k, int_t k0, int_t *usub, double *uval, gridinfo_t *, int *msgcnt, MPI_Request *, int *ToSendD, int)
getNsupers
int getNsupers(int, Glu_persist_t *)
Definition: trfAux.c:42
dfreeDiagFactBufsArr
int dfreeDiagFactBufsArr(int_t mxLeafNode, ddiagFactBufs_t **dFBufs)
pdlangs
double pdlangs(char *, SuperMatrix *, gridinfo_t *)
Definition: pdlangs.c:65
daxpy_
int daxpy_(const int *n, const double *alpha, const double *x, const int *incx, double *y, const int *incy)
dinitTrf3Dpartition
trf3Dpartition_t * dinitTrf3Dpartition(int_t nsupers, superlu_dist_options_t *options, dLUstruct_t *LUstruct, gridinfo3d_t *grid3d)
dGenCOOLblocks
void dGenCOOLblocks(int, int_t, gridinfo_t *, Glu_persist_t *, dLocalLU_t *, int_t **, int_t **, double **, int_t *, int_t *)
Definition: dutil_dist.c:794
dSchurComplementSetup
int_t dSchurComplementSetup(int_t k, int *msgcnt, Ublock_info_t *, Remain_info_t *, uPanelInfo_t *, lPanelInfo_t *, int_t *, int_t *, int_t *, double *bigU, int_t *Lsub_buf, double *Lval_buf, int_t *Usub_buf, double *Uval_buf, gridinfo_t *, dLUstruct_t *)
dDestroy_A3d_gathered_on_2d
void dDestroy_A3d_gathered_on_2d(dSOLVEstruct_t *, gridinfo3d_t *)
Definition: pdutil.c:894
dzSendLPanel
int_t dzSendLPanel(int_t k, int_t receiver, dLUstruct_t *LUstruct, gridinfo3d_t *grid3d, SCT_t *SCT)
dp3dScatter
int_t dp3dScatter(int_t n, dLUstruct_t *LUstruct, gridinfo3d_t *grid3d)
dCreate_CompRowLoc_Matrix_dist
void dCreate_CompRowLoc_Matrix_dist(SuperMatrix *, int_t, int_t, int_t, int_t, int_t, double *, int_t *, int_t *, Stype_t, Dtype_t, Mtype_t)
Free_HyP
void Free_HyP(HyP_t *HyP)
Definition: sec_structs.c:594
dlsum_bmod_inv
void dlsum_bmod_inv(double *, double *, double *, double *, int, int_t, int *bmod, int_t *, Ucb_indptr_t **, int_t **, int_t *, gridinfo_t *, dLocalLU_t *, SuperLUStat_t **, int_t *, int_t *, int_t, int_t, int, int)
Definition: pdgstrs_lsum.c:1355
dblock_gemm_scatterTopLeft
int_t dblock_gemm_scatterTopLeft(int_t lb, int_t j, double *bigV, int_t knsupc, int_t klst, int_t *lsub, int_t *usub, int_t ldt, int *indirect, int *indirect2, HyP_t *HyP, dLUstruct_t *, gridinfo_t *, SCT_t *SCT, SuperLUStat_t *)
dlsum_fmod
void dlsum_fmod(double *, double *, double *, double *, int, int, int_t, int *fmod, int_t, int_t, int_t, int_t *, gridinfo_t *, dLocalLU_t *, MPI_Request[], SuperLUStat_t *)
Definition: pdgstrs_lsum.c:62
dIRecv_LDiagBlock
int_t dIRecv_LDiagBlock(int_t k0, double *L_blk_ptr, int_t size, int_t src, MPI_Request *, gridinfo_t *, SCT_t *, int)
dScalePermstructFree
void dScalePermstructFree(dScalePermstruct_t *)
Deallocate ScalePermstruct.
Definition: dutil_dist.c:409
dzSendUPanel
int_t dzSendUPanel(int_t k, int_t receiver, dLUstruct_t *LUstruct, gridinfo3d_t *grid3d, SCT_t *SCT)
pdPermute_Dense_Matrix
int pdPermute_Dense_Matrix(int_t, int_t, int_t[], int_t[], double[], int, double[], int, int, gridinfo_t *)
Permute the distributed dense matrix: B <= perm(X). perm[i] = j means the i-th row of X is in the j-t...
Definition: pdutil.c:294
dRgather_U
void dRgather_U(int_t k, int_t jj0, int_t *usub, double *uval, double *bigU, gEtreeInfo_t *, Glu_persist_t *, gridinfo_t *, HyP_t *, int_t *myIperm, int_t *iperm_c_supno, int_t *perm_u)
dlsum_bmod_inv_master
void dlsum_bmod_inv_master(double *, double *, double *, double *, int, int_t, int *bmod, int_t *, Ucb_indptr_t **, int_t **, int_t *, gridinfo_t *, dLocalLU_t *, SuperLUStat_t **, int_t, int_t, int, int)
Definition: pdgstrs_lsum.c:1819
checkRecvUDiag
int_t checkRecvUDiag(int_t k, commRequests_t *comReqs, gridinfo_t *grid, SCT_t *SCT)
Definition: treeFactorization.c:401
pdgsmv_AXglobal
int pdgsmv_AXglobal(int_t, int_t[], double[], int_t[], double[], double[])
Definition: pdgsmv_AXglobal.c:259
dAllocLlu_3d
int_t dAllocLlu_3d(int_t nsupers, dLUstruct_t *LUstruct, gridinfo3d_t *grid3d)
pdgstrs2
void pdgstrs2(int_t m, int_t k0, int_t k, Glu_persist_t *Glu_persist, gridinfo_t *grid, dLocalLU_t *Llu, SuperLUStat_t *stat)
dgemv_
void dgemv_(const char *, const int *, const int *, const double *, const double *a, const int *, const double *, const int *, const double *, double *, const int *)
dinitDiagFactBufsArr
ddiagFactBufs_t ** dinitDiagFactBufsArr(int_t mxLeafNode, int_t ldt, gridinfo_t *grid)
dIBcastRecvLPanel
int_t dIBcastRecvLPanel(int_t k, int_t k0, int *msgcnt, MPI_Request *, MPI_Request *, int_t *Lsub_buf, double *Lval_buf, int_t *factored, gridinfo_t *, dLUstruct_t *, SCT_t *, int tag_ub)
dinf_norm_error_dist
void dinf_norm_error_dist(int_t, int_t, double *, int_t, double *, int_t, gridinfo_t *)
Check the inf-norm of the error vector.
Definition: dutil_dist.c:522
ddistribute
float ddistribute(fact_t, int_t, SuperMatrix *, Glu_freeable_t *, dLUstruct_t *, gridinfo_t *)
Definition: ddistribute.c:65
pdgsmv_AXglobal_setup
int pdgsmv_AXglobal_setup(SuperMatrix *, Glu_persist_t *, gridinfo_t *, int_t *, int_t *[], double *[], int_t *[], int_t[])
pdgsrfs_ABXglobal
void pdgsrfs_ABXglobal(int_t, SuperMatrix *, double, dLUstruct_t *, gridinfo_t *, double *, int_t, double *, int_t, int, double *, SuperLUStat_t *, int *)
Definition: pdgsrfs_ABXglobal.c:125
dGatherNRformat_loc3d
void dGatherNRformat_loc3d(fact_t Fact, NRformat_loc *A, double *B, int ldb, int nrhs, gridinfo3d_t *grid3d, NRformat_loc3d **)
dcollect3dLpanels
int_t dcollect3dLpanels(int_t layer, int_t nsupers, dLUstruct_t *LUstruct, gridinfo3d_t *grid3d)
dcreate_matrix_postfix3d
int dcreate_matrix_postfix3d(SuperMatrix *A, int nrhs, double **rhs, int *ldb, double **x, int *ldx, FILE *fp, char *postfix, gridinfo3d_t *grid3d)
Definition: dcreate_matrix3d.c:72
dreadMM_dist
void dreadMM_dist(FILE *, int_t *, int_t *, int_t *, double **, int_t **, int_t **)
Definition: dreadMM.c:38
dInit_HyP
void dInit_HyP(HyP_t *HyP, dLocalLU_t *Llu, int_t mcb, int_t mrb)
dlsum_fmod_inv_master
void dlsum_fmod_inv_master(double *, double *, double *, double *, int, int, int_t, int *fmod, int_t, int_t *, gridinfo_t *, dLocalLU_t *, SuperLUStat_t **, int_t, int_t, int_t, int_t, int, int)
Definition: pdgstrs_lsum.c:956
freePackLUInfo
int freePackLUInfo(packLUInfo_t *packLUInfo)
Definition: treeFactorization.c:376
dDestroy_Tree
void dDestroy_Tree(int_t, gridinfo_t *, dLUstruct_t *)
Destroy broadcast and reduction trees used in triangular solve.
Definition: pdutil.c:437
dblock_gemm_scatter
void dblock_gemm_scatter(int_t lb, int_t j, Ublock_info_t *Ublock_info, Remain_info_t *Remain_info, double *L_mat, int ldl, double *U_mat, int ldu, double *bigV, int_t knsupc, int_t klst, int_t *lsub, int_t *usub, int_t ldt, int_t thread_id, int *indirect, int *indirect2, int_t **Lrowind_bc_ptr, double **Lnzval_bc_ptr, int_t **Ufstnz_br_ptr, double **Unzval_br_ptr, int_t *xsup, gridinfo_t *, SuperLUStat_t *)
dscatter3dUPanels
int_t dscatter3dUPanels(int_t nsupers, dLUstruct_t *LUstruct, gridinfo3d_t *grid3d)
dDestroy_LU
void dDestroy_LU(int_t, gridinfo_t *, dLUstruct_t *)
Destroy distributed L & U matrices.
Definition: pdutil.c:485
pdgssvx3d
void pdgssvx3d(superlu_dist_options_t *, SuperMatrix *, dScalePermstruct_t *, double B[], int ldb, int nrhs, gridinfo3d_t *, dLUstruct_t *, dSOLVEstruct_t *, double *berr, SuperLUStat_t *, int *info)
Definition: pdgssvx3d.c:501
dWait_LRecv
int_t dWait_LRecv(MPI_Request *, int *msgcnt, int *msgcntsU, gridinfo_t *, SCT_t *)
dCopy_CompRowLoc_Matrix_dist
void dCopy_CompRowLoc_Matrix_dist(SuperMatrix *, SuperMatrix *)
Definition: dutil_dist.c:324
dCreate_CompCol_Matrix_dist
void dCreate_CompCol_Matrix_dist(SuperMatrix *, int_t, int_t, int_t, double *, int_t *, int_t *, Stype_t, Dtype_t, Mtype_t)
dgatherAllFactoredLUFr
int_t dgatherAllFactoredLUFr(int_t *myZeroTrIdxs, sForest_t *sForests, dLUstruct_t *LUstruct, gridinfo3d_t *grid3d, SCT_t *SCT)
dTrs2_ScatterU
int_t dTrs2_ScatterU(int_t iukp, int_t rukp, int_t klst, int_t nsupc, int_t ldu, int_t *usub, double *uval, double *tempv)
Definition: pdgstrf2.c:694
dDiagFactIBCast
int_t dDiagFactIBCast(int_t k, int_t k0, double *BlockUFactor, double *BlockLFactor, int_t *IrecvPlcd_D, MPI_Request *, MPI_Request *, MPI_Request *, MPI_Request *, gridinfo_t *, superlu_dist_options_t *, double thresh, dLUstruct_t *LUstruct, SuperLUStat_t *, int *info, SCT_t *, int tag_ub)
dsparseTreeFactor
int_t dsparseTreeFactor(int_t nnodes, int_t *perm_c_supno, treeTopoInfo_t *treeTopoInfo, commRequests_t *comReqs, dscuBufs_t *scuBufs, packLUInfo_t *packLUInfo, msgs_t *msgs, dLUValSubBuf_t *LUvsb, ddiagFactBufs_t *dFBuf, factStat_t *factStat, factNodelists_t *fNlists, superlu_dist_options_t *options, int_t *gIperm_c_supno, int_t ldt, dLUstruct_t *LUstruct, gridinfo3d_t *grid3d, SuperLUStat_t *stat, double thresh, SCT_t *SCT, int *info)
dPrint_Dense_Matrix_dist
void dPrint_Dense_Matrix_dist(SuperMatrix *)
dblock_gemm_scatterTopRight
int_t dblock_gemm_scatterTopRight(int_t lb, int_t j, double *bigV, int_t knsupc, int_t klst, int_t *lsub, int_t *usub, int_t ldt, int *indirect, int *indirect2, HyP_t *HyP, dLUstruct_t *, gridinfo_t *, SCT_t *SCT, SuperLUStat_t *)
dRgather_L
void dRgather_L(int_t k, int_t *lsub, double *lusup, gEtreeInfo_t *, Glu_persist_t *, gridinfo_t *, HyP_t *, int_t *myIperm, int_t *iperm_c_supno)
dblock_gemm_scatterBottomLeft
int_t dblock_gemm_scatterBottomLeft(int_t lb, int_t j, double *bigV, int_t knsupc, int_t klst, int_t *lsub, int_t *usub, int_t ldt, int *indirect, int *indirect2, HyP_t *HyP, dLUstruct_t *, gridinfo_t *, SCT_t *SCT, SuperLUStat_t *)
pdgstrs
void pdgstrs(int_t, dLUstruct_t *, dScalePermstruct_t *, gridinfo_t *, double *, int_t, int_t, int_t, int, dSOLVEstruct_t *, SuperLUStat_t *, int *)
Definition: pdgstrs.c:1030
dScatter_B3d
int dScatter_B3d(NRformat_loc3d *A3d, gridinfo3d_t *grid3d)
dgather_u
void dgather_u(int_t num_u_blks, Ublock_info_t *Ublock_info, int_t *usub, double *uval, double *bigU, int_t ldu, int_t *xsup, int_t klst)
file_dPrint_CompRowLoc_Matrix_dist
int file_dPrint_CompRowLoc_Matrix_dist(FILE *fp, SuperMatrix *A)
dIBcastRecvUPanel
int_t dIBcastRecvUPanel(int_t k, int_t k0, int *msgcnt, MPI_Request *, MPI_Request *, int_t *Usub_buf, double *Uval_buf, gridinfo_t *, dLUstruct_t *, SCT_t *, int tag_ub)
dZeroLblocks
void dZeroLblocks(int, int, gridinfo_t *, dLUstruct_t *)
Sets all entries of matrix L to zero.
Definition: dutil_dist.c:619
sp_dtrsv_dist
int sp_dtrsv_dist(char *, char *, char *, SuperMatrix *, SuperMatrix *, double *, int *)
Definition: dsp_blas2_dist.c:95
dldperm_dist
int dldperm_dist(int, int, int_t, int_t[], int_t[], double[], int_t *, double[], double[])
Definition: dldperm_dist.c:96
dLpanelUpdate
int_t dLpanelUpdate(int_t off0, int_t nsupc, double *ublk_ptr, int_t ld_ujrow, double *lusup, int_t nsupr, SCT_t *)
duser_malloc_dist
void * duser_malloc_dist(int_t, int_t)
Definition: dmemory_dist.c:30
Local_Dgstrf2
void Local_Dgstrf2(superlu_dist_options_t *options, int_t k, double thresh, double *BlockUFactor, Glu_persist_t *, gridinfo_t *, dLocalLU_t *, SuperLUStat_t *, int *info, SCT_t *)
Definition: pdgstrf2.c:402
LDiagBlockRecvWait
int_t LDiagBlockRecvWait(int_t k, int_t *factored_U, MPI_Request *, gridinfo_t *)
Definition: communication_aux.c:218
pdgstrs_Bglobal
void pdgstrs_Bglobal(int_t, dLUstruct_t *, gridinfo_t *, double *, int_t, int, SuperLUStat_t *, int *)
Definition: pdgstrs_Bglobal.c:104
dWait_URecv
int_t dWait_URecv(MPI_Request *, int *msgcnt, SCT_t *)
dLUstructFree
void dLUstructFree(dLUstruct_t *)
Deallocate LUstruct.
Definition: pdutil.c:419
dPrint_CompRowLoc_Matrix_dist
int dPrint_CompRowLoc_Matrix_dist(SuperMatrix *)
dCreate_SuperNode_Matrix_dist
void dCreate_SuperNode_Matrix_dist(SuperMatrix *, int_t, int_t, int_t, double *, int_t *, int_t *, int_t *, int_t *, int_t *, Stype_t, Dtype_t, Mtype_t)
pdgssvx_ABglobal
void pdgssvx_ABglobal(superlu_dist_options_t *, SuperMatrix *, dScalePermstruct_t *, double *, int, int, gridinfo_t *, dLUstruct_t *, double *, SuperLUStat_t *, int *)
dPrint_CompCol_Matrix_dist
void dPrint_CompCol_Matrix_dist(SuperMatrix *)
dtrsv_
int dtrsv_(char *, char *, char *, int *, double *, int *, double *, int *)
dIrecv_LPanel
int_t dIrecv_LPanel(int_t k, int_t k0, int_t *Lsub_buf, double *Lval_buf, gridinfo_t *, MPI_Request *, dLocalLU_t *, int)
pdgstrf2_trsm
void pdgstrf2_trsm(superlu_dist_options_t *options, int_t k0, int_t k, double thresh, Glu_persist_t *, gridinfo_t *, dLocalLU_t *, MPI_Request *, int tag_ub, SuperLUStat_t *, int *info)
Definition: pdgstrf2.c:143
pdgstrs_init
int_t pdgstrs_init(int_t, int_t, int_t, int_t, int_t[], int_t[], gridinfo_t *grid, Glu_persist_t *, dSOLVEstruct_t *)
Definition: pdutil.c:657
treeFactor
int_t treeFactor(int_t nnnodes, int_t *perm_c_supno, commRequests_t *comReqs, dscuBufs_t *scuBufs, packLUInfo_t *packLUInfo, msgs_t *msgs, dLUValSubBuf_t *LUvsb, ddiagFactBufs_t *dFBuf, factStat_t *factStat, factNodelists_t *fNlists, superlu_dist_options_t *options, int_t *gIperm_c_supno, int_t ldt, dLUstruct_t *LUstruct, gridinfo3d_t *grid3d, SuperLUStat_t *stat, double thresh, SCT_t *SCT, int *info)
dISend_LDiagBlock
int_t dISend_LDiagBlock(int_t k0, double *lblk_ptr, int_t size, MPI_Request *, gridinfo_t *, int)
dgatherFactoredLU
int_t dgatherFactoredLU(int_t sender, int_t receiver, int_t nnodes, int_t *nodeList, dLUValSubBuf_t *LUvsb, dLUstruct_t *LUstruct, gridinfo3d_t *grid3d, SCT_t *SCT)
dinitDiagFactBufs
int_t dinitDiagFactBufs(int_t ldt, ddiagFactBufs_t *dFBuf)
pdgstrf
int_t pdgstrf(superlu_dist_options_t *, int, int, double anorm, dLUstruct_t *, gridinfo_t *, SuperLUStat_t *, int *)
Definition: pdgstrf.c:242
dgetBigU
double * dgetBigU(int_t, gridinfo_t *, dLUstruct_t *)
dSchurComplementSetupGPU
int_t dSchurComplementSetupGPU(int_t k, msgs_t *msgs, packLUInfo_t *, int_t *, int_t *, int_t *, gEtreeInfo_t *, factNodelists_t *, dscuBufs_t *, dLUValSubBuf_t *LUvsb, gridinfo_t *, dLUstruct_t *, HyP_t *)
dAllocGlu_3d
int dAllocGlu_3d(int_t n, int_t nsupers, dLUstruct_t *)
Definition: dutil_dist.c:432
dBcast_UPanel
int_t dBcast_UPanel(int_t k, int_t k0, int_t *usub, double *uval, gridinfo_t *, int *msgcnt, int *ToSendD, SCT_t *, int)
file_Printdouble5
int file_Printdouble5(FILE *, char *, int_t, double *)
Definition: dutil_dist.c:555
dreadhb_dist
void dreadhb_dist(int, FILE *, int_t *, int_t *, int_t *, double **, int_t **, int_t **)
Definition: dreadhb.c:107
dTrs2_GatherU
int_t dTrs2_GatherU(int_t iukp, int_t rukp, int_t klst, int_t nsupc, int_t ldu, int_t *usub, double *uval, double *tempv)
Definition: pdgstrf2.c:669
pdgstrs2_omp
void pdgstrs2_omp(int_t k0, int_t k, Glu_persist_t *, gridinfo_t *, dLocalLU_t *, Ublock_info_t *, SuperLUStat_t *)
Definition: pdgstrf2.c:762
superlu_dgemv
int superlu_dgemv(const char *trans, const int m, const int n, const double alpha, const double *a, const int lda, const double *x, const int incx, const double beta, double *y, const int incy)
dCompRow_to_CompCol_dist
void dCompRow_to_CompCol_dist(int_t, int_t, int_t, double *, int_t *, int_t *, double **, int_t **, int_t **)
dscatter_u
void dscatter_u(int ib, int jb, int nsupc, int_t iukp, int_t *xsup, int klst, int nbrow, int_t lptr, int temp_nbrow, int_t *lsub, int_t *usub, double *tempv, int_t **Ufstnz_br_ptr, double **Unzval_br_ptr, gridinfo_t *grid)
dgsequ_dist
void dgsequ_dist(SuperMatrix *, double *, double *, double *, double *, double *, int_t *)
Definition: dgsequ_dist.c:84
dRecv_UDiagBlock
int_t dRecv_UDiagBlock(int_t k0, double *ublk_ptr, int_t size, int_t src, gridinfo_t *, SCT_t *, int)
dLUstructInit
void dLUstructInit(const int_t, dLUstruct_t *)
Allocate storage in LUstruct.
Definition: pdutil.c:405
dCreate_Dense_Matrix_dist
void dCreate_Dense_Matrix_dist(SuperMatrix *, int_t, int_t, double *, int_t, Stype_t, Dtype_t, Mtype_t)
dallocateA_dist
void dallocateA_dist(int_t, int_t, double **, int_t **, int_t **)
Definition: dmemory_dist.c:147
dlsum_fmod_inv
void dlsum_fmod_inv(double *, double *, double *, double *, int, int_t, int *fmod, int_t *, gridinfo_t *, dLocalLU_t *, SuperLUStat_t **, int_t *, int_t *, int_t, int_t, int_t, int_t, int, int)
Definition: pdgstrs_lsum.c:416
dcreate_matrix_rb
int dcreate_matrix_rb(SuperMatrix *, int, double **, int *, double **, int *, FILE *, gridinfo_t *)
pddistribute
float pddistribute(fact_t, int_t, SuperMatrix *, dScalePermstruct_t *, Glu_freeable_t *, dLUstruct_t *, gridinfo_t *)
Definition: pddistribute.c:328
superlu_dtrsm
int superlu_dtrsm(const char *sideRL, const char *uplo, const char *transa, const char *diag, const int m, const int n, const double alpha, const double *a, const int lda, double *b, const int ldb)
pxgstrs_finalize
void pxgstrs_finalize(pxgstrs_comm_t *)
Definition: util.c:266
dlsum_bmod
void dlsum_bmod(double *, double *, double *, int, int_t, int *bmod, int_t *, Ucb_indptr_t **, int_t **, int_t *, gridinfo_t *, dLocalLU_t *, MPI_Request[], SuperLUStat_t *)
Definition: pdgstrs_lsum.c:246
pdinf_norm_error
void pdinf_norm_error(int, int_t, int_t, double[], int_t, double[], int_t, MPI_Comm)
Check the inf-norm of the error vector.
Definition: pdutil.c:921
dread_binary
int dread_binary(FILE *, int_t *, int_t *, int_t *, double **, int_t **, int_t **)
Definition: dbinary_io.c:4
doubleMalloc_dist
double * doubleMalloc_dist(int_t)
Definition: dmemory_dist.c:155
pdCompRow_loc_to_CompCol_global
int pdCompRow_loc_to_CompCol_global(int_t, SuperMatrix *, gridinfo_t *, SuperMatrix *)
Gather A from the distributed compressed row format to global A in compressed column format.
Definition: pdutil.c:32
dSolveInit
int dSolveInit(superlu_dist_options_t *, SuperMatrix *, int_t[], int_t[], int_t, dLUstruct_t *, gridinfo_t *, dSOLVEstruct_t *)
Initialize the data structure for the solution phase.
Definition: pdutil.c:786
dIBcast_LPanel
int_t dIBcast_LPanel(int_t k, int_t k0, int_t *lsub, double *lusup, gridinfo_t *, int *msgcnt, MPI_Request *, int **ToSendR, int_t *xsup, int)
pdgstrf2
void pdgstrf2(superlu_dist_options_t *, int_t nsupers, int_t k0, int_t k, double thresh, Glu_persist_t *, gridinfo_t *, dLocalLU_t *, MPI_Request *, int, SuperLUStat_t *, int *)
dPrintUblocks
void dPrintUblocks(int, int_t, gridinfo_t *, Glu_persist_t *, dLocalLU_t *)
Print the blocks in the factored matrix U.
Definition: dutil_dist.c:1208
dreduceAllAncestors3d
int dreduceAllAncestors3d(int_t ilvl, int_t *myNodeCount, int_t **treePerm, dLUValSubBuf_t *LUvsb, dLUstruct_t *LUstruct, gridinfo3d_t *grid3d, SCT_t *SCT)
dzeroSetLU
int_t dzeroSetLU(int_t nnodes, int_t *nodeList, dLUstruct_t *, gridinfo3d_t *)
dWaitU
int_t dWaitU(int_t k, int *msgcnt, MPI_Request *, MPI_Request *, gridinfo_t *, dLUstruct_t *, SCT_t *)
dcreate_matrix_postfix
int dcreate_matrix_postfix(SuperMatrix *, int, double **, int *, double **, int *, FILE *, char *, gridinfo_t *)
Definition: dcreate_matrix.c:76
dUPanelTrSolve
int_t dUPanelTrSolve(int_t k, double *BlockLFactor, double *bigV, int_t ldt, Ublock_info_t *, gridinfo_t *, dLUstruct_t *, SuperLUStat_t *, SCT_t *)
superlu_dgemm
int superlu_dgemm(const char *transa, const char *transb, int m, int n, int k, double alpha, double *a, int lda, double *b, int ldb, double beta, double *c, int ldc)
superlu_dscal
int superlu_dscal(const int n, const double alpha, double *x, const int incx)
dDeAllocGlu_3d
int dDeAllocGlu_3d(dLUstruct_t *)
Definition: dutil_dist.c:442
dClone_CompRowLoc_Matrix_dist
void dClone_CompRowLoc_Matrix_dist(SuperMatrix *, SuperMatrix *)
dlangs_dist
double dlangs_dist(char *, SuperMatrix *)
Definition: dlangs_dist.c:62
dLluBufFreeArr
int dLluBufFreeArr(int_t numLA, dLUValSubBuf_t **LUvsbs)
ancestorFactor
int_t ancestorFactor(int_t ilvl, sForest_t *sforest, commRequests_t **comReqss, dscuBufs_t *scuBufs, packLUInfo_t *packLUInfo, msgs_t **msgss, dLUValSubBuf_t **LUvsbs, ddiagFactBufs_t **dFBufs, factStat_t *factStat, factNodelists_t *fNlists, gEtreeInfo_t *gEtreeInfo, superlu_dist_options_t *options, int_t *gIperm_c_supno, int_t ldt, HyP_t *HyP, dLUstruct_t *LUstruct, gridinfo3d_t *grid3d, SuperLUStat_t *stat, double thresh, SCT_t *SCT, int tag_ub, int *info)
pdgsmv_finalize
void pdgsmv_finalize(pdgsmv_comm_t *)
Definition: pdgsmv.c:371
dScalePermstructInit
void dScalePermstructInit(const int_t, const int_t, dScalePermstruct_t *)
Allocate storage in ScalePermstruct.
Definition: dutil_dist.c:398
dlaqgs_dist
void dlaqgs_dist(SuperMatrix *, double *, double *, double, double, double, char *)
Definition: dlaqgs_dist.c:81
ddist_psymbtonum
float ddist_psymbtonum(fact_t, int_t, SuperMatrix *, dScalePermstruct_t *, Pslu_freeable_t *, dLUstruct_t *, gridinfo_t *)
Definition: pdsymbfact_distdata.c:1188
d_c2cpp_GetHWPM
int d_c2cpp_GetHWPM(SuperMatrix *, gridinfo_t *, dScalePermstruct_t *)
Definition: d_c2cpp_GetHWPM.cpp:53
dgemm_
int dgemm_(const char *, const char *, const int *, const int *, const int *, const double *, const double *, const int *, const double *, const int *, const double *, double *, const int *)
dger_
void dger_(const int *, const int *, const double *, const double *, const int *, const double *, const int *, double *, const int *)
superlu_dtrsv
int superlu_dtrsv(char *uplo, char *trans, char *diag, int n, double *a, int lda, double *x, int incx)
dCopy_Dense_Matrix_dist
void dCopy_Dense_Matrix_dist(int_t, int_t, double *, int_t, double *, int_t)
dcreate_matrix_dat
int dcreate_matrix_dat(SuperMatrix *, int, double **, int *, double **, int *, FILE *, gridinfo_t *)
dBcast_LPanel
int_t dBcast_LPanel(int_t k, int_t k0, int_t *lsub, double *lusup, gridinfo_t *, int *msgcnt, int **ToSendR, int_t *xsup, SCT_t *, int)
dLPanelTrSolve
int_t dLPanelTrSolve(int_t k, int_t *factored_L, double *BlockUFactor, gridinfo_t *, dLUstruct_t *)
dreadtriple_dist
void dreadtriple_dist(FILE *, int_t *, int_t *, int_t *, double **, int_t **, int_t **)
Definition: dreadtriple.c:35
pdgstrf3d
int_t pdgstrf3d(superlu_dist_options_t *, int m, int n, double anorm, trf3Dpartition_t *, SCT_t *, dLUstruct_t *, gridinfo3d_t *, SuperLUStat_t *, int *)
Definition: pdgstrf3d.c:121
duser_free_dist
void duser_free_dist(int_t, int_t)
Definition: dmemory_dist.c:49
dZero_CompRowLoc_Matrix_dist
void dZero_CompRowLoc_Matrix_dist(SuperMatrix *)
Sets all entries of a matrix to zero, A_{i,j}=0, for i,j=1,..,n.
Definition: dutil_dist.c:339
dDestroy_trf3Dpartition
void dDestroy_trf3Dpartition(trf3Dpartition_t *trf3Dpartition, gridinfo3d_t *grid3d)
dzRecvUPanel
int_t dzRecvUPanel(int_t k, int_t sender, double alpha, double beta, double *Uval_buf, dLUstruct_t *LUstruct, gridinfo3d_t *grid3d, SCT_t *SCT)
dinit3DLUstruct
int_t dinit3DLUstruct(int_t *myTreeIdxs, int_t *myZeroTrIdxs, int_t *nodeCount, int_t **nodeList, dLUstruct_t *LUstruct, gridinfo3d_t *grid3d)
checkRecvLDiag
int_t checkRecvLDiag(int_t k, commRequests_t *comReqs, gridinfo_t *, SCT_t *)
Definition: treeFactorization.c:422
dgetBigV
double * dgetBigV(int_t, int_t)
dFillRHS_dist
void dFillRHS_dist(char *, int_t, double *, int_t, SuperMatrix *, double *, int_t)
Let rhs[i] = sum of i-th row of A, so the solution vector is all 1's.
Definition: dutil_dist.c:499
dreadtriple_noheader
void dreadtriple_noheader(FILE *, int_t *, int_t *, int_t *, double **, int_t **, int_t **)
Definition: dreadtriple_noheader.c:35
dComputeLevelsets
void dComputeLevelsets(int, int_t, gridinfo_t *, Glu_persist_t *, dLocalLU_t *, int_t *)
Definition: dutil_dist.c:752
dsparseTreeFactor_ASYNC
int_t dsparseTreeFactor_ASYNC(sForest_t *sforest, commRequests_t **comReqss, dscuBufs_t *scuBufs, packLUInfo_t *packLUInfo, msgs_t **msgss, dLUValSubBuf_t **LUvsbs, ddiagFactBufs_t **dFBufs, factStat_t *factStat, factNodelists_t *fNlists, gEtreeInfo_t *gEtreeInfo, superlu_dist_options_t *options, int_t *gIperm_c_supno, int_t ldt, HyP_t *HyP, dLUstruct_t *LUstruct, gridinfo3d_t *grid3d, SuperLUStat_t *stat, double thresh, SCT_t *SCT, int tag_ub, int *info)
doubleCalloc_dist
double * doubleCalloc_dist(int_t)
Definition: dmemory_dist.c:162
dUPanelUpdate
int_t dUPanelUpdate(int_t k, int_t *factored_U, MPI_Request *, double *BlockLFactor, double *bigV, int_t ldt, Ublock_info_t *, gridinfo_t *, dLUstruct_t *, SuperLUStat_t *, SCT_t *)
Printdouble5
void Printdouble5(char *, int_t, double *)
Definition: dutil_dist.c:543
dp3dCollect
int_t dp3dCollect(int_t layer, int_t n, dLUstruct_t *LUstruct, gridinfo3d_t *grid3d)
dPrintLblocks
void dPrintLblocks(int, int_t, gridinfo_t *, Glu_persist_t *, dLocalLU_t *)
Print the blocks in the factored matrix L.
Definition: dutil_dist.c:570
dLluBufInitArr
dLUValSubBuf_t ** dLluBufInitArr(int_t numLA, dLUstruct_t *LUstruct)
sp_dgemm_dist
int sp_dgemm_dist(char *, int, double, SuperMatrix *, double *, int, double, double *, int)
Definition: dsp_blas3_dist.c:124
pdlaqgs
void pdlaqgs(SuperMatrix *, double *, double *, double, double, double, char *)
Definition: pdlaqgs.c:85
dblock_gemm_scatterBottomRight
int_t dblock_gemm_scatterBottomRight(int_t lb, int_t j, double *bigV, int_t knsupc, int_t klst, int_t *lsub, int_t *usub, int_t ldt, int *indirect, int *indirect2, HyP_t *HyP, dLUstruct_t *, gridinfo_t *, SCT_t *SCT, SuperLUStat_t *)
dIrecv_UPanel
int_t dIrecv_UPanel(int_t k, int_t k0, int_t *Usub_buf, double *, dLocalLU_t *, gridinfo_t *, MPI_Request *, int)
dLPanelUpdate
int_t dLPanelUpdate(int_t k, int_t *IrecvPlcd_D, int_t *factored_L, MPI_Request *, double *BlockUFactor, gridinfo_t *, dLUstruct_t *, SCT_t *)
dCopy_CompCol_Matrix_dist
void dCopy_CompCol_Matrix_dist(SuperMatrix *, SuperMatrix *)
dfill_dist
void dfill_dist(double *, int_t, double)
Fills a double precision array with a given value.
Definition: dutil_dist.c:512
pdgssvx
void pdgssvx(superlu_dist_options_t *, SuperMatrix *, dScalePermstruct_t *, double *, int, int, gridinfo_t *, dLUstruct_t *, dSOLVEstruct_t *, double *, SuperLUStat_t *, int *)
NBUFFERS
#define NBUFFERS
Definition: superlu_defs.h:194
int_t
int int_t
Definition: superlu_defs.h:114
DiagScale_t
DiagScale_t
Definition: superlu_enum_consts.h:35
fact_t
fact_t
Definition: superlu_enum_consts.h:30
Mtype_t
Mtype_t
Definition: supermatrix.h:42
Dtype_t
Dtype_t
Definition: supermatrix.h:35
Stype_t
Stype_t
Definition: supermatrix.h:22
pdconvertU
void pdconvertU(superlu_dist_options_t *options, gridinfo_t *grid, dLUstruct_t *LUstruct, SuperLUStat_t *stat, int_t n)
Definition: pdgssvx.c:1639
superlupara_mod::lsub
integer, parameter, public lsub
Definition: superlupara.f90:35
superlupara_mod::trans
integer, parameter, public trans
Definition: superlupara.f90:35
superlupara_mod::factored
integer, parameter, public factored
Definition: superlupara.f90:35
superlupara_mod::lusup
integer, parameter, public lusup
Definition: superlupara.f90:35
superlupara_mod::usub
integer, parameter, public usub
Definition: superlupara.f90:35
C_Tree
Definition: superlu_defs.h:1142
Glu_freeable_t
Definition: superlu_defs.h:490
Glu_persist_t
Definition: superlu_defs.h:435
HyP_t
Definition: superlu_ddefs.h:329
NRformat_loc3d
Definition: supermatrix.h:194
NRformat_loc
Definition: supermatrix.h:176
Pslu_freeable_t
Definition: psymbfact.h:57
Remain_info_t
Definition: superlu_defs.h:770
SCT_t
Definition: util_dist.h:172
SuperLUStat_t
Definition: util_dist.h:95
SuperMatrix
Definition: supermatrix.h:54
Ublock_info_t
Definition: superlu_defs.h:760
Ucb_indptr_t
Definition: superlu_defs.h:752
commRequests_t
Definition: superlu_defs.h:924
dLUValSubBuf_t
Definition: superlu_ddefs.h:357
dLUstruct_t
Definition: superlu_ddefs.h:254
dLocalLU_t
Definition: superlu_ddefs.h:97
dLocalLU_t::Unzval_bc_offset
int64_t * Unzval_bc_offset
Definition: superlu_ddefs.h:132
dLocalLU_t::Ucolind_bc_dat
int_t * Ucolind_bc_dat
Definition: superlu_ddefs.h:126
dLocalLU_t::Uindval_loc_bc_cnt
int64_t Uindval_loc_bc_cnt
Definition: superlu_ddefs.h:138
dLocalLU_t::Unzval_bc_cnt
int64_t Unzval_bc_cnt
Definition: superlu_ddefs.h:133
dLocalLU_t::Ucolind_bc_ptr
int_t ** Ucolind_bc_ptr
Definition: superlu_ddefs.h:125
dLocalLU_t::Ucolind_bc_offset
int64_t * Ucolind_bc_offset
Definition: superlu_ddefs.h:127
dLocalLU_t::Uindval_loc_bc_offset
int64_t * Uindval_loc_bc_offset
Definition: superlu_ddefs.h:137
dLocalLU_t::Ucolind_bc_cnt
int64_t Ucolind_bc_cnt
Definition: superlu_ddefs.h:128
dLocalLU_t::Unzval_bc_ptr
double ** Unzval_bc_ptr
Definition: superlu_ddefs.h:130
dLocalLU_t::Uindval_loc_bc_ptr
int_t ** Uindval_loc_bc_ptr
Definition: superlu_ddefs.h:135
dLocalLU_t::Unzval_bc_dat
double * Unzval_bc_dat
Definition: superlu_ddefs.h:131
dLocalLU_t::Uindval_loc_bc_dat
int_t * Uindval_loc_bc_dat
Definition: superlu_ddefs.h:136
dSOLVEstruct_t
Definition: superlu_ddefs.h:284
dScalePermstruct_t
Definition: superlu_ddefs.h:76
ddiagFactBufs_t
Definition: superlu_ddefs.h:391
dscuBufs_t
Definition: superlu_ddefs.h:385
dtrf3Dpartition_t
Definition: superlu_ddefs.h:406
dtrf3Dpartition_t::myNodeCount
int_t * myNodeCount
Definition: superlu_ddefs.h:409
dtrf3Dpartition_t::sForests
sForest_t ** sForests
Definition: superlu_ddefs.h:413
dtrf3Dpartition_t::supernode2treeMap
int_t * supernode2treeMap
Definition: superlu_ddefs.h:414
dtrf3Dpartition_t::myZeroTrIdxs
int_t * myZeroTrIdxs
Definition: superlu_ddefs.h:411
dtrf3Dpartition_t::treePerm
int_t ** treePerm
Definition: superlu_ddefs.h:412
dtrf3Dpartition_t::LUvsb
dLUValSubBuf_t * LUvsb
Definition: superlu_ddefs.h:415
dtrf3Dpartition_t::gEtreeInfo
gEtreeInfo_t gEtreeInfo
Definition: superlu_ddefs.h:407
dtrf3Dpartition_t::myTreeIdxs
int_t * myTreeIdxs
Definition: superlu_ddefs.h:410
dtrf3Dpartition_t::iperm_c_supno
int_t * iperm_c_supno
Definition: superlu_ddefs.h:408
factNodelists_t
Definition: superlu_defs.h:937
factStat_t
Definition: superlu_defs.h:839
gEtreeInfo_t
Definition: superlu_defs.h:890
gridinfo3d_t
Definition: superlu_defs.h:398
gridinfo_t
Definition: superlu_defs.h:388
lPanelInfo_t
Definition: superlu_ddefs.h:317
msgs_t
Definition: superlu_defs.h:947
packLUInfo_t
Definition: superlu_ddefs.h:397
pdgsmv_comm_t
Definition: superlu_ddefs.h:263
pdgsmv_comm_t::val_torecv
void * val_torecv
Definition: superlu_ddefs.h:308
pdgsmv_comm_t::val_tosend
void * val_tosend
Definition: superlu_ddefs.h:307
pxgstrs_comm_t
Definition: superlu_defs.h:551
sForest_t
Definition: superlu_defs.h:901
superlu_dist_mem_usage_t
Definition: superlu_defs.h:744
superlu_dist_options_t
Definition: superlu_defs.h:712
treeTopoInfo_t
Definition: superlu_defs.h:882
uPanelInfo_t
Definition: superlu_ddefs.h:308
MAX_LOOKAHEADS
#define MAX_LOOKAHEADS
Definition: superlu_ddefs.h:96
superlu_defs.h
Definitions which are precision-neutral.