general_matrix_3x1.h

//
// Copyright (c) 2014 Limit Point Systems, Inc. 
//
// Licensed under the Apache License, Version 2.0 (the "License");
// you may not use this file except in compliance with the License.
// You may obtain a copy of the License at
//
//   http://www.apache.org/licenses/LICENSE-2.0
//
// Unless required by applicable law or agreed to in writing, software
// distributed under the License is distributed on an "AS IS" BASIS,
// WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
// See the License for the specific language governing permissions and
// limitations under the License.
//

// Interface for general_matrix_3x1 classes.

#ifndef GENERAL_MATRIX_3X1_H
#define GENERAL_MATRIX_3X1_H

#ifndef SHEAF_DLL_SPEC_H
#include "SheafSystem/sheaf_dll_spec.h"
#endif

#ifndef STD_IOSTREAM_H
#include "SheafSystem/std_iostream.h"
#endif

namespace fiber_bundle
{

template <typename T> class general_matrix_1x2;
template <typename T> class general_matrix_1x3;
template <typename T> class general_matrix_3x2;
template <typename T> class general_matrix_3x3;

template <typename T> class jcb_e13_row_dofs_type;

//==============================================================================
//==============================================================================

template <typename T>
class general_matrix_3x1
{
public:

  operator jcb_e13_row_dofs_type<T>& () const;

  static int number_of_rows();

  static int number_of_columns();

  static int d();

  T* operator[](int xrow);

  const T* operator[](int xrow) const;

  operator T* ();

  operator const T* () const;

  int row_index(int xrow) const;

  T components[3];

  //============================================================================
  //============================================================================

  void assign(const T& xvalue);

  void multiply(const T& xscalar, general_matrix_3x1<T>& xresult) const;

  general_matrix_3x1<T> multiply(const T& xscalar) const;

  void multiply(const general_matrix_1x2<T>& xother,
                general_matrix_3x2<T>& xresult) const;

  general_matrix_3x2<T> multiply(const general_matrix_1x2<T>& xother) const;

  void multiply(const general_matrix_1x3<T>& xother,
                general_matrix_3x3<T>& xresult) const;

  general_matrix_3x3<T> multiply(const general_matrix_1x3<T>& xother) const;

  void transpose(general_matrix_1x3<T>& xresult) const;

  general_matrix_1x3<T> transpose() const;


protected:
private:

};

// =============================================================================
// NON-MEMBER FUNCTIONS
// =============================================================================

template <typename T>
std::ostream& operator<<(std::ostream& xos, const general_matrix_3x1<T>& xm);


} // namespace fiber_bundle

#endif // ifndef GENERAL_MATRIX_3X1_H