Added Michael Shirts long-range dispersion and repulsion corrections that are correct...
[gromacs.git] / src / gmxlib / repfirst.c
blobe2c463be2101b4762c8dddf2de4f65a111b3760c
1 /*
2 * $Id$
3 *
4 * This source code is part of
5 *
6 * G R O M A C S
7 *
8 * GROningen MAchine for Chemical Simulations
9 *
10 * VERSION 3.2.0
11 * Written by David van der Spoel, Erik Lindahl, Berk Hess, and others.
12 * Copyright (c) 1991-2000, University of Groningen, The Netherlands.
13 * Copyright (c) 2001-2004, The GROMACS development team,
14 * check out http://www.gromacs.org for more information.
16 * This program is free software; you can redistribute it and/or
17 * modify it under the terms of the GNU General Public License
18 * as published by the Free Software Foundation; either version 2
19 * of the License, or (at your option) any later version.
21 * If you want to redistribute modifications, please consider that
22 * scientific software is very special. Version control is crucial -
23 * bugs must be traceable. We will be happy to consider code for
24 * inclusion in the official distribution, but derived work must not
25 * be called official GROMACS. Details are found in the README & COPYING
26 * files - if they are missing, get the official version at www.gromacs.org.
28 * To help us fund GROMACS development, we humbly ask that you cite
29 * the papers on the package - you can find them in the top README file.
31 * For more info, check our website at http://www.gromacs.org
33 * And Hey:
34 * GROningen Mixture of Alchemy and Childrens' Stories
36 /* This file is completely threadsafe - keep it that way! */
38 #include "string2.h"
39 #include "typedefs.h"
40 #include "smalloc.h"
41 #include "macros.h"
42 #include "replace.h"
43 #include "strdb.h"
45 bool replace1(char *string,char *buf,char *search,char *replace)
47 char *ptr,*bufptr;
48 int blen,stringlen,slen,rlen;
49 int i,j,tmp;
51 slen=strlen(search);
52 stringlen=strlen(string);
53 if ((string == NULL) || (slen == 0) || (stringlen == 0))
54 return FALSE;
56 rlen=strlen(replace);
57 sprintf(buf,"%s",string);
59 bufptr=buf;
60 if ((ptr=strstr(bufptr,search)) != NULL) {
61 if (rlen <= slen) {
62 for(i=0; (i<rlen); i++)
63 ptr[i]=replace[i];
64 if (rlen < slen) {
65 while (ptr[i+slen-rlen] != '\0') {
66 ptr[i]=ptr[i+slen-rlen];
67 i++;
69 ptr[i]='\0';
72 else {
73 tmp=strlen(ptr);
74 for(j=tmp; (j>=slen); j--)
75 ptr[rlen-slen+j]=ptr[j];
76 for(i=0; (i<rlen); i++)
77 ptr[i]=replace[i];
79 bufptr=ptr+rlen;
81 return TRUE;
84 return FALSE;
87 void main(int argc,char *argv[])
89 #define MYBUFS 10240
90 int i,nstr;
91 char **str;
92 char **rep;
93 bool *bRep;
94 char buffer[MYBUFS];
95 char newbuf[MYBUFS];
97 nstr=get_strings("index.gmx",&str);
98 snew(bRep,nstr);
99 snew(rep,nstr);
100 for(i=0; (i<nstr); i++) {
101 sprintf(buffer,"\\myindex{%s}",str[i]);
102 rep[i]=strdup(buffer);
105 while(fgets2(buffer,MYBUFS-1,stdin) != NULL) {
106 newbuf[0]='\0';
107 for(i=0; (i<nstr); i++)
108 if (!bRep[i]) {
109 bRep[i]=replace1(buffer,newbuf,str[i],rep[i]);
110 if (bRep[i])
111 strcpy(buffer,newbuf);
113 printf("%s\n",newbuf);