/usr/include/eclib/mvector.h is in libec-dev 20160101-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 | // mvector.h: declarations of multiprecision integer vector class
//////////////////////////////////////////////////////////////////////////
//
// Copyright 1990-2012 John Cremona
//
// This file is part of the eclib package.
//
// eclib is free software; you can redistribute it and/or modify it
// under the terms of the GNU General Public License as published by the
// Free Software Foundation; either version 2 of the License, or (at your
// option) any later version.
//
// eclib 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 General Public License
// for more details.
//
// You should have received a copy of the GNU General Public License
// along with eclib; if not, write to the Free Software Foundation,
// Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301, USA
//
//////////////////////////////////////////////////////////////////////////
#if !defined(_MVECTOR_H)
#define _MVECTOR_H 1 //flags that this file has been included
#include "vector.h"
class vec_m {
friend class mat_m;
friend class msubspace;
public:
// constructors
vec_m(long n=0);
vec_m(long n, bigint* array);
vec_m(const vec_m&); // copy constructor
vec_m(const vec_i&); // copy constructor
vec_m(const vec_l&); // copy constructor
~vec_m(); // destructor
// member functions & operators
void init(long n=0); // (re)-initializes
vec_m& operator=(const vec_m&); // assignment
bigint& operator[](long i) const; // the i'th component
vec_m& operator+=(const vec_m&);
void addmodp(const vec_m&, const bigint&);
vec_m& operator-=(const vec_m&);
vec_m& operator*=(const bigint&);
vec_m& operator/=(const bigint&);
vec_m slice(long,long=-1) const; // returns subvec_m
vec_m operator[](const vec_i&) const; // subscript composition
vec_m operator[](const vec_l&) const; // subscript composition
void set(long i, const bigint& x); // sets v[i]=x
void add(long i, const bigint& x); // v[i]+=x
bigint sub(long i) const; // same as v[i] (no ref)
//converts to a vector of ints or longs if possible.
//the parameter here is a dummy just to distinguish these
vec_i shorten(int) const;
vec_l shorten(long) const;
// non-member (friend) functions and operators
friend long dim(const vec_m&); // the dimension
friend bigint operator*(const vec_m&, const vec_m&); // dot product
friend vec_m operator*(const mat_m& m, const vec_m& v);
friend int operator==(const vec_m&, const vec_m&);
friend int operator!=(const vec_m&, const vec_m&);
friend int trivial(const vec_m&); // v==zerovec_m?
friend ostream& operator<< (ostream&s, const vec_m&);
friend istream& operator>> (istream&s, vec_m&);
friend bigint mvecgcd(const vec_m&);
friend void swapvec(vec_m& v, vec_m& w);
friend int member(const bigint& a, const vec_m& v);//tests if a=v[i] for some i
friend mat_m restrict_mat(const mat_m& m, const msubspace& s);
friend mat_m prestrict(const mat_m& m, const msubspace& s, const bigint& pr);
// Implementation
private:
long d;
bigint * entries;
};
// Declaration of non-member, non-friend functions
vec_m operator+(const vec_m&); // unary
vec_m operator-(const vec_m&); // unary
vec_m operator+(const vec_m&, const vec_m&);
vec_m addmodp(const vec_m&, const vec_m&, const bigint&);
vec_m operator-(const vec_m&, const vec_m&);
vec_m operator*(const bigint&, const vec_m&); // componentwise
vec_m operator*(long, const vec_m&); // componentwise
vec_m operator/(const vec_m&, const bigint&); // componentwise
void makeprimitive(vec_m& v);
void elim(const vec_m& a, vec_m& b, long pos);
void elim1(const vec_m& a, vec_m& b, long pos);
void elim2(const vec_m& a, vec_m& b, long pos, const bigint& lastpivot);
vec_m express(const vec_m& v, const vec_m& v1, const vec_m& v2);
vec_m lift(const vec_m& v, const bigint& pr);
//lifts a mod-p vec_m to a rational and scales to a primitive vec in Z.
int liftok(vec_m& v, const bigint& pr);
//same, in place; returns success
bigint dotmodp(const vec_m& v1, const vec_m& v2, const bigint& pr);
// extra function definitions
long dim(const vec_m& v);
int operator!=(const vec_m& v, const vec_m& w);
vec_m operator+(const vec_m& v);
vec_m operator-(const vec_m& v);
vec_m operator+(const vec_m& v1, const vec_m& v2);
vec_m addmodp(const vec_m& v1, const vec_m& v2, const bigint& pr);
vec_m operator-(const vec_m& v1, const vec_m& v2);
vec_m operator*(const bigint& scal, const vec_m& v);
vec_m operator*(long scal, const vec_m& v);
vec_m operator/(const vec_m& v, const bigint& scal);
void makeprimitive(vec_m& v);
void elim(const vec_m& a, vec_m& b, long pos);
void elim1(const vec_m& a, vec_m& b, long pos);
void elim2(const vec_m& a, vec_m& b, long pos, const bigint& lastpivot);
#endif
|