xuMRTConstants readXuMRTConstants( char* fname ) {
xuMRTConstants param;
char* msg;
uint status;
// Relaxation times
uint n;
status = vstring(&msg, "%s/Lambda", fname);
if( !lookUpVectorEntry("properties/macroProperties", msg, ¶m.Lambda, &n) ) {
status = vstring(&msg, "Ill defined entry %s/Lambda", fname);
if(status)
errorMsg( msg );
}
// Sigma
status = vstring(&msg, "%s/sigma", fname);
if( !lookUpScalarEntry("properties/macroProperties", msg, 0, ¶m.sigma) ) {
status = vstring(&msg, "Ill defined entry %s/sigma", fname);
if(status)
errorMsg( msg );
}
// Surface tension
status = vstring(&msg, "%s/SurfaceTension", fname);
char* aux;
if( lookUpStringEntry("properties/macroProperties", msg, &aux, "none") ) {
if( strcmp(aux, "none") == 0 ) {
param.surfaceTension = noXuSurf;
}
else {
if( strcmp(aux, "xuSurfTen") == 0 ) {
param.surfaceTension = xuSurfTen;
// Look up kappa
status = vstring(&msg, "%s/kappa_st", fname);
if( !lookUpScalarEntry("properties/macroProperties", msg, -1, ¶m.kappa_st) ) {
status = vstring(&msg, "Unable to find %s/kappa_st", fname);
errorMsg( msg );
}
}
}
}
else {
param.surfaceTension = noXuSurf;
}
return param;
}
static void
getremcap(char *host)
{
char **p, ***q;
char *bp;
char *rempath;
int stat;
rempath = getenv("REMOTE");
if (rempath != NULL) {
if (*rempath != '/')
/* we have an entry */
cgetset(rempath);
else { /* we have a path */
db_array[1] = rempath;
db_array[2] = _PATH_REMOTE;
}
}
if ((stat = cgetent(&bp, db_array, host)) < 0) {
if (DV ||
(host[0] == '/' && access(DV = host, R_OK | W_OK) == 0)) {
CU = DV;
HO = host;
HW = 1;
DU = 0;
if (!BR)
BR = DEFBR;
FS = BUFSIZ;
return;
}
switch(stat) {
case -1:
warnx("unknown host %s", host);
break;
case -2:
warnx("can't open host description file");
break;
case -3:
warnx("possible reference loop in host description file");
break;
}
exit(3);
}
for (p = capstrings, q = caps; *p != NULL; p++, q++)
if (**q == NULL)
cgetstr(bp, *p, *q);
if (!BR && (cgetnum(bp, "br", &BR) == -1))
BR = DEFBR;
if (cgetnum(bp, "fs", &FS) == -1)
FS = BUFSIZ;
if (DU < 0)
DU = 0;
else
DU = cgetflag("du");
if (DV == NULL) {
fprintf(stderr, "%s: missing device spec\n", host);
exit(3);
}
if (DU && CU == NULL)
CU = DV;
if (DU && PN == NULL) {
fprintf(stderr, "%s: missing phone number\n", host);
exit(3);
}
HD = cgetflag("hd");
/*
* This effectively eliminates the "hw" attribute
* from the description file
*/
if (!HW)
HW = (CU == NULL) || (DU && equal(DV, CU));
HO = host;
/*
If login script, verify access
*/
if (LI != NULL) {
if (*LI == '~')
(void) expand_tilde (&LI, NULL);
if (access (LI, F_OK | X_OK) != 0) {
printf("tip (warning): can't open login script \"%s\"\n", LI);
LI = NULL;
}
}
/*
If logout script, verify access
*/
if (LO != NULL) {
if (*LO == '~')
(void) expand_tilde (&LO, NULL);
if (access (LO, F_OK | X_OK) != 0) {
printf("tip (warning): can't open logout script \"%s\"\n", LO);
LO = NULL;
}
}
/*
* see if uppercase mode should be turned on initially
*/
if (cgetflag("ra"))
boolean(value(RAISE)) = 1;
if (cgetflag("ec"))
boolean(value(ECHOCHECK)) = 1;
if (cgetflag("be"))
boolean(value(BEAUTIFY)) = 1;
if (cgetflag("nb"))
boolean(value(BEAUTIFY)) = 0;
if (cgetflag("sc"))
boolean(value(SCRIPT)) = 1;
if (cgetflag("tb"))
boolean(value(TABEXPAND)) = 1;
if (cgetflag("vb"))
boolean(value(VERBOSE)) = 1;
if (cgetflag("nv"))
boolean(value(VERBOSE)) = 0;
if (cgetflag("ta"))
boolean(value(TAND)) = 1;
if (cgetflag("nt"))
boolean(value(TAND)) = 0;
if (cgetflag("rw"))
boolean(value(RAWFTP)) = 1;
if (cgetflag("hd"))
boolean(value(HALFDUPLEX)) = 1;
if (RE == NULL)
RE = (char *)"tip.record";
if (EX == NULL)
EX = (char *)"\t\n\b\f";
if (ES != NULL)
vstring("es", ES);
if (FO != NULL)
vstring("fo", FO);
if (PR != NULL)
vstring("pr", PR);
if (RC != NULL)
vstring("rc", RC);
if (cgetnum(bp, "dl", &DL) == -1)
DL = 0;
if (cgetnum(bp, "cl", &CL) == -1)
CL = 0;
if (cgetnum(bp, "et", &ET) == -1)
ET = 10;
}
static void
getremcap(char *host)
{
char **p, ***q, *bp, *rempath;
int stat;
rempath = getenv("REMOTE");
if (rempath != NULL) {
if (*rempath != '/')
/* we have an entry */
cgetset(rempath);
else { /* we have a path */
db_array[1] = rempath;
db_array[2] = _PATH_REMOTE;
}
}
if ((stat = cgetent(&bp, db_array, host)) < 0) {
if ((DV != NULL) ||
(host[0] == '/' && access(DV = host, R_OK | W_OK) == 0)) {
CU = DV;
HO = host;
HW = 1;
DU = 0;
if (!BR)
BR = DEFBR;
FS = DEFFS;
return;
}
switch (stat) {
case -1:
fprintf(stderr, "%s: unknown host %s\n", __progname,
host);
break;
case -2:
fprintf(stderr,
"%s: can't open host description file\n",
__progname);
break;
case -3:
fprintf(stderr,
"%s: possible reference loop in host description file\n", __progname);
break;
}
exit(3);
}
for (p = capstrings, q = caps; *p != NULL; p++, q++)
if (**q == NULL)
cgetstr(bp, *p, *q);
if (!BR && (cgetnum(bp, "br", &BR) == -1))
BR = DEFBR;
if (!LD && (cgetnum(bp, "ld", &LD) == -1))
LD = TTYDISC;
if (cgetnum(bp, "fs", &FS) == -1)
FS = DEFFS;
if (DU < 0)
DU = 0;
else
DU = cgetflag("du");
if (DV == NOSTR) {
fprintf(stderr, "%s: missing device spec\n", host);
exit(3);
}
if (DU && CU == NOSTR)
CU = DV;
if (DU && PN == NOSTR) {
fprintf(stderr, "%s: missing phone number\n", host);
exit(3);
}
if (DU && AT == NOSTR) {
fprintf(stderr, "%s: missing acu type\n", host);
exit(3);
}
HD = cgetflag("hd");
/*
* This effectively eliminates the "hw" attribute
* from the description file
*/
if (!HW)
HW = (CU == NOSTR) || (DU && equal(DV, CU));
HO = host;
/*
* see if uppercase mode should be turned on initially
*/
if (cgetflag("ra"))
setboolean(value(RAISE), 1);
if (cgetflag("ec"))
setboolean(value(ECHOCHECK), 1);
if (cgetflag("be"))
setboolean(value(BEAUTIFY), 1);
if (cgetflag("nb"))
setboolean(value(BEAUTIFY), 0);
if (cgetflag("sc"))
setboolean(value(SCRIPT), 1);
if (cgetflag("tb"))
setboolean(value(TABEXPAND), 1);
if (cgetflag("vb"))
setboolean(value(VERBOSE), 1);
if (cgetflag("nv"))
setboolean(value(VERBOSE), 0);
if (cgetflag("ta"))
setboolean(value(TAND), 1);
if (cgetflag("nt"))
setboolean(value(TAND), 0);
if (cgetflag("rw"))
setboolean(value(RAWFTP), 1);
if (cgetflag("hd"))
setboolean(value(HALFDUPLEX), 1);
if (cgetflag("dc"))
setboolean(value(DC), 1);
if (cgetflag("hf"))
setboolean(value(HARDWAREFLOW), 1);
if (RE == NOSTR)
RE = (char *)"tip.record";
if (EX == NOSTR)
EX = (char *)"\t\n\b\f";
if (ES != NOSTR)
vstring("es", ES);
if (FO != NOSTR)
vstring("fo", FO);
if (PR != NOSTR)
vstring("pr", PR);
if (RC != NOSTR)
vstring("rc", RC);
if (cgetnum(bp, "dl", &DL) == -1)
DL = 0;
if (cgetnum(bp, "cl", &CL) == -1)
CL = 0;
if (cgetnum(bp, "et", &ET) == -1)
ET = 10;
}
static void
getremcap(char *host)
{
int stat;
char tbuf[BUFSIZ];
static char buf[BUFSIZ/2];
char *bp = buf;
char **p, ***q;
if ((stat = rgetent(tbuf, host, sizeof (tbuf))) <= 0) {
if (DV ||
host[0] == '/' && access(DV = host, R_OK | W_OK) == 0) {
/*
* If the user specifies a device on the commandline,
* don't trust it.
*/
if (host[0] == '/')
trusted_device = 0;
CU = DV;
HO = host;
HW = 1;
DU = 0;
if (!BR)
BR = DEFBR;
FS = DEFFS;
RE = (char *)"tip.record";
EX = (char *)"\t\n\b\f";
DL = 0;
CL = 0;
ET = 10;
return;
}
(void) fprintf(stderr, stat == 0 ?
"tip: unknown host %s\n" :
"tip: can't open host description file\n", host);
exit(3);
}
for (p = capstrings, q = caps; *p != NULL; p++, q++)
if (**q == NULL)
**q = rgetstr(*p, &bp);
if (!BR && (BR = rgetnum("br")) < 0)
BR = DEFBR;
if ((FS = rgetnum("fs")) < 0)
FS = DEFFS;
if (DU < 0)
DU = 0;
else
DU = rgetflag("du");
if (DV == NOSTR) {
(void) fprintf(stderr, "%s: missing device spec\n", host);
exit(3);
}
if (DU && CU == NOSTR)
CU = DV;
if (DU && PN == NOSTR) {
(void) fprintf(stderr, "%s: missing phone number\n", host);
exit(3);
}
DB = rgetflag("db");
/*
* This effectively eliminates the "hw" attribute
* from the description file
*/
if (!HW)
HW = (CU == NOSTR) || (DU && equal(DV, CU));
HO = host;
/*
* see if uppercase mode should be turned on initially
*/
if (rgetflag("ra"))
boolean(value(RAISE)) = 1;
if (rgetflag("ec"))
boolean(value(ECHOCHECK)) = 1;
if (rgetflag("be"))
boolean(value(BEAUTIFY)) = 1;
if (rgetflag("nb"))
boolean(value(BEAUTIFY)) = 0;
if (rgetflag("sc"))
boolean(value(SCRIPT)) = 1;
if (rgetflag("tb"))
boolean(value(TABEXPAND)) = 1;
if (rgetflag("vb"))
boolean(value(VERBOSE)) = 1;
if (rgetflag("nv"))
boolean(value(VERBOSE)) = 0;
if (rgetflag("ta"))
boolean(value(TAND)) = 1;
if (rgetflag("nt"))
boolean(value(TAND)) = 0;
if (rgetflag("rw"))
boolean(value(RAWFTP)) = 1;
if (rgetflag("hd"))
boolean(value(HALFDUPLEX)) = 1;
if (rgetflag("hf"))
boolean(value(HARDWAREFLOW)) = 1;
if (RE == NULL)
RE = (char *)"tip.record";
if (EX == NULL)
EX = (char *)"\t\n\b\f";
if (ES != NOSTR)
(void) vstring("es", ES);
if (FO != NOSTR)
(void) vstring("fo", FO);
if (PR != NOSTR)
(void) vstring("pr", PR);
if (RC != NOSTR)
(void) vstring("rc", RC);
if ((DL = rgetnum("dl")) < 0)
DL = 0;
if ((CL = rgetnum("cl")) < 0)
CL = 0;
if ((ET = rgetnum("et")) < 0)
ET = 10;
}
void sphereBoundedScalarField( latticeMesh* mesh, scalar** field, char* fname ) {
// Allocate memory
*field = (scalar*)malloc( mesh->mesh.nPoints * sizeof(scalar) );
// Look for values
char* aux;
scalar in = 0,
out = 0;
scalar* centre;
scalar radius;
uint status,
n;
// Centre
if( vstring(&aux, "%s/internalField/centre", fname) )
status = lookUpVectorEntry("start/initialFields", aux, ¢re, &n);
// Radius
char* ins;
if( vstring(&ins, "%s/internalField/radius", fname) )
status = lookUpScalarEntry("start/initialFields", ins, 0, &radius);
// Inside value
if( vstring(&ins, "%s/internalField/inside", fname) )
status = lookUpScalarEntry("start/initialFields", ins, 0, &in);
// Outside value
if( vstring(&ins, "%s/internalField/outside", fname) )
status = lookUpScalarEntry("start/initialFields", ins, 0, &out);
if (status) {}
uint ii;
for( ii = 0 ; ii < mesh->mesh.nPoints ; ii++ ) {
scalar r = ( mesh->mesh.points[ii][0] - centre[0] ) * ( mesh->mesh.points[ii][0] - centre[0] ) +
( mesh->mesh.points[ii][1] - centre[1] ) * ( mesh->mesh.points[ii][1] - centre[1] ) +
( mesh->mesh.points[ii][2] - centre[2] ) * ( mesh->mesh.points[ii][2] - centre[2] );
if( r < radius * radius ) {
field[0][ii] = in;
}
else {
/* if(mesh->mesh.points[ii][1] <= (centre[1] - radius)) { */
/* field[0][ii] = in; */
/* } */
/* else{ */
field[0][ii] = out;
/* } */
}
}
free(aux);
free(centre);
}
