[Okta Webinar] Learn how to a build a cloud-first strategyRegister Now

x
?
Solved

Reading and Writing Boost multi_array's to binary files

Posted on 2004-08-20
3
Medium Priority
?
653 Views
Last Modified: 2010-08-05
I am having problems writing a boost multi_array to a binary file and then reading it back.

I have a 3d multi_array which contains integers:

typedef boost::multi_array<int, 3> i3dArray;
typedef i3dArray::index iIndex;

Start Program
{
   i3dArray contactMatrix(boost::extents[NOATOMS][NOATOMS][NO_DISTANCES_BINS]);  // [167][167][150]

   loop through files of atomic structures
   {
        count the number of atom types i see at different distances and increment the relevant array positions
         i3dArray[i][j][k]++
    }

      // Verify values - this works (i.e. the i3dArray has correct values)
     for(iIndex i = 0; i != 50; ++i)
        for(iIndex j = 0; j != 50; ++j)
           for(iIndex k = 0; k != 50; ++k)
              cout << contact[i][j][k];

     // Then i want to write this object to a binary file.
    ofstream outfile( outname.c_str(), ios::out | ios::binary);
      outfile.write( (char*)&contact, sizeof(contact) );
    outfile.close();    

    // later i want to read in this object
    ifstream Source(file, ios::in | ios::binary );
    if (Source.fail())
      exit(1);
    Source.read( (char*)&contact, sizeof(contact) );
    Source.close();

      // Printing gives me a Segmentation Fault
     for(iIndex i = 0; i != 50; ++i)
        for(iIndex j = 0; j != 50; ++j)
           for(iIndex k = 0; k != 50; ++k)
              cout << contact[i][j][k];
}

I know this is a crude and incompetent method of writing an object but more correct ideas would be appreciated.
0
Comment
Question by:keepthursdayspecial
  • 2
3 Comments
 
LVL 86

Accepted Solution

by:
jkr earned 1500 total points
ID: 11851593
You need to access the multi_array's data area by calling the 'data()' meber function, i.e.

   ofstream outfile( outname.c_str(), ios::out | ios::binary);
     outfile.write( (char*)contact.data(), sizeof(int) * contacts.size());
   outfile.close();    

   ifstream Source(file, ios::in | ios::binary );
    Source.read( (char*)contact.data(), sizeof(int) * contacts.size() );
   Source.close();

Just taking the multi_array's address will not work, since this is not just 'flat' data...


0
 

Author Comment

by:keepthursdayspecial
ID: 11852059
Oh, I see.

However, this doesn't write anything to the binary file.
0
 

Author Comment

by:keepthursdayspecial
ID: 11852231
     string outname = "out/pdb"+pdbfiles[t]+".ent.chem";
      ofstream outfile( outname.c_str(), ios::out | ios::binary );
      outfile.write( (char*)contactMat.data(), sizeof(int)*contactMat.size() );
      outfile.close();

when i check the sizeof(int)*matrix i get 72 which is the same size as the output file.  

      cout << sizeof(int) << " " << sizeof(int) * contactMat.size() << " " << contactMat.size() << endl;

So 4 72 8.

The multi_array is [18] [18] [150] is this instance.

When I read the file back into a new multi_array of the same size i get an array of 0's.
0

Featured Post

Important Lessons on Recovering from Petya

In their most recent webinar, Skyport Systems explores ways to isolate and protect critical databases to keep the core of your company safe from harm.

Question has a verified solution.

If you are experiencing a similar issue, please ask a related question

Errors will happen. It is a fact of life for the programmer. How and when errors are detected have a great impact on quality and cost of a product. It is better to detect errors at compile time, when possible and practical. Errors that make their wa…
Container Orchestration platforms empower organizations to scale their apps at an exceptional rate. This is the reason numerous innovation-driven companies are moving apps to an appropriated datacenter wide platform that empowers them to scale at a …
The goal of the tutorial is to teach the user how to use functions in C++. The video will cover how to define functions, how to call functions and how to create functions prototypes. Microsoft Visual C++ 2010 Express will be used as a text editor an…
The viewer will learn how to clear a vector as well as how to detect empty vectors in C++.
Suggested Courses

872 members asked questions and received personalized solutions in the past 7 days.

Join the community of 500,000 technology professionals and ask your questions.

Join & Ask a Question