1*b75c6efcSBarry Smith /* 2*b75c6efcSBarry Smith 3*b75c6efcSBarry Smith This is the equivalent of MATLAB's fwrite() only on sockets instead of 4*b75c6efcSBarry Smith binary files. 5*b75c6efcSBarry Smith */ 6*b75c6efcSBarry Smith 7*b75c6efcSBarry Smith #include <petscsys.h> 8*b75c6efcSBarry Smith #include <../src/sys/classes/viewer/impls/socket/socket.h> 9*b75c6efcSBarry Smith #include <mex.h> 10*b75c6efcSBarry Smith 11*b75c6efcSBarry Smith PetscErrorCode PetscBinaryWrite(int, const void *p, int, PetscDataType); 12*b75c6efcSBarry Smith 13*b75c6efcSBarry Smith #define PETSC_MEX_ERROR(a) \ 14*b75c6efcSBarry Smith { \ 15*b75c6efcSBarry Smith fprintf(stdout, "sread: %s \n", a); \ 16*b75c6efcSBarry Smith return; \ 17*b75c6efcSBarry Smith } 18*b75c6efcSBarry Smith /*-----------------------------------------------------------------*/ 19*b75c6efcSBarry Smith /* */ 20*b75c6efcSBarry Smith /*-----------------------------------------------------------------*/ 21*b75c6efcSBarry Smith PETSC_EXTERN void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]) 22*b75c6efcSBarry Smith { 23*b75c6efcSBarry Smith int i, fd, cnt, dt; 24*b75c6efcSBarry Smith PetscErrorCode ierr; 25*b75c6efcSBarry Smith 26*b75c6efcSBarry Smith /* check output parameters */ 27*b75c6efcSBarry Smith if (nrhs != 3) PETSC_MEX_ERROR("Receive requires three input arguments."); 28*b75c6efcSBarry Smith fd = (int)mxGetScalar(prhs[0]); 29*b75c6efcSBarry Smith cnt = mxGetNumberOfElements(prhs[1]); 30*b75c6efcSBarry Smith dt = (PetscDataType)mxGetScalar(prhs[2]); 31*b75c6efcSBarry Smith 32*b75c6efcSBarry Smith if (dt == PETSC_DOUBLE) { 33*b75c6efcSBarry Smith ierr = PetscBinaryWrite(fd, mxGetPr(prhs[1]), cnt, (PetscDataType)dt); 34*b75c6efcSBarry Smith if (ierr) PETSC_MEX_ERROR("Unable to send double items."); 35*b75c6efcSBarry Smith } else if (dt == PETSC_INT) { 36*b75c6efcSBarry Smith int *tmp = (int *)mxMalloc((cnt + 5) * sizeof(int)); 37*b75c6efcSBarry Smith double *t = mxGetPr(prhs[1]); 38*b75c6efcSBarry Smith for (i = 0; i < cnt; i++) tmp[i] = (int)t[i]; 39*b75c6efcSBarry Smith ierr = PetscBinaryWrite(fd, tmp, cnt, (PetscDataType)dt); 40*b75c6efcSBarry Smith if (ierr) PETSC_MEX_ERROR("Unable to send int items."); 41*b75c6efcSBarry Smith mxFree(tmp); 42*b75c6efcSBarry Smith } else if (dt == PETSC_CHAR) { 43*b75c6efcSBarry Smith char *tmp = (char *)mxMalloc((cnt + 5) * sizeof(char)); 44*b75c6efcSBarry Smith mxGetNChars(prhs[1], tmp, cnt + 1); 45*b75c6efcSBarry Smith ierr = PetscBinaryWrite(fd, tmp, cnt, (PetscDataType)dt); 46*b75c6efcSBarry Smith if (ierr) PETSC_MEX_ERROR("Unable to send char items."); 47*b75c6efcSBarry Smith mxFree(tmp); 48*b75c6efcSBarry Smith } else PETSC_MEX_ERROR("Unknown datatype."); 49*b75c6efcSBarry Smith return; 50*b75c6efcSBarry Smith } 51*b75c6efcSBarry Smith 52*b75c6efcSBarry Smith int main(int argc, char **argv) 53*b75c6efcSBarry Smith { 54*b75c6efcSBarry Smith return 0; 55*b75c6efcSBarry Smith } 56