이것은 type mapMPI에서 중요하지만 혼란스러운 개념입니다. 유형 맵을 표시하거나 인쇄하는 루틴을 원합니다.
예를 들어 (MPI-3 표준에서 가져옴),
MPI_TYPE_CREATE_RESIZED(MPI_INT, -3, 9, type1)
typemap 결과
{(lb_marker, -3), (int, 0), (ub_marker, 6)}.
해당 유형을 다시 사용하십시오.
MPI_TYPE_CONTIGUOUS(2, type1, type2)
그리고 typemap은
{(lb_marker, -3), (int, 0), (int,9), (ub_marker, 15)}
해당 typemap을 자동으로 표시하는 방법을 원합니다.
물론 하나는 사용할 수 MPI_Type_get_contents및 MPI_Type_get_envelope재귀 유형 내장 타격까지 내려. 이것은 오히려 큰 고통이며 나는 이것을 위해 이것을 할 수있는 어떤 도구가 20 년 동안 존재할 것이라고 생각했을 것입니다.
유망하지만 제대로 작동하지 않는 일부 도구 :
-
나는 ~ 2001 MPImap을 발견했다 여기 . 첫째, 최신 Tcl / TK 용으로 업데이트해야하며 일부 메모리 오류를 해결하기 위해 패치를 적용해야합니다. 응답하지 않는 GUI가 표시됩니다. 대신 런타임에 호출 할 수있는 라이브러리 / 루틴을 찾고 있습니다.
-
MPIDU_Datatype_deubgMPICH 특정 내부 유형 덤프 루틴입니다. 유형 맵을 표시하지 않습니다 (데이터 루프 표현을 표시하고 다시 닫음). -
한때 MPI 유형 맵을 표시하는 기능을 나열하는 XMPI라는 디버거가있었습니다. 이 디버거는 LAM-MPI에 특정한 것으로 보이며 get_contents / get_envelope를 사용하지 않습니다.
답변
마찬가지로 롭 레이 텀이 말했다, 좋은 기존 솔루션이 없다. tim 이 제공 한 링크의 도움 으로 Github에서 사용할 수있는이 기능을 만들었습니다 . 연속 + 크기 조정 테스트 ( 여기 )에 대한 예제를 사용 했으며 출력은
contiguous + resize
"(LB, -3), (MPI_INT, 0), (MPI_INT, 9), (UB, 15)"
이 기능을 사용하면 printMapDatatype(mydatatype). 나는 이것이 당신이 찾고 있던 것이기를 바랍니다.
다음과 같은 경우에 기능이 있습니다.
MPI_Aint printdatatype( MPI_Datatype datatype, MPI_Aint prevExtentTot ) {
int *array_of_ints;
MPI_Aint *array_of_adds;
MPI_Datatype *array_of_dtypes;
int num_ints, num_adds, num_dtypes, combiner;
int i, j;
MPI_Type_get_envelope( datatype, &num_ints, &num_adds, &num_dtypes, &combiner );
array_of_ints = (int *) malloc( num_ints * sizeof(int) );
array_of_adds = (MPI_Aint *) malloc( num_adds * sizeof(MPI_Aint) );
array_of_dtypes = (MPI_Datatype *) malloc( num_dtypes * sizeof(MPI_Datatype) );
MPI_Aint extent, subExtent;
MPI_Type_extent(datatype, &extent);
switch (combiner) {
case MPI_COMBINER_NAMED:
// To print the specific type, we can match against the predefined forms.
if (datatype == MPI_BYTE) printf( "(MPI_BYTE, %ld)", prevExtentTot);
else if (datatype == MPI_LB) printf( "(MPI_LB, %ld)", prevExtentTot);
else if (datatype == MPI_PACKED) printf( "(MPI_PACKED, %ld)", prevExtentTot);
else if (datatype == MPI_UB) printf( "(MPI_UB, %ld)", prevExtentTot);
else if (datatype == MPI_CHAR) printf( "(MPI_CHAR, %ld)", prevExtentTot);
else if (datatype == MPI_DOUBLE) printf( "(MPI_DOUBLE, %ld)", prevExtentTot);
else if (datatype == MPI_FLOAT) printf( "(MPI_FLOAT, %ld)", prevExtentTot);
else if (datatype == MPI_INT) printf( "(MPI_INT, %ld)", prevExtentTot );
else if (datatype == MPI_LONG) printf( "(MPI_LONG, %ld)", prevExtentTot);
else if (datatype == MPI_LONG_DOUBLE) printf( "(MPI_LONG_DOUBLE, %ld)", prevExtentTot);
else if (datatype == MPI_LONG_LONG) printf( "(MPI_LONG_LONG, %ld)", prevExtentTot);
else if (datatype == MPI_LONG_LONG_INT) printf( "(MPI_LONG_LONG_INT, %ld)", prevExtentTot);
else if (datatype == MPI_SHORT) printf( "(MPI_SHORT, %ld)", prevExtentTot);
else if (datatype == MPI_SIGNED_CHAR) printf( "(MPI_SIGNED_CHAR, %ld)", prevExtentTot);
else if (datatype == MPI_UNSIGNED) printf( "(MPI_UNSIGNED, %ld)", prevExtentTot);
else if (datatype == MPI_UNSIGNED_CHAR) printf( "(MPI_UNSIGNED_CHAR, %ld)", prevExtentTot);
else if (datatype == MPI_UNSIGNED_LONG) printf( "(MPI_UNSIGNED_LONG, %ld)", prevExtentTot);
else if (datatype == MPI_UNSIGNED_LONG_LONG)printf( "(MPI_UNSIGNED_LONG_LONG, %ld)", prevExtentTot);
else if (datatype == MPI_UNSIGNED_SHORT) printf( "(MPI_UNSIGNED_SHORT, %ld)", prevExtentTot);
else if (datatype == MPI_WCHAR) printf( "(MPI_WCHAR, %ld)", prevExtentTot);
free( array_of_ints );
free( array_of_adds );
free( array_of_dtypes );
return prevExtentTot;
break;
case MPI_COMBINER_DUP:
MPI_Type_get_contents( datatype, num_ints, num_adds, num_dtypes, array_of_ints, array_of_adds, array_of_dtypes );
printdatatype( array_of_dtypes[0], prevExtentTot);
printf(", \n");
break;
case MPI_COMBINER_CONTIGUOUS:
MPI_Type_get_contents( datatype, num_ints, num_adds, num_dtypes, array_of_ints, array_of_adds, array_of_dtypes );
MPI_Type_extent(array_of_dtypes[0], &subExtent); // no need to do in loop because same type
for (i=0; i < array_of_ints[0]; i++) {
prevExtentTot = printdatatype( array_of_dtypes[0], prevExtentTot);
prevExtentTot += subExtent;
printf(", ");
}
break;
case MPI_COMBINER_VECTOR:
MPI_Type_get_contents( datatype, num_ints, num_adds, num_dtypes, array_of_ints, array_of_adds, array_of_dtypes );
MPI_Type_extent(array_of_dtypes[0], &subExtent); // no need to do in loop because same type
printf("[");
for (i = 0; i < array_of_ints[0]; i++) { //count
printf( "BL : %d - ", array_of_ints[1]);
for (j = 0; j < array_of_ints[2]; j++) { // stride
if (j < array_of_ints[1]) { // if in blocklength
prevExtentTot = printdatatype( array_of_dtypes[0], prevExtentTot);
printf(", ");
}
prevExtentTot += subExtent;
}
}
printf("], ");
break;
case MPI_COMBINER_HVECTOR:
case MPI_COMBINER_HVECTOR_INTEGER:{
MPI_Aint backupPrevExtent = prevExtentTot;
MPI_Type_get_contents( datatype, num_ints, num_adds, num_dtypes, array_of_ints, array_of_adds, array_of_dtypes );
MPI_Type_extent(array_of_dtypes[0], &subExtent); // no need to do in loop because same type
printf("[");
for (i = 0; i < array_of_ints[0]; i++) { //count
printf( "BL : %d - ", array_of_ints[1]);
for (j = 0; j < array_of_ints[1]; j++) { // blocklength
prevExtentTot = printdatatype( array_of_dtypes[0], prevExtentTot);
printf(", ");
prevExtentTot += subExtent;
}
prevExtentTot = backupPrevExtent + array_of_adds[0]; // + stride un byte
}
printf("], ");
break;
}
case MPI_COMBINER_INDEXED:{
MPI_Aint tmpPrevExtent;
int count, blocklength;
MPI_Type_get_contents( datatype, num_ints, num_adds, num_dtypes, array_of_ints, array_of_adds, array_of_dtypes );
MPI_Type_extent(array_of_dtypes[0], &subExtent); // no need to do in loop because same type
printf("<");
count = array_of_ints[0];
for (i = 0; i < count; i++) { // count
blocklength = array_of_ints[i + 1]; // array of blocklength
tmpPrevExtent = prevExtentTot;
tmpPrevExtent += array_of_ints[count + 1 + i] * subExtent; // + displacement * size of block
printf( "BL : %d - ", blocklength);
for (j = 0; j < blocklength; j++) { // blocklength
tmpPrevExtent = printdatatype( array_of_dtypes[0], tmpPrevExtent);
printf(", ");
tmpPrevExtent += subExtent;
}
}
printf(">, ");
prevExtentTot = tmpPrevExtent;
break;
}
case MPI_COMBINER_HINDEXED:
case MPI_COMBINER_HINDEXED_INTEGER:{
MPI_Aint tmpPrevExtent;
int count, blocklength;
MPI_Type_get_contents( datatype, num_ints, num_adds, num_dtypes, array_of_ints, array_of_adds, array_of_dtypes );
MPI_Type_extent(array_of_dtypes[0], &subExtent); // no need to do in loop because same type
printf("<");
count = array_of_ints[0];
for (i = 0; i < count; i++) { // count
blocklength = array_of_ints[i + 1]; // array of blocklength
tmpPrevExtent = prevExtentTot;
tmpPrevExtent += array_of_adds[i]; // + displacement in byte
printf( "BL : %d - ", blocklength);
for (j = 0; j < blocklength; j++) {
tmpPrevExtent = printdatatype( array_of_dtypes[0], tmpPrevExtent);
printf(", ");
tmpPrevExtent += subExtent;
}
}
printf(">, ");
prevExtentTot = tmpPrevExtent;
break;
}
case MPI_COMBINER_INDEXED_BLOCK:{
MPI_Aint tmpPrevExtent;
int count;
MPI_Type_get_contents( datatype, num_ints, num_adds, num_dtypes, array_of_ints, array_of_adds, array_of_dtypes );
MPI_Type_extent(array_of_dtypes[0], &subExtent); // no need to do in loop because same type
printf("<");
count = array_of_ints[0];
for (i = 0; i < count; i++) { // count
tmpPrevExtent = prevExtentTot;
tmpPrevExtent += array_of_ints[i + 2] * subExtent; // + displacement * size of block
printf( "BL : %d - ", array_of_ints[i + 1]);
for (j = 0; j < array_of_ints[1]; j++) { // blocklength
tmpPrevExtent = printdatatype( array_of_dtypes[0], tmpPrevExtent);
printf(", ");
tmpPrevExtent += subExtent;
}
}
printf(">, ");
prevExtentTot = tmpPrevExtent;
break;
}
case MPI_COMBINER_STRUCT:
case MPI_COMBINER_STRUCT_INTEGER:{
MPI_Aint tmpPrevExtent;
MPI_Type_get_contents( datatype, num_ints, num_adds, num_dtypes, array_of_ints, array_of_adds, array_of_dtypes );
printf( "{");
for (i = 0; i < array_of_ints[0]; i++) { // count
tmpPrevExtent = prevExtentTot + array_of_adds[i]; // origin + displacement
printf( "BL : %d - ", array_of_ints[i + 1]);
tmpPrevExtent = printdatatype( array_of_dtypes[i], tmpPrevExtent);
tmpPrevExtent += subExtent;
printf(", ");
}
printf("}, ");
prevExtentTot = tmpPrevExtent;
break;
}
case MPI_COMBINER_SUBARRAY:
// I don't know what is interresting to display here...
printf("... subarray not handled ...");
break;
case MPI_COMBINER_DARRAY:
// Same
printf("... darray not handled ...");
break;
case MPI_COMBINER_RESIZED:
MPI_Type_get_contents( datatype, num_ints, num_adds, num_dtypes, array_of_ints, array_of_adds, array_of_dtypes );
prevExtentTot = printdatatype( array_of_dtypes[0], prevExtentTot);
printf(", \n");
break;
default:
printf( "Unrecognized combiner type\n" );
}
free( array_of_ints );
free( array_of_adds );
free( array_of_dtypes );
return prevExtentTot;
}
void printMapDatatype(MPI_Datatype datatype) {
MPI_Aint lb, ub;
MPI_Type_lb(datatype, &lb);
MPI_Type_ub(datatype, &ub);
printf("\"(LB, %ld), ", lb);
printdatatype(datatype, 0);
printf("(UB, %ld)\"\n", ub);
}
