 |
|
Пример. Функция MPI_Pack/MPI_Unpack
#include <stdio.h>
#include <mpi.h>
void main( int argc, char *argv[] ) {
int ii, iii=2;
double jj, jjj=5.3;
int position, count;
char buff_pack[100], buff_unpack[100];
MPI_Status stat;
if (rank == 0) {
position = 0;
MPI_Pack(&iii,1,MPI_INT,buff_pack,100, &position, MPI_COMM_WORLD);
MPI_Pack(&jjj,1,MPI_DOUBLE,buff_pack,100,&position, MPI_COMM_WORLD);
MPI_Send(&buff_pack[0], position, MPI_PACKED, 1, 0, MPI_COMM_WORLD);
}
else {
MPI_Probe(0, 0, MPI_COMM_WORLD, &stat);
MPI_Get_count(&stat, MPI_PACKED, &count);
MPI_Recv( &buff_unpack[0], count, MPI_PACKED, 0, 0, MPI_COMM_WORLD, &stat);
position = 0;
MPI_Unpack(&buff_unpack[0],count,&position, &ii,1, MPI_INT, MPI_COMM_WORLD);
MPI_Unpack(&buff_unpack[0],count,&position, &jj,1, MPI_DOUBLE,
MPI_COMM_WORLD);
printf("ii=%d jj=%f\n",ii,jj);
}
MPI_Finalize();
}
MPI содержит механизм, который позволяет передавать данные различной структуры и размера, выполняя передачу данные непосредственно из памяти, т.е. без дополнительного копирования. Аналогично тому, как создаются производные типы данных в Си, производные типы данных MPI, конструируются из базовых типов с использованием специальных конструкторов.
В общем случае описываемые значения не обязательно непрерывно располагаются в памяти. Каждый тип данных MPI имеет карту типа (type map), состоящую изпоследовательности пар значений, первое задает базовых тип данных, второе - смещение адреса в памяти относительно начала базового адреса, т.е.
Часть карты типа с указанием только базовых типов называется сигнатурой типа:
Сигнатура типа описывает, какие базовые типы данных образуют производный тип данных MPI. Смещения карты типа определяют, где находятся значения данных.
Например, пусть в сообщение должны входить значения переменных:
double a; /* адрес 24 */
double b; /* адрес 40 */
int n; /* адрес 48 */
Тогда производный тип для описания таких данных должен иметь карту типа следующего вида: { (MPI_DOUBLE,0), (MPI_DOUBLE,16), (MPI_INT,24) }
Для характеристики типа данных в MPI используется следующий ряд понятий:
- нижняя граница типа 
- верхняя граница типа 
- протяженность типа extent(TypeMap) = ub(TypeMap)-lb(TypeMap).
Нижняя граница – определяет смещение для первого байта значений типа данных.
Верхняя граница – определяет смещение для байта, расположенного за последним элементом типа данных. При этом величина смещения для верхней границы может быть округлена вверх с учетом требований выравнивания адресов.
Понятие протяженности отличается от понятия размер типа. Протяженность – это размер памяти в байтах, который нужно отводить для одного элемента производного типа. Размер типа данных - это число байтов, которые занимают в памяти данные (разность между адресами последнего и первого байтов данных). Отсюда - различие в значениях протяженности и размера зависит от величины округления при выравнивании адресов. Так, в рассматриваемом примере размер типа равен 28, а протяженность – 32 (предполагается, что выравнивание выставлено на 8 байт).
Итак, создание общего типа данных определяется:
- последовательностью базовых типов данных,
- последовательностью смещений байтов.
Непрерывный.
Конструирует новый тип данных, размножением старого типа данных в непрерывную область памяти.
MPI_Type_contiguous(count, oldtype, newtype)
IN countчисло повторений (неотрицательное целое)
IN oldtypeстарый тип данных (дескриптор)
OUT newtypeновый тип данных (дескриптор)
newtype это тип данных, полученный размножением count копий oldtype.
Сконструированный объект типа данных должен быть зафиксирован перед тем, как он может быть использован в обменах. Функция, которая фиксирует новый тип данных имеет следующий вид:
MPI_Type_commit(datatype)
| INOUT | Datatype | тип данных, который фиксируется (handle) |
После того, как тип данных был объявлен, он может быть использован для выполнения обменов и других операций с типами данных. На его основе может быть создан уже новый тип данных.