nextnano_postprocessing.c

Go to the documentation of this file.

/*
** nextnano_postprocessing.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 17 14:41:47 2009 Johnny Q. Hacker
** Last update Sun May 12 01:17:25 2002 Speed Blue
*/

#include <stdlib.h>
#include <glib.h>
#include <libnextnano/nextnano_postprocessing.h>
#include <libnextnano/read_nextnano_inputfile.h>
#include <libcommon/gcc_hints.h>

void nextnano_region_d1_free(struct nextnano_region_d1 *r) {
  if(r->material != NULL) {
    g_string_free(r->material, TRUE);
  }
  free(r);
}

void nextnano_region_free(struct nextnano_region *r) {
  if(r->storage != NULL) {
    switch(r->type){
    case NNR_D1:
      nextnano_region_d1_free(r->storage);
      break;
    }
  }
  free(r);
}

struct nextnano_region_d1* nextnano_region_d1_new() {
  struct nextnano_region_d1 *r;
  if(unlikely((r = (struct nextnano_region_d1*)malloc(sizeof(struct nextnano_region_d1))) == NULL)) return NULL;
  r->material = NULL;
  return r;
}

struct nextnano_region* nextnano_region_new(enum nextnano_region_type type) {
  struct nextnano_region *r;
  if(unlikely((r = (struct nextnano_region*)malloc(sizeof(struct nextnano_region))) == NULL)) return NULL;
  r->type = type;
  switch(type){
  case NNR_D1:
    r->storage = nextnano_region_d1_new();
    break;
  }
  if(unlikely((r->storage) == NULL)) {
    free(r);
  }
  return r;
}

/*ERROR TYPES
 *  0: Success
 *  1: Unsupported number of dimensions
 *  2: Failed to get new region
 *  3: failed to get new nextnano region.
 */
struct nanostructure* get_nanostructure(GHashTable* ifile, int* err) {
  int lerr;
  char *val;
  double *d;
  int *i;
  GList *regions = NULL;
  *err = 0;
  struct nanostructure *s;
  if(unlikely((s=nanostructure_new(NS_D1)) == NULL)) {
    *err = 2;
    return NULL;
  }
  struct nextnano_region *r;
  /*First, figure out what regions we have*/
  i = read_nextnano_extract_item_typed(ifile, "simulation-dimension", 1, "dimension", &lerr, INT_TYPE);
  if((*i) != 1) {
    *err = 1;
  }
  for(int c=0; (i = read_nextnano_extract_item_typed(ifile, "regions", 1, "region-number", &lerr, INT_TYPE)) != NULL; c++) {
    if(unlikely((r=nextnano_region_new(NNR_D1)) == NULL)) {
      *err = 3;
      /*Continue; we're going to be missing a region; let the caller
  decide what to do*/
      continue;
    }
    
  }
  /*Any unclaimed sections of the overall nanostructure size are in
    the default material.*/
  return s;
}