source: ThirdParty/mpqc_open/src/lib/util/psi3/libpsio/tocread.cc@ 00f983

Action_Thermostats Add_AtomRandomPerturbation Add_RotateAroundBondAction Add_SelectAtomByNameAction Adding_Graph_to_ChangeBondActions Adding_MD_integration_tests Adding_StructOpt_integration_tests AutomationFragmentation_failures Candidate_v1.6.0 Candidate_v1.6.1 ChangeBugEmailaddress ChangingTestPorts ChemicalSpaceEvaluator Combining_Subpackages Debian_Package_split Debian_package_split_molecuildergui_only Disabling_MemDebug Docu_Python_wait EmpiricalPotential_contain_HomologyGraph_documentation Enable_parallel_make_install Enhance_userguide Enhanced_StructuralOptimization Enhanced_StructuralOptimization_continued Example_ManyWaysToTranslateAtom Exclude_Hydrogens_annealWithBondGraph FitPartialCharges_GlobalError Fix_ChronosMutex Fix_StatusMsg Fix_StepWorldTime_single_argument Fix_Verbose_Codepatterns ForceAnnealing_goodresults ForceAnnealing_oldresults ForceAnnealing_tocheck ForceAnnealing_with_BondGraph ForceAnnealing_with_BondGraph_continued ForceAnnealing_with_BondGraph_continued_betteresults ForceAnnealing_with_BondGraph_contraction-expansion GeometryObjects Gui_displays_atomic_force_velocity IndependentFragmentGrids_IntegrationTest JobMarket_RobustOnKillsSegFaults JobMarket_StableWorkerPool JobMarket_unresolvable_hostname_fix ODR_violation_mpqc_open PartialCharges_OrthogonalSummation PythonUI_with_named_parameters QtGui_reactivate_TimeChanged_changes Recreated_GuiChecks RotateToPrincipalAxisSystem_UndoRedo StoppableMakroAction Subpackage_levmar Subpackage_vmg ThirdParty_MPQC_rebuilt_buildsystem TremoloParser_IncreasedPrecision TremoloParser_MultipleTimesteps Ubuntu_1604_changes stable
Last change on this file since 00f983 was 860145, checked in by Frederik Heber <heber@…>, 8 years ago

Merge commit '0b990dfaa8c6007a996d030163a25f7f5fc8a7e7' as 'ThirdParty/mpqc_open'
  • Property mode set to 100644
File size: 2.3 KB
Line 
1/*!
2 \file tocread.c
3 \ingroup (PSIO)
4*/
5
6#include <unistd.h>
7#include <stdlib.h>
8#include <util/psi3/libpsio/psio.h>
9
10namespace psi3 {
11namespace libpsio {
12
13/*!
14** PSIO_TOCREAD(): Read the table of contents for file number 'unit'.
15**
16** \ingroup (PSIO)
17*/
18int psio_tocread(unsigned int unit)
19{
20 unsigned int i;
21 int errcod, stream, volume, entry_size;
22 psio_ud *this_unit;
23 psio_tocentry *last_entry, *this_entry;
24
25 this_unit = &(psio_unit[unit]);
26 entry_size = sizeof(psio_tocentry) - 2*sizeof(psio_tocentry *);
27
28 /* Check that this unit is actually open */
29 if(this_unit->vol[0].stream == -1) return(0);
30
31 /* Seek vol[0] to its beginning */
32 stream = this_unit->vol[0].stream;
33 errcod = lseek(stream, 0L, SEEK_SET);
34 if(errcod == -1) psio_error(unit,PSIO_ERROR_LSEEK);
35
36 /* Read the TOC from disk */
37 errcod = read(stream, (char *) &(this_unit->tocaddress.page),
38 sizeof(ULI));
39 if(errcod != sizeof(ULI)) return(1);
40 errcod = read(stream, (char *) &(this_unit->tocaddress.offset),
41 sizeof(ULI));
42 if(errcod != sizeof(ULI)) return(1);
43 errcod = read(stream, (char *) &(this_unit->toclen), sizeof(ULI));
44 if(errcod != sizeof(ULI)) return(1);
45
46 /* Malloc room for the TOC */
47 this_unit->toc = (psio_tocentry *) malloc(sizeof(psio_tocentry));
48 this_entry = this_unit->toc;
49 this_entry->last = NULL;
50 for(i=1; i < this_unit->toclen; i++) {
51 last_entry = this_entry;
52 this_entry = (psio_tocentry *) malloc(sizeof(psio_tocentry));
53 last_entry->next = this_entry;
54 this_entry->last = last_entry;
55 }
56 this_entry->next = NULL;
57
58 /* Seek the TOC volume to the correct position */
59 volume = (this_unit->tocaddress.page) % (this_unit->numvols);
60 errcod = psio_volseek(&(this_unit->vol[volume]), this_unit->tocaddress.page,
61 this_unit->tocaddress.offset, this_unit->numvols);
62 if(errcod == -1) psio_error(unit,PSIO_ERROR_LSEEK);
63
64 /* Read the TOC entry-by-entry */
65 this_entry = this_unit->toc;
66 for(i=0; i < this_unit->toclen; i++) {
67
68 /* This read() assumes a fixed ordering on the members of this_entry */
69 errcod = read(this_unit->vol[volume].stream, (char *) this_entry,
70 entry_size);
71 if(errcod != entry_size) psio_error(unit,PSIO_ERROR_READ);
72
73 this_entry = this_entry->next;
74 }
75
76 return(0);
77}
78
79}
80}
Note: See TracBrowser for help on using the repository browser.