NAMD
Functions
common.C File Reference
#include <unistd.h>
#include <errno.h>
#include <string.h>
#include <fcntl.h>
#include <sys/stat.h>
#include <ctype.h>
#include "common.h"
#include "Node.h"
#include "SimParameters.h"
#include "InfoStream.h"
#include "charm++.h"

Go to the source code of this file.

Functions

const char * NAMD_version_string ()
 
void NAMD_version (int *major, int *minor, int *patchLevel, int *releaseType)
 
int NAMD_gcd (int a, int b)
 
char * NAMD_stringdup (const char *s)
 
void NAMD_quit (const char *err_msg)
 
void NAMD_die (const char *err_msg)
 
void NAMD_err (const char *err_msg)
 
void NAMD_bug (const char *err_msg)
 
int NAMD_file_exists (const char *filename)
 
void NAMD_backup_file (const char *filename, const char *extension)
 
int NAMD_open_text (const char *fname, int append)
 
void NAMD_write (int fd, const char *buf, size_t count, const char *fname)
 
void NAMD_close (int fd, const char *fname)
 
FILE * Fopen (const char *filename, const char *mode)
 
int Fclose (FILE *fout)
 

Function Documentation

◆ Fclose()

int Fclose ( FILE *  fout)

Definition at line 435 of file common.C.

Referenced by ConfigList::ConfigList(), Ambertoppar::genclose(), PDB::PDB(), read_binary_file(), and Parameters::read_parameter_file().

436 {
437  int rc = -1;
438 #if !defined(NOCOMPRESSED)
439  rc = pclose(fout);
440 #endif
441  if (rc == -1) // stream not associated with a popen()
442  {
443  rc = fclose(fout);
444  }
445  return rc;
446 } /* Fclose() */

◆ Fopen()

FILE* Fopen ( const char *  filename,
const char *  mode 
)

Definition at line 341 of file common.C.

References endi(), and iout.

Referenced by ConfigList::ConfigList(), Ambertoppar::genopen(), GridforceFullMainGrid::initialize(), PDB::PDB(), read_binary_file(), and Parameters::read_parameter_file().

342 {
343  struct stat buf;
344  // check if basic filename exists (and not a directory)
345 
346 #if defined(NOCOMPRESSED)
347  if (!stat(filename,&buf))
348  {
349  FILE *rval;
350  while ( ! (rval = fopen(filename,mode)) ) {
351  if ( errno != EINTR ) break;
352  }
353  return(rval);
354  }
355 #else
356  if (!stat(filename,&buf))
357  {
358  if (!S_ISDIR(buf.st_mode)) {
359  FILE *rval;
360  while ( ! (rval = fopen(filename,mode)) ) {
361  if ( errno != EINTR ) break;
362  }
363  return(rval);
364  }
365  }
366  // check for a compressed file
367  char *realfilename;
368  char *command;
369  FILE *fout;
370  command = (char *)malloc(strlen(filename)+25);
371  // check for .Z (unix compress)
372  sprintf(command,"zcat %s.Z",filename);
373  realfilename = command+5;
374  iout << "Command = " << command << "\n" << endi;
375  iout << "Filename.Z = " << realfilename << "\n" << endi;
376  if (!stat(realfilename,&buf))
377  {
378  if (!S_ISDIR(buf.st_mode))
379  {
380  fout = popen(command,mode);
381  // on HP-UX, the first character(s) out of pipe may be
382  // garbage! (Argh!)
383  int C;
384  do
385  {
386  C = fgetc(fout);
387  // iout << "C is " << C << "\n" << endi;
388  if (isalnum(C) || isspace(C))
389  {
390  ungetc(C,fout);
391  C = -1; // outta loop
392  }
393  } while(C != -1);
394  free(command);
395  return(fout);
396  }
397  }
398  // check for .gz (gzip)
399  sprintf(command,"gzip -d -c %s.gz",filename);
400  realfilename = command+11;
401  iout << "Command = " << command << "\n" << endi;
402  iout << "Filename.gz = " << realfilename << "\n" << endi;
403  if (!stat(realfilename,&buf))
404  {
405  if (!S_ISDIR(buf.st_mode))
406  {
407  fout = popen(command,mode);
408  // on HP-UX, the first character(s) out of pipe may be
409  // garbage! (Argh!)
410  int C;
411  do
412  {
413  C = fgetc(fout);
414  // iout << "C is " << C << "\n" << endi;
415  if (isalnum(C) || isspace(C))
416  {
417  ungetc(C,fout);
418  C = -1; // outta loop
419  }
420  } while(C != -1);
421  free(command);
422  return(fout);
423  }
424  }
425  free(command);
426 #endif /* !defined(NOCOMPRESSED) */
427 
428  return(NULL);
429 } /* Fopen() */
std::ostream & endi(std::ostream &s)
Definition: InfoStream.C:54
#define iout
Definition: InfoStream.h:51

◆ NAMD_backup_file()

void NAMD_backup_file ( const char *  filename,
const char *  extension 
)

Definition at line 235 of file common.C.

References endi(), iERROR(), iout, and NAMD_file_exists().

Referenced by Controller::adaptTempInit(), colvarproxy_namd::backup_file(), ofstream_namd::open(), open_dcd_write(), Controller::outputExtendedSystem(), Controller::outputFepEnergy(), Controller::outputTiEnergy(), CollectionMaster::receiveDataStream(), and Controller::write_accelMDG_rest_file().

236 {
237  if (NAMD_file_exists(filename)) {
238  if ( ! extension ) extension = ".BAK";
239  char *backup = new char[strlen(filename)+strlen(extension)+1];
240  strcpy(backup, filename);
241  strcat(backup, extension);
242 #if defined(WIN32) && !defined(__CYGWIN__)
243  if ( remove(backup) ) if ( errno != ENOENT ) {
244  char *sys_err_msg = strerror(errno);
245  if ( ! sys_err_msg ) sys_err_msg = "(unknown error)";
246  iout << iERROR << "Error on removing file "
247  << backup << ": " << sys_err_msg << "\n" << endi;
248  fflush(stdout);
249  }
250 #endif
251  while ( rename(filename,backup) )
252  {
253  if ( errno == EINTR ) continue;
254  const char *sys_err_msg = strerror(errno);
255  if ( ! sys_err_msg ) sys_err_msg = "(unknown error)";
256  iout << iERROR << "Error on renaming file " << filename
257  << " to " << backup << ": " << sys_err_msg << "\n" << endi;
258  fflush(stdout);
259  if ( errno == EXDEV ) continue;
260  break;
261  // char errmsg[256];
262  // sprintf(errmsg, "Error on renaming file %s to %s",filename,backup);
263  // NAMD_err(errmsg);
264  }
265  delete [] backup;
266  }
267 }
std::ostream & endi(std::ostream &s)
Definition: InfoStream.C:54
int NAMD_file_exists(const char *filename)
Definition: common.C:219
#define iout
Definition: InfoStream.h:51
std::ostream & iERROR(std::ostream &s)
Definition: InfoStream.C:83

◆ NAMD_bug()

void NAMD_bug ( const char *  err_msg)

Definition at line 195 of file common.C.

References Node::Object().

Referenced by Sequencer::algorithm(), Controller::algorithm(), CollectProxyVectorInstance::append(), CudaComputeNonbonded::assignPatches(), CudaComputeNonbonded::assignPatchesOnPe(), PmePencilXYZ::backwardDone(), LdbCoordinator::barrier(), HomePatch::buildSpanningTree(), ReductionMgr::buildSpanTree(), ComputeNonbondedUtil::calc_error(), Controller::calcPressure(), ComputePmeMgr::chargeGridReady(), checkSum(), ofstream_namd::close(), Box< Patch, CompAtom >::close(), Controller::compareChecksums(), ComputeCUDAMgr::ComputeCUDAMgr(), ComputePmeCUDAMgr::ComputePmeCUDAMgr(), ComputeTclBC::ComputeTclBC(), CudaPmeTranspose::copyDataDeviceToDevice(), CudaPmeTranspose::copyDataDeviceToHost(), CudaPmeTranspose::copyDataHostToDevice(), WorkDistrib::createAtomLists(), ComputeCUDAMgr::createCudaComputeNonbonded(), ComputePmeCUDAMgr::createDevicesAndAtomFiler(), CsmAckHandler(), CsmHandler(), cudaNAMD_bug(), CudaPmeRealSpaceCompute::CudaPmeRealSpaceCompute(), MigrateAtomsCombinedMsg::distribute(), ComputeGridForce::doForce(), HomePatch::doGroupSizeCheck(), HomePatch::doGroupSizeCheck_SOA(), ComputePme::doWork(), CudaComputeNonbonded::doWork(), WorkDistrib::enqueueAngles(), WorkDistrib::enqueueAniso(), WorkDistrib::enqueueBonds(), WorkDistrib::enqueueCrossterms(), WorkDistrib::enqueueDihedrals(), WorkDistrib::enqueueExcls(), WorkDistrib::enqueueGromacsPair(), WorkDistrib::enqueueImpropers(), WorkDistrib::enqueueLCPO(), WorkDistrib::enqueuePme(), WorkDistrib::enqueueSelfA1(), WorkDistrib::enqueueSelfA2(), WorkDistrib::enqueueSelfA3(), WorkDistrib::enqueueSelfB1(), WorkDistrib::enqueueSelfB2(), WorkDistrib::enqueueSelfB3(), WorkDistrib::enqueueThole(), WorkDistrib::enqueueWork(), WorkDistrib::enqueueWorkA1(), WorkDistrib::enqueueWorkA2(), WorkDistrib::enqueueWorkA3(), WorkDistrib::enqueueWorkB1(), WorkDistrib::enqueueWorkB2(), WorkDistrib::enqueueWorkB3(), WorkDistrib::enqueueWorkC(), ScriptTcl::eval(), ComputeMap::extendPtrs(), findRecipEvirPe(), CudaComputeNonbonded::finishReductions(), ofstream_namd::flush(), PmeAtomStorage::getAtomIndex(), getBlockDim(), CudaPmeTranspose::getBuffer(), ComputeCUDAMgr::getComputeCUDAMgr(), ComputeCUDAMgr::getCudaComputeNonbonded(), ComputePmeCUDAMgr::getDevice(), ComputePmeCUDAMgr::getDevicePencilY(), ComputePmeCUDAMgr::getDevicePencilZ(), ParseOptions::getfromptr(), ComputePmeCUDAMgr::getHomePencil(), ComputeLjPmeSerialMgr::getLJparameters(), Communicate::getMessage(), ComputePmeCUDAMgr::getNode(), getPencilDim(), Controller::getTotalPotentialEnergy(), globalAtomIDToSOAID(), GridforceFullBaseGrid::grid_index(), FFTCompute::init(), ComputePmeCUDA::initialize(), ComputePatch::initialize(), ComputeHomePatch::initialize(), ComputePme::initialize(), LdbCoordinator::initialize(), ParallelIOMgr::initialize(), ComputePmeMgr::initialize(), ComputePmeMgr::initialize_computes(), IRSet::insert(), ParseOptions::issetfromptr(), ParseOptions::istruefromptr(), CudaComputeNonbonded::launchWork(), LdbCoordinator::LdbCoordinator(), ScriptTcl::load(), NamdState::loadStructure(), PatchMap::makePatches(), WorkDistrib::messageEnqueueWork(), CudaComputeNonbonded::messageEnqueueWork(), WorkDistrib::messageFinishCUDA(), WorkDistrib::messageFinishMIC(), HomePatch::minimize_rattle2(), Sequencer::multigratorPressure(), NAMD_version(), NAMD_write(), Node::Node(), PmeZPencil::node_process_untrans(), PmeYPencil::node_process_untrans(), ofstream_namd::open(), Box< Patch, CompAtom >::open(), CudaComputeNonbonded::openBoxesOnPe(), PatchMap::pack(), PatchMgr::PatchMgr(), CudaComputeNonbonded::patchReady(), PmePencilX::PmePencilX(), PmePencilXY::PmePencilXY(), PmePencilXYZ::PmePencilXYZ(), PmePencilY::PmePencilY(), PmePencilZ::PmePencilZ(), PmeTranspose::PmeTranspose(), ComputePmeMgr::pollChargeGridReady(), ComputePmeMgr::pollForcesReady(), PmePencilXYZMap::populateInitial(), PmePencilXMap::populateInitial(), PmePencilXYMap::populateInitial(), PmePencilMap::populateInitial(), Controller::printEnergies(), Rebalancer::printLoads(), PmePencilXMap::procNum(), PmePencilXYMap::procNum(), ProxyMgr::ProxyMgr(), HomePatch::rattle2(), Sequencer::reassignVelocities(), ProxyPatch::receiveAll(), ProxyPatch::receiveData(), ParallelIOMgr::receiveForces(), ParallelIOMgr::receivePositions(), ParallelIOMgr::receiveVelocities(), recursive_bisect_with_curve(), recvAck(), ComputePmeMgr::recvAck(), ComputePmeCUDADevice::recvAtoms(), ComputePmeCUDADevice::recvAtomsFromNeighbor(), recvBcast(), BroadcastMgr::recvBroadcast(), PatchMgr::recvCheckpointAck(), HomePatch::recvCheckpointLoad(), PatchMgr::recvCheckpointReq(), WorkDistrib::recvComputeMapChanges(), ComputeMgr::recvCudaGlobalMasterUpdateResultMsg(), GlobalMasterServer::recvData(), recvData(), ComputePmeCUDAMgr::recvDevices(), recvEvalCommand(), recvEvalResult(), PatchMgr::recvExchangeReq(), ComputePmeCUDADevice::recvForcesFromNeighbor(), ComputePmeMgr::recvGrid(), ComputeQMMgr::recvPartQM(), ComputePmeMgr::recvRecipEvir(), Output::recvReplicaDcdData(), ComputeGlobal::recvResults(), ProxyMgr::recvResults(), ComputePmeMgr::recvUngrid(), NodePmeMgr::recvUngrid(), ReductionSet::ReductionSet(), AtomMapper::registerIDsCompAtomExt(), AtomMapper::registerIDsFullAtom(), registerUserEventsForAllComputeObjs(), Node::reloadGridforceGrid(), ReductionMgr::remoteSubmit(), IRSet::remove(), replica_bcast(), replica_min_double(), CollectProxyVectorInstance::reset(), CollectionMgr::CollectVectorInstance::reset(), CollectionMaster::CollectVectorInstance::reset(), ComputeGlobal::saveTotalForces(), Sequencer::scalePositionsVelocities(), ComputeNonbondedUtil::select(), PmeZPencil::send_ungrid(), PatchMgr::sendAtoms(), NodePmeMgr::sendDataHelper(), ComputePmeMgr::sendDataPart(), Communicate::sendMessage(), PatchMgr::sendMovePatches(), PatchMgr::sendOneHomePatch(), NodePmeMgr::sendPencilsHelper(), ComputePmeMgr::sendPencilsPart(), CudaPmeTranspose::setDataPtrsYZX(), CudaPmeTranspose::setDataPtrsZXY(), ComputePmeCUDADevice::setPencilProxy(), ComputePmeCUDAMgr::setupPencils(), SimpleBroadcastObject< Tensor >::SimpleBroadcastObject(), Box< Patch, CompAtom >::skip(), CudaComputeNonbonded::skipPatchesOnPe(), CudaPmeKSpaceCompute::solve(), Node::startup(), ComputeMap::storeCompute(), ComputePmeCUDA::storePmeForceMsg(), PmePencilXYZ::submitReductions(), PmePencilZ::submitReductions(), ComputePmeMgr::submitReductions(), Sync::Sync(), TopoManagerWrapper::TopoManagerWrapper(), CudaPmeTranspose::transposeXYZtoYZX(), CudaPmeTranspose::transposeXYZtoZXY(), IRSet::unchecked_insert(), ComputePme::ungridForces(), ComputeMap::unpack(), GridforceGrid::unpack_grid(), CudaComputeNonbonded::unregisterBoxesOnPe(), AtomMapper::unregisterIDsCompAtomExt(), AtomMapper::unregisterIDsFullAtom(), ComputeCUDAMgr::update(), ComputeMgr::updateComputes(), Node::updateGridScale(), WorkDistrib::WorkDistrib(), AtomMapper::~AtomMapper(), and Rebalancer::~Rebalancer().

196 {
197  if ( ! err_msg ) err_msg = "(unknown error)";
198  const char *bug_msg =
199  "FATAL ERROR: See http://www.ks.uiuc.edu/Research/namd/bugreport.html";
200  CkPrintf("FATAL ERROR: %s\n%s\n",err_msg, bug_msg);
201  fflush(stdout);
202  char repstr[24] = "";
203  if (CmiNumPartitions() > 1 ||
204  (Node::Object()->simParameters != 0 &&
205  Node::Object()->simParameters->CUDASOAintegrateMode)) {
206  sprintf(repstr,"REPLICA %d ", CmiMyPartition());
207  // CkAbort ensures that all replicas die
208  CkAbort("%sFATAL ERROR: %s\n%s\n", repstr, err_msg, bug_msg);
209  }
210  CkError("%sFATAL ERROR: %s\n%s\n", repstr, err_msg, bug_msg);
211 #if CHARM_VERSION < 61000
212  CkExit();
213 #else
214  CkExit(2);
215 #endif
216 }
static Node * Object()
Definition: Node.h:86

◆ NAMD_close()

void NAMD_close ( int  fd,
const char *  fname 
)

Definition at line 315 of file common.C.

315  {
316 #ifdef WIN32
317  while ( _close(fd) ) {
318 #else
319  while ( close(fd) ) {
320 #endif
321  if ( errno != EINTR ) {
322  char errmsg[1024];
323  sprintf(errmsg, "Error on closing file %s", fname);
324  NAMD_err(errmsg);
325  }
326  }
327 }
void NAMD_err(const char *err_msg)
Definition: common.C:170

◆ NAMD_die()

void NAMD_die ( const char *  err_msg)

Definition at line 147 of file common.C.

References Node::Object().

Referenced by Controller::adaptTempInit(), ConfigList::add_element(), ComputeMsmMgr::addPotential(), after_backend_init(), Sequencer::algorithm(), Controller::algorithm(), AnisoElem::AnisoElem(), msm::PatchData::anterpolation(), msm::PatchData::anterpolationC1Hermite(), Parameters::assign_angle_index(), Parameters::assign_bond_index(), Parameters::assign_crossterm_index(), Parameters::assign_dihedral_index(), Parameters::assign_improper_index(), Parameters::assign_vdw_index(), WorkDistrib::assignNodeToPatch(), Molecule::build_go_arrays(), Molecule::build_go_params(), Molecule::build_go_sigmas(), Molecule::build_go_sigmas2(), Molecule::build_gridforce_params(), buildBondData(), HomePatch::buildRattleList(), ReductionMgr::buildSpanTree(), ParallelIOMgr::calcAtomsInEachPatch(), ComputeQMMgr::calcMOPAC(), ComputeQMMgr::calcORCA(), GlobalMasterTest::calculate(), GlobalMasterTcl::calculate(), GlobalMaster::calculate(), ComputeQMMgr::calcUSR(), Calloc(), GlobalMaster::check(), GroupRestraintParam::CheckParam(), ComputeMsmMgr::compute(), PmeRealSpace::compute_forces(), LjPmeRealSpace::compute_scaledForces(), AngleElem::computeForce(), LjPmeCompute::computeLJpotential(), ComputeTclBC::ComputeTclBC(), ConfigList::ConfigList(), NamdState::configListInit(), ConfigList::ConfigListNode::ConfigListNode(), Sequencer::correctMomentum(), Controller::correctMomentum(), WorkDistrib::createAtomLists(), cuda_errcheck(), cuda_getargs(), CudaComputeNonbonded::CudaComputeNonbonded(), cudaDie(), CudaGlobalMasterClient::CudaGlobalMasterClient(), CudaGlobalMasterServer::CudaGlobalMasterServer(), curandDie(), ComputeMsmMgr::d_stencil_1d(), Molecule::delete_qm_bonded(), ComputeGridForce::doForce(), ComputeNonbondedSelf::doForce(), HomePatch::doMarginCheck(), HomePatch::doMarginCheck_SOA(), Parameters::done_reading_files(), ComputeFullDirect::doWork(), ComputeTclBC::doWork(), ARestraint::EarlyExit(), Node::earlyExit(), GlobalMasterMisc::easy_calc(), msm::Array< msm::PatchDiagram >::elem(), msm::GridFixed< T, N >::elem(), msm::Grid< Vtype >::elem(), colvarproxy_namd::error(), ScriptTcl::eval(), PmeZPencil::fft_init(), PmeYPencil::fft_init(), PmeXPencil::fft_init(), LjPmeRealSpace::fill_charges(), PmeRealSpace::fill_charges(), Ambertoppar::get(), Molecule::get_atom_from_index_in_residue(), Molecule::get_atom_from_name(), Molecule::get_atomtype(), Parameters::get_dihedral_params(), Molecule::get_fep_bonded_type(), Parameters::get_improper_params(), Molecule::get_residue_size(), Parameters::get_vdw_params(), ComputeLjPmeSerialMgr::getLJparameters(), AngleElem::getMoleculePointers(), ImproperElem::getMoleculePointers(), CrosstermElem::getMoleculePointers(), DihedralElem::getMoleculePointers(), TholeElem::getMoleculePointers(), AnisoElem::getMoleculePointers(), BondElem::getMoleculePointers(), ExclElem::getMoleculePointers(), GromacsPairElem::getMoleculePointers(), TholeElem::getTupleInfo(), AnisoElem::getTupleInfo(), GlobalMasterIMD::GlobalMasterIMD(), GlobalMasterSymmetry::GlobalMasterSymmetry(), GromacsTopFile::GromacsTopFile(), BackEnd::init(), LjPmeCompute::initialize(), LjPmeMgr::initialize(), LdbCoordinator::initialize(), DeviceCUDA::initialize(), GridforceFullMainGrid::initialize(), GridforceFullSubGrid::initialize(), GridforceLiteGrid::initialize(), ComputeMsmMgr::initialize(), ComputePmeMgr::initialize(), SimParameters::initialize_config_data(), msm::PatchData::interpolation(), msm::PatchData::interpolationC1Hermite(), LjPmeKSpace::LjPmeKSpace(), ScriptTcl::load(), NamdState::loadStructure(), ludcmp(), main::main(), Malloc(), Node::mallocTest(), WorkDistrib::mapComputes(), HomePatch::minimize_rattle2(), Molecule::Molecule(), HomePatch::mollyAverage(), HomePatch::mollyMollify(), NAMD_new_handler(), NAMD_read_int(), NAMD_read_line(), NAMD_seek(), NAMD_write(), ComputeMsmMgr::ndsplitting(), LjPmeMgr::optimizeFFT(), Node::outputPatchComputeMaps(), MGridforceParamsList::pack_data(), WorkDistrib::patchMapInit(), PDB::PDB(), PDBUnknown::PDBUnknown(), Ambertoppar::preadln(), Molecule::prepare_qm(), ProblemParsing(), MsmBlockProxyMsg::put(), MsmC1HermiteBlockProxyMsg::put(), HomePatch::rattle1old(), HomePatch::rattle2(), Molecule::read_alch_unpert_angles(), Molecule::read_alch_unpert_bonds(), Molecule::read_alch_unpert_dihedrals(), read_binary_file(), Parameters::read_charmm_parameter_file(), Parameters::read_ener_table(), Parameters::read_energy_type(), Parameters::read_energy_type_bothcubspline(), Parameters::read_energy_type_cubspline(), Molecule::read_go_file(), Parameters::read_parameter_file(), Parameters::read_parm(), SimParameters::readExtendedSystem(), AmberParm7Reader::readparm(), GridforceFullBaseGrid::readSubgridHierarchy(), Realloc(), RecBisection::RecBisection(), Molecule::receive_GoMolecule(), Parameters::receive_Parameters(), CollectionMaster::receiveDataStream(), ParallelIOMgr::recvAtomsCntPerPatch(), Controller::recvCheckpointReq(), HomePatch::recvCheckpointReq(), ComputeMgr::recvComputeDPMEData(), ComputeMgr::recvComputeDPMEResults(), ComputeMgr::recvComputeEwaldData(), ComputeMgr::recvComputeEwaldResults(), ComputeMgr::recvComputeGlobalData(), ComputeMgr::recvComputeGlobalResults(), ComputeMsmSerialMgr::recvCoord(), ComputeExtMgr::recvCoord(), ComputeLjPmeSerialMgr::recvCoord(), ComputeFmmSerialMgr::recvCoord(), ComputeMgr::recvCudaGlobalMasterCreateMsg(), ComputeMgr::recvCudaGlobalMasterRemoveMsg(), ComputeMgr::recvCudaGlobalMasterUpdateMsg(), ComputeQMMgr::recvPartQM(), ComputeQMMgr::recvPntChrg(), Output::recvReplicaDcdData(), Node::reloadCharges(), Node::reloadGridforceGrid(), ReductionMgr::remoteRegister(), ReductionMgr::remoteUnregister(), ReductionSet::removeData(), Sequencer::rescaleaccelMD(), Controller::rescaleaccelMD(), ScriptTcl::run(), ComputeMsm::saveResults(), SimParameters::scriptSet(), ComputeNonbondedUtil::select(), Molecule::send_GoMolecule(), Parameters::send_Parameters(), ComputeMgr::sendComputeDPMEData(), ComputeMgr::sendComputeEwaldData(), StringList::set(), GroupRestraintParam::SetGroup1AtomFileIndices(), GroupRestraintParam::SetGroup1AtomListIndices(), GroupRestraintParam::SetGroup1RefPosition(), GroupRestraintParam::SetGroup2AtomFileIndices(), GroupRestraintParam::SetGroup2AtomListIndices(), GroupRestraintParam::SetGroupName(), msm::Array< msm::PatchDiagram >::setmax(), GroupRestraintParam::SetResCenter(), GroupRestraintParam::SetResDirection(), ComputeMsmMgr::setup_hgrid_1d(), ComputeMsmMgr::setup_periodic_blocksize(), SimParameters::setupIDWS(), PatchMap::sizeGrid(), ComputeMsmMgr::splitting(), ComputeMsmMgr::stencil_1d(), NamdCentLB::Strategy(), StringList::StringList(), ScriptTcl::tclmain(), TestArray_write(), TholeElem::TholeElem(), ComputeLonepairsCUDA::updateAtoms(), Node::updateGridScale(), ReductionMgr::willRequire(), and ReductionMgr::willSubmit().

148 {
149  if ( ! err_msg ) err_msg = "(unknown error)";
150  CkPrintf("FATAL ERROR: %s\n", err_msg);
151  fflush(stdout);
152  char repstr[24] = "";
153  if (CmiNumPartitions() > 1 ||
154  (Node::Object()->simParameters != 0 &&
155  Node::Object()->simParameters->CUDASOAintegrateMode)) {
156  sprintf(repstr,"REPLICA %d ", CmiMyPartition());
157  // CkAbort ensures that all replicas die
158  CkAbort("%sFATAL ERROR: %s\n", repstr, err_msg);
159  }
160  CkError("%sFATAL ERROR: %s\n", repstr, err_msg);
161 #if CHARM_VERSION < 61000
162  CkExit();
163 #else
164  CkExit(1);
165 #endif
166 }
static Node * Object()
Definition: Node.h:86

◆ NAMD_err()

void NAMD_err ( const char *  err_msg)

Definition at line 170 of file common.C.

References Node::Object().

Referenced by after_backend_init(), ComputeQMMgr::calcMOPAC(), ComputeQMMgr::calcORCA(), ComputeQMMgr::calcUSR(), close_dcd_write(), colvarproxy_namd::error(), NAMD_close(), NAMD_open(), NAMD_open_text(), NAMD_seek(), NAMD_write(), Controller::outputExtendedSystem(), PDB::PDB(), AmberParm7Reader::read_amber_parm_stage1(), ComputeQMMgr::recvPartQM(), Output::recvReplicaDcdData(), and PDB::write().

171 {
172  if ( ! err_msg ) err_msg = "(unknown error)";
173  const char *sys_err_msg = strerror(errno);
174  if ( ! sys_err_msg ) sys_err_msg = "(unknown error)";
175  CkPrintf("FATAL ERROR: %s: %s\n", err_msg, sys_err_msg);
176  fflush(stdout);
177  char repstr[24] = "";
178  if (CmiNumPartitions() > 1 ||
179  (Node::Object()->simParameters != 0 &&
180  Node::Object()->simParameters->CUDASOAintegrateMode)) {
181  sprintf(repstr,"REPLICA %d ", CmiMyPartition());
182  // CkAbort ensures that all replicas die
183  CkAbort("%sFATAL ERROR: %s: %s\n", repstr, err_msg, sys_err_msg);
184  }
185  CkError("%sFATAL ERROR: %s: %s\n", repstr, err_msg, sys_err_msg);
186 #if CHARM_VERSION < 61000
187  CkExit();
188 #else
189  CkExit(1);
190 #endif
191 }
static Node * Object()
Definition: Node.h:86

◆ NAMD_file_exists()

int NAMD_file_exists ( const char *  filename)

Definition at line 219 of file common.C.

References endi(), iERROR(), and iout.

Referenced by NAMD_backup_file().

219  {
220  int rval;
221  do {
222  rval = access(filename, F_OK);
223  } while ( rval != 0 && errno == EINTR );
224  if ( rval != 0 && errno != ENOENT ) {
225  const char *sys_err_msg = strerror(errno);
226  if ( ! sys_err_msg ) sys_err_msg = "(unknown error)";
227  iout << iERROR << "Error on checking file "
228  << filename << ": " << sys_err_msg << "\n" << endi;
229  fflush(stdout);
230  }
231  return ! rval;
232 }
std::ostream & endi(std::ostream &s)
Definition: InfoStream.C:54
#define iout
Definition: InfoStream.h:51
std::ostream & iERROR(std::ostream &s)
Definition: InfoStream.C:83

◆ NAMD_gcd()

int NAMD_gcd ( int  a,
int  b 
)

Definition at line 102 of file common.C.

References NAMD_gcd().

Referenced by Sequencer::integrate(), Sequencer::integrate_SOA(), and NAMD_gcd().

103 {
104  if(a == 0 || b == 0) return 0;
105  else if(a == b) return a;
106  else if(a > b) return NAMD_gcd(a-b,b);
107  return NAMD_gcd(a, b-a);
108 }
int NAMD_gcd(int a, int b)
Definition: common.C:102

◆ NAMD_open_text()

int NAMD_open_text ( const char *  fname,
int  append 
)

Definition at line 270 of file common.C.

References NAMD_err().

Referenced by ofstream_namd::open().

270  {
271  int fd;
272 
273  // open the file and die if the open fails
274 #ifdef WIN32
275  while ( (fd = _open(fname, O_WRONLY|(append?O_APPEND:O_EXCL)|O_CREAT|O_TEXT,_S_IREAD|_S_IWRITE)) < 0) {
276 #else
277 #ifdef NAMD_NO_O_EXCL
278  while ( (fd = open(fname, O_WRONLY|(append?O_APPEND:O_TRUNC)|O_CREAT,
279 #else
280  while ( (fd = open(fname, O_WRONLY|(append?O_APPEND:O_EXCL)|O_CREAT,
281 #endif
282  S_IRUSR|S_IWUSR|S_IRGRP|S_IROTH)) < 0) {
283 #endif
284  if ( errno != EINTR ) {
285  char errmsg[1024];
286  sprintf(errmsg, "Unable to open text file %s", fname);
287  NAMD_err(errmsg);
288  }
289  }
290 
291  return fd;
292 }
void NAMD_err(const char *err_msg)
Definition: common.C:170

◆ NAMD_quit()

void NAMD_quit ( const char *  err_msg)

Definition at line 125 of file common.C.

References Node::Object().

Referenced by GlobalMasterIMD::get_vmd_forces(), and Controller::printTiming().

126 {
127  if ( ! err_msg ) err_msg = "(unknown error)";
128  CkPrintf("EXITING: %s\n", err_msg);
129  fflush(stdout);
130  char repstr[24] = "";
131  if (CmiNumPartitions() > 1
132  || Node::Object()->simParameters->CUDASOAintegrateMode) {
133  sprintf(repstr,"REPLICA %d ", CmiMyPartition());
134  // CkAbort ensures that all replicas die
135  CkAbort("%sEXITING: %s\n", repstr, err_msg);
136  }
137  CkError("%sEXITING: %s\n", repstr, err_msg);
138 #if CHARM_VERSION < 61000
139  CkExit();
140 #else
141  CkExit(0);
142 #endif
143 }
static Node * Object()
Definition: Node.h:86

◆ NAMD_stringdup()

char* NAMD_stringdup ( const char *  s)

Definition at line 111 of file common.C.

111  {
112  char *rs;
113 
114  if(!s)
115  return NULL;
116 
117  rs = new char[strlen(s) + 1];
118  strcpy(rs,s);
119 
120  return rs;
121 }

◆ NAMD_version()

void NAMD_version ( int *  major,
int *  minor,
int *  patchLevel,
int *  releaseType 
)

Definition at line 56 of file common.C.

References NAMD_ALPHA_RELEASE, NAMD_BETA_RELEASE, NAMD_bug(), NAMD_FINAL_RELEASE, NAMD_UNKNOWN_RELEASE, and NAMD_version_string().

57 {
58  int NAMD_major = 0;
59  int NAMD_minor = 0;
60  int NAMD_patchLevel = 0;
61  int NAMD_releaseType = NAMD_UNKNOWN_RELEASE;
62 
63  // have to parse the version string to determine numbers
64  const char *vstr = NAMD_version_string();
65  char dot[4] = ""; // expect to read "."
66  char release[8] = ""; // expect to read "" or "."
67  // or "alpha" or "beta" or "final" or "release"
68  char extra[4] = ""; // bad version string if we read an extra character
69  int n = sscanf(vstr, "%d%3[^0-9]%d%7[^0-9]%d%c",
70  &NAMD_major, dot, &NAMD_minor, release, &NAMD_patchLevel, extra);
71  if (n < 3 || n > 5) {
72  NAMD_bug("Unable to parse NAMD_VERSION from Makefile");
73  }
74  if (n >= 4) {
75  if (strcmp(release, "a") == 0 || strcmp(release, "alpha") == 0) {
76  NAMD_releaseType = NAMD_ALPHA_RELEASE;
77  }
78  else if (strcmp(release, "b") == 0 || strcmp(release, "beta") == 0) {
79  NAMD_releaseType = NAMD_BETA_RELEASE;
80  }
81  else if (strcmp(release, ".") == 0 ||
82  strcmp(release, "final") == 0 ||
83  strcmp(release, "release") == 0 ||
84  strcmp(release, "patch") == 0) {
85  NAMD_releaseType = NAMD_FINAL_RELEASE;
86  }
87  else {
88  NAMD_bug("Unable to parse NAMD_VERSION from Makefile");
89  }
90  }
91  else {
92  NAMD_releaseType = NAMD_FINAL_RELEASE;
93  }
94  // now return values to caller
95  if (major != NULL) *major = NAMD_major;
96  if (minor != NULL) *minor = NAMD_minor;
97  if (patchLevel != NULL) *patchLevel = NAMD_patchLevel;
98  if (releaseType != NULL) *releaseType = NAMD_releaseType;
99 }
const char * NAMD_version_string()
Definition: common.C:51
void NAMD_bug(const char *err_msg)
Definition: common.C:195

◆ NAMD_version_string()

const char* NAMD_version_string ( void  )

Copyright (c) 1995, 1996, 1997, 1998, 1999, 2000 by The Board of Trustees of the University of Illinois. All rights reserved.

Definition at line 51 of file common.C.

Referenced by NAMD_version().

51  {
52  return NAMD_VERSION; // macro defined in toplevel Makefile
53 }

◆ NAMD_write()

void NAMD_write ( int  fd,
const char *  buf,
size_t  count,
const char *  fname 
)

Definition at line 295 of file common.C.

295  {
296  while ( count ) {
297 #if defined(WIN32) && !defined(__CYGWIN__)
298  long retval = _write(fd,buf,count);
299 #else
300  ssize_t retval = write(fd,buf,count);
301 #endif
302  if ( retval < 0 && errno == EINTR ) retval = 0;
303  if ( retval < 0 ) {
304  char errmsg[1024];
305  sprintf(errmsg, "Error on writing to file %s", fname);
306  NAMD_err(errmsg);
307  }
308  if ( retval > count ) NAMD_bug("extra bytes written in NAMD_write()");
309  buf += retval;
310  count -= retval;
311  }
312 }
void NAMD_err(const char *err_msg)
Definition: common.C:170
void NAMD_bug(const char *err_msg)
Definition: common.C:195