Sparsity.h

Go to the documentation of this file.

// Copyright (C) 2003 Johan Hoffman and Anders Logg.
// Licensed under the GNU GPL Version 2.

#ifndef __SPARSITY_H
#define __SPARSITY_H

#include <dolfin/NewArray.h>

namespace dolfin
{

  class ODE;
  class Vector;
  class Matrix;

  class Sparsity
  {
  public:
    
    // Constructor
    Sparsity(unsigned int N);
    
    // Destructor
    ~Sparsity();
    
    void clear();
    
    void clear(unsigned int i);

    void setsize(unsigned int i, unsigned int size);
    
    void set(unsigned int i, unsigned int j);
    
    void set(const Matrix& A);

    void transp(const Sparsity& sparsity);
    
    void detect(ODE& ode);

    inline bool sparse() const { return !pattern.empty(); }

    NewArray<unsigned int>& row(unsigned int i);

    const NewArray<unsigned int>& row(unsigned int i) const;

    void show() const;
    
  private:
    
    // Check given dependency
    bool checkdep(ODE& ode, Vector& u, real f0, unsigned int i, unsigned int j);

    // Number of components
    unsigned int N;
    
    // Increment of automatic detection of sparsity
    real increment;

    // The sparsity pattern
    NewArray< NewArray<unsigned int> > pattern;

  };

}

#endif

Documentation automatically generated with Doxygen on 10 Sep 2004.