コード例 #1
0
ファイル: util.c プロジェクト: josephcourtney/compass
static PyObject *calculate_MHD(PyObject *self, PyObject *args) {
    PyObject *list_a, *list_b;
    double **a = NULL;
    double **b = NULL;
    int a_len, b_len, a_dim, b_dim;
    if (!PyArg_ParseTuple(args, "OO", &list_a, &list_b)) {
        printf("Problem Parsing Arguments!\n");
        return NULL;
    }
    if( convert_2d_pylist( list_a, &a, &a_len, &a_dim ) != 0 ) {
        printf("Problem converting list #1\n");
        return NULL;
    }    

    if( convert_2d_pylist( list_b, &b, &b_len, &b_dim) != 0 ) {
        printf("Problem converting list #2\n");
        return NULL;
    }

    if( a_dim != b_dim) {
        printf("Dimension mismatch!\n");
        return NULL;
    }
    return Py_BuildValue("d", mhd( a, b, a_len, b_len, a_dim ) );
}
コード例 #2
0
ファイル: Atom.cpp プロジェクト: andreyto/imp-fork-proddl
Atom get_atom(Residue rd, AtomType at) {
  Hierarchy mhd(rd.get_particle());
  for (unsigned int i=0; i< mhd.get_number_of_children(); ++i) {
    Atom a(mhd.get_child(i));
    if (a.get_atom_type() == at) return a;
  }
  IMP_LOG(VERBOSE, "Atom not found " << at << std::endl);
  return Atom();
}
コード例 #3
0
ファイル: Atom.cpp プロジェクト: andreyto/imp-fork-proddl
Residue get_residue(Atom d, bool nothrow) {
  Hierarchy mhd(d.get_particle());
  do {
    mhd= mhd.get_parent();
    if (mhd== Hierarchy()) {
      if (nothrow) return Residue();
      else {
        IMP_THROW("Atom is not the child of a residue "  << d,
                  ValueException);
      }
    }
  } while (!Residue::particle_is_instance(mhd.get_particle()));
  Residue rd(mhd.get_particle());
  return rd;
}
コード例 #4
0
ファイル: main.cpp プロジェクト: chadvoegele/ledger-rest
int main(int argc, char** argv) {
  ledger_rest::args args(argc, argv);

  ledger_rest::stderr_logger logger(args.get_log_level());
  ledger_rest::ledger_rest_runnable ledger(args, logger);

  ledger_rest::mhd mhd(args, logger, ledger);

  std::list<ledger_rest::runnable*> runners{ &mhd, &ledger };
  ledger_rest::runner runner(logger, runners);

  ledger_rest::set_runner(&runner);
  if (std::signal(SIGINT, ledger_rest::stop_runner) == SIG_ERR) {
    logger.log(5, "Error setting signal handler.");
    exit(EXIT_FAILURE);
  }

  runner.run();
  return 0;
}