read_nextnano_density.c

Go to the documentation of this file.

/*
** read_nextnano_density.c
** 
** Made by (Johnny Q. Hacker)
** Login   <solarion@borkborkbork>

    Copyright (C) 2009 Joseph Pingenot

    This program is free software: you can redistribute it and/or modify
    it under the terms of the GNU Affero General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.

    This program is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU Affero General Public License for more details.

    You should have received a copy of the GNU Affero General Public License
    along with this program.  If not, see <http://www.gnu.org/licenses/>.

** Started on  Tue Mar  3 13:05:05 2009 Johnny Q. Hacker
** Last update Sun May 12 01:17:25 2002 Speed Blue
*/

#include <stdio.h>
#include <libnextnano/read_nextnano_inputfile.h>
#include <libnextnano/read_nextnano_density.h>
#include <libcommon/textfile_common.h>

GList *read_nextnano_density_el(GFile* base, GHashTable *input) {
  int err;
  GError *error = NULL;
  GList *titles = g_list_append(NULL, "x");
  titles = g_list_append(titles, "total");
  titles = g_list_append(titles, "quantum");
  titles = g_list_append(titles, "classical");
  char *density_dir_str = read_nextnano_extract_item(input, "output-densities", 0, "destination-directory", &err);
  if(err != 0) {
    fprintf(stderr, "Error extracting output-densities destination-directory from nextnano input.\n");
    g_list_free(titles);
    return NULL;
  }
  GFile *density_dir = g_file_get_child(base, density_dir_str);
  GFile *densityfile = g_file_get_child(density_dir, "density1Del.dat");
  GFileInputStream *density_stream = g_file_read(densityfile, NULL, &error);
  if(error != NULL) {
    fprintf(stderr, "Error opening the density file for reading (error %s)\n", error->message);
    g_list_free(titles);
    return NULL;
  }
  GList *density = read_titled_whitespaced_data(G_INPUT_STREAM(density_stream), &error, &titles, '!', DOUBLE_TYPE);
  g_list_free(titles);
  return density;
}