00001 /* 00002 This file belongs to Aeneas. Aeneas is a GNU package released under GPL 3 license. 00003 This code is a simulator for Submicron 3D Semiconductor Devices. 00004 It implements the Monte Carlo transport in 3D tetrahedra meshes 00005 for the simulation of the semiclassical Boltzmann equation for both electrons. 00006 It also includes all the relevant quantum effects for nanodevices. 00007 00008 Copyright (C) 2007 Jean Michel Sellier <sellier@dmi.unict.it> 00009 00010 This program is free software; you can redistribute it and/or modify 00011 it under the terms of the GNU General Public License as published by 00012 the Free Software Foundation; either version 3, or (at your option) 00013 any later version. 00014 00015 This program is distributed in the hope that it will be useful, 00016 but WITHOUT ANY WARRANTY; without even the implied warranty of 00017 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 00018 GNU General Public License for more details. 00019 00020 You should have received a copy of the GNU General Public License 00021 along with this program. If not, see <http://www.gnu.org/licenses/>. 00022 */ 00023 00024 // Created on : 06 june 2007, Siracusa, Jean Michel Sellier 00025 // Last modified : 02 august 2007, Siracusa, Jean Michel Sellier 00026 00027 // Multiply a matrix in row-index sparse storage arrays sa and ija by a vector 00028 // x[1..n], giving a vector b[1..n] 00029 00030 void dsprsax(double xx[], double bb[], int n) 00031 { 00032 int i,k; 00033 00034 if(ija[1]!=n+2){ 00035 printf("dsprsax error : mismatched vector and matrix!"); 00036 exit(0); 00037 } 00038 for(i=1;i<=n;i++){ 00039 bb[i]=sa[i]*xx[i]; // start with diagonal term 00040 for(k=ija[i];k<=ija[i+1]-1;k++) // loop over off-diagonal terms 00041 bb[i]+=sa[k]*xx[ija[k]]; 00042 } 00043 }