128-bit AVX2 SIMD for AMD Ryzen
[gromacs.git] / src / gromacs / imd / imdsocket.cpp
blob06008c35bb87c9922b34cca31953eeaa46c4eeae
1 /*
2 * This file is part of the GROMACS molecular simulation package.
4 * Copyright (c) 2014,2015,2016,2017, by the GROMACS development team, led by
5 * Mark Abraham, David van der Spoel, Berk Hess, and Erik Lindahl,
6 * and including many others, as listed in the AUTHORS file in the
7 * top-level source directory and at http://www.gromacs.org.
9 * GROMACS is free software; you can redistribute it and/or
10 * modify it under the terms of the GNU Lesser General Public License
11 * as published by the Free Software Foundation; either version 2.1
12 * of the License, or (at your option) any later version.
14 * GROMACS is distributed in the hope that it will be useful,
15 * but WITHOUT ANY WARRANTY; without even the implied warranty of
16 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
17 * Lesser General Public License for more details.
19 * You should have received a copy of the GNU Lesser General Public
20 * License along with GROMACS; if not, see
21 * http://www.gnu.org/licenses, or write to the Free Software Foundation,
22 * Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
24 * If you want to redistribute modifications to GROMACS, please
25 * consider that scientific software is very special. Version
26 * control is crucial - bugs must be traceable. We will be happy to
27 * consider code for inclusion in the official distribution, but
28 * derived work must not be called official GROMACS. Details are found
29 * in the README & COPYING files - if they are missing, get the
30 * official version at http://www.gromacs.org.
32 * To help us fund GROMACS development, we humbly ask that you cite
33 * the research papers on the package. Check out http://www.gromacs.org.
36 /*! \internal \file
38 * \brief
39 * Implements functions of imdsocket.h.
41 * This file re-implements vmdsock.c functions from original IMD API from scratch,
42 * see imdsocket.h for references to the IMD API.
44 * \author Martin Hoefling, Carsten Kutzner <ckutzne@gwdg.de>
46 * \ingroup module_imd
49 #include "gmxpre.h"
51 #include "imdsocket.h"
53 #include "config.h"
55 #include <errno.h>
56 #include <string.h>
58 #include "gromacs/imd/imd.h"
59 #include "gromacs/utility/fatalerror.h"
60 #include "gromacs/utility/smalloc.h"
62 #if GMX_NATIVE_WINDOWS
63 #ifdef GMX_HAVE_WINSOCK
64 /*! \brief Define socklen type on Windows. */
65 typedef int socklen_t;
67 /*! \brief Define a function to initialize winsock. */
68 extern int imdsock_winsockinit()
70 int ret = -1;
73 WSADATA wsd;
75 /* We use winsock 1.1 compatibility for now. Though I guess no one will try on Windows 95. */
76 ret = WSAStartup(MAKEWORD(1, 1), &wsd);
77 return ret;
79 #endif
80 #else
81 /* On UNIX, we can use nice errors from errno.h */
82 #include <unistd.h>
83 #ifdef GMX_IMD
84 #include <time.h>
86 #include <sys/select.h>
87 #include <sys/time.h>
88 #endif
89 #endif
92 /*! \brief Simple error handling. */
93 #if GMX_NATIVE_WINDOWS
94 #define ERR_ARGS __FILE__, __LINE__, NULL
95 #else
96 #define ERR_ARGS __FILE__, __LINE__, strerror(errno)
97 #endif
100 /*! \brief Currently only 1 client connection is supported. */
101 #define MAXIMDCONNECTIONS 1
104 /*! \brief Print a nice error message on UNIX systems, using errno.h. */
105 static void print_IMD_error(const char *file, int line, char *msg)
107 fprintf(stderr, "%s Error in file %s on line %d.\n", IMDstr, file, line);
109 if (nullptr != msg)
111 fprintf(stderr, "%s\n", msg);
115 /*! \brief Byte swap in case we are little-endian */
116 static uint16_t gmx_htons(uint16_t src)
118 uint16_t num = 1;
120 if (*(char *)&num == 1)
122 return src;
124 else
126 uint16_t dest = 0;
128 dest |= (src & 0x0000FF00) >> 8;
129 dest |= (src & 0x000000FF) << 8;
131 return dest;
135 /*! \brief Byte-unswap 16 bit word in case we are little-endian */
136 static uint16_t gmx_ntohs(uint16_t src)
138 return gmx_htons(src);
142 extern IMDSocket* imdsock_create()
144 IMDSocket *sock = nullptr;
147 #ifdef GMX_IMD
148 snew(sock, 1);
149 /* Try to create socket: */
150 if ((sock->sockfd = socket(PF_INET, SOCK_STREAM, 0)) == -1)
152 print_IMD_error(ERR_ARGS);
153 sfree(sock);
155 return nullptr;
157 else
158 #endif
160 return sock;
165 extern int imdsock_bind(IMDSocket *sock, int port)
167 int ret;
170 #ifdef GMX_IMD
171 memset(&(sock->address), 0, sizeof(sock->address));
172 sock->address.sin_family = PF_INET;
173 sock->address.sin_port = gmx_htons(port);
175 /* Try to bind to address and port ...*/
176 ret = bind(sock->sockfd, (struct sockaddr *) &sock->address, sizeof(sock->address));
177 #else
178 ret = -1;
179 #endif
181 if (ret)
183 print_IMD_error(ERR_ARGS);
186 return ret;
190 extern int imd_sock_listen(IMDSocket *sock)
192 int ret;
195 #ifdef GMX_IMD
196 /* Try to set to listening state */
197 ret = listen(sock->sockfd, MAXIMDCONNECTIONS);
198 #else
199 ret = -1;
200 #endif
202 if (ret)
204 print_IMD_error(ERR_ARGS);
207 return ret;
211 extern IMDSocket* imdsock_accept(IMDSocket *sock)
213 int ret;
215 #ifdef GMX_IMD
216 socklen_t length;
219 length = sizeof(sock->address);
220 ret = accept(sock->sockfd, (struct sockaddr *) &sock->address, &length);
222 /* successful, redirect to distinct clientsocket */
223 if (ret >= 0)
225 IMDSocket *newsock;
227 snew(newsock, 1);
228 newsock->address = sock->address;
229 newsock->sockfd = ret;
231 return newsock;
233 else
234 #endif
236 print_IMD_error(ERR_ARGS);
238 return nullptr;
243 extern int imdsock_getport(IMDSocket *sock, int *port)
245 int ret;
246 #ifdef GMX_IMD
247 socklen_t len;
250 len = sizeof(struct sockaddr_in);
251 ret = getsockname(sock->sockfd, (struct sockaddr *) &(sock->address), &len);
252 if (ret)
254 fprintf(stderr, "%s getsockname failed with error %d.\n", IMDstr, ret);
255 print_IMD_error(ERR_ARGS);
257 else
259 *port = gmx_ntohs(sock->address.sin_port);
261 #else
262 gmx_incons("imdsock_getport called without IMD support.");
263 ret = -1;
264 #endif
266 return ret;
270 extern int imdsock_write(IMDSocket *sock, const char *buffer, int length)
272 /* No read and write on windows, we have to use send and recv instead... */
273 #if GMX_NATIVE_WINDOWS
274 #ifdef GMX_HAVE_WINSOCK
275 return send(sock->sockfd, (const char *) buffer, length, NOFLAGS);
276 #else
277 return -1;
278 #endif
279 #else
280 return write(sock->sockfd, buffer, length);
281 #endif
285 extern int imdsock_read(IMDSocket *sock, char *buffer, int length)
287 /* See above... */
288 #if GMX_NATIVE_WINDOWS
289 #ifdef GMX_HAVE_WINSOCK
290 return recv(sock->sockfd, (char *) buffer, length, NOFLAGS);
291 #else
292 return -1;
293 #endif
294 #else
295 return read(sock->sockfd, buffer, length);
296 #endif
300 extern void imdsock_shutdown(IMDSocket *sock)
302 int ret = -1;
305 /* is the socket already NULL? */
306 if (sock == nullptr)
308 return;
311 #ifdef GMX_IMD
312 /* If not, try to properly shut down. */
313 ret = shutdown(sock->sockfd, 1);
314 #endif
316 if (ret == -1)
318 fprintf(stderr, "%s Failed to shutdown socket. Did the client already disconnect?\n", IMDstr);
319 print_IMD_error(ERR_ARGS);
324 extern int imdsock_destroy(IMDSocket *sock)
326 int ret = -1;
329 if (sock == nullptr)
331 return 1;
334 #if GMX_NATIVE_WINDOWS
335 /* On Windows, this function is called closesocket */
336 #ifdef GMX_HAVE_WINSOCK
337 ret = closesocket(sock->sockfd);
338 #endif
339 #else
340 ret = close(sock->sockfd);
341 #endif
343 if (ret == -1)
345 sfree(sock);
346 print_IMD_error(ERR_ARGS);
348 return 0;
351 return 1;
355 extern int imdsock_tryread(IMDSocket *sock, int timeoutsec, int timeoutusec)
357 int ret = -1;
360 #ifdef GMX_IMD
361 fd_set readfds;
362 /* Create new time structure with sec and usec. */
363 struct timeval *tval;
366 snew(tval, 1);
368 /* clear the set */
369 FD_ZERO(&readfds);
370 /* add the socket to the read set */
371 FD_SET(sock->sockfd, &readfds);
373 /* set the timeout */
374 tval->tv_sec = timeoutsec;
375 tval->tv_usec = timeoutusec;
378 /* check the set for read readiness. */
379 ret = select(sock->sockfd + 1, &readfds, nullptr, nullptr, tval);
380 /* redo on system interrupt */
382 while (ret < 0 && errno == EINTR);
384 sfree(tval);
385 #endif
387 if (ret < 0)
389 print_IMD_error(ERR_ARGS);
392 return ret;