FilesExpand file tree

 master

/

ComplexMatrixOperation.js

Copy path

Blame

More file actions

Blame

More file actions

Latest commit

 

History

History

History

311 lines (273 loc) · 7.82 KB

 master

/

ComplexMatrixOperation.js

Copy path

Top

File metadata and controls

• Code
• Blame

311 lines (273 loc) · 7.82 KB

Raw

Copy raw file

Download raw file

Open symbols panel

Edit and raw actions

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

120

121

122

123

124

125

126

127

128

129

130

131

132

133

134

135

136

137

138

139

140

141

142

143

144

145

146

147

148

149

150

151

152

153

154

155

156

157

158

159

160

161

162

163

164

165

166

167

168

169

170

171

172

173

174

175

176

177

178

179

180

181

182

183

184

185

186

187

188

189

190

191

192

193

194

195

196

197

198

199

200

201

202

203

204

205

206

207

208

209

210

211

212

213

214

215

216

217

218

219

220

221

222

223

224

225

226

227

228

229

230

231

232

233

234

235

236

237

238

239

240

241

242

243

244

245

246

247

248

249

250

251

252

253

254

255

256

257

258

259

260

261

262

263

264

265

266

267

268

269

270

271

272

273

274

275

276

277

278

279

280

281

282

283

284

285

286

287

288

289

290

291

292

293

294

295

296

297

298

299

300

301

302

303

304

305

306

307

308

309

310

311

var miscellaneous = require('./Miscellaneous');

var clone = require('matrix_deep_clone');

var elementMultiply = require('./ElementMultiply.js');

var complexMatrixOperations = function() {

function transpose(data) {

function stackColumns(data) {

var dim = miscellaneous.get_Dimensions(data);

var columns = dim[dim.length - 1];

var accessPattern = formPattern(data);

var newData = [];

for (var i = 0; i < columns; i++) {

var yak = [];

var counter = -1;

var temp = getColumns(data, i);

for (var v of accessPattern) {

yak[++counter] = temp[v];

}

newData = newData.concat(yak);

}

return newData;

}

function getColumns(data, index = 0) {

function innerGetColumns(inner_data, i) {

if (typeof(inner_data) == 'number' && i == index) {

newData[++counter] = inner_data;

} else {

for (var j in inner_data) {

innerGetColumns(inner_data[j], j);

}

}

}

var newData = [];

var counter = -1;

innerGetColumns(data, index);

return newData;

}

function formPattern(data) {

var lastAccessed = 0;

var counter = -1;

var temp = getColumns(data, 0);

var traverse = 0;

var yak = [];

var dim = miscellaneous.get_Dimensions(data);

var rows = dim[dim.length - 2];

var columns = dim[dim.length - 1];

for (var k = 0; k < temp.length; k++) {

yak[++counter] = traverse;

traverse = traverse + rows;

if (traverse >= temp.length) {

traverse = lastAccessed + 1;

lastAccessed++;

}

}

return yak;

}

var finalData = stackColumns(data);

var dim = miscellaneous.get_Dimensions(data).reverse();

return miscellaneous.reshape(finalData, dim);

}

function findDeterminant(data) {

function isSquare(data) {

var dim = miscellaneous.get_Dimensions(data);

var rows = dim[dim.length - 2];

var columns = dim[dim.length - 1];

if (rows == columns && dim.length == 2) {

return true;

} else {

return false;

}

}

function nCrossNDeterminant(data) {

function twoCrossTwoDeterminant(data) {

return (data[0][0] * data[1][1]) - (data[0][1] * data[1][0]);

}

var dim = miscellaneous.get_Dimensions(data);

var rows = dim[dim.length - 2];

var column = dim[dim.length - 1];

if (rows == 2 && column == 2) {

return twoCrossTwoDeterminant(data)

} else {

var toOperate = data[0];

var yui = -1;

var timepass = [];

var ans = 0;

var sign = 1;

for (i in toOperate) {

// starting j from one since at zero toOperate data is there

var newData = [];

var tig = -1;

for (j = 1; j < data.length; j++) {

var temp = [];

var counter = -1;

var yak = data[j];

for (v in yak) {

if (v != i) {

temp[++counter] = yak[v];

} else {

}

}

newData[++tig] = temp;

}

ans = ans + sign * toOperate[i] * nCrossNDeterminant(newData);

sign = sign * -1;

timepass[i] = newData;

}

return ans;

}

}

var dim = miscellaneous.get_Dimensions(data);

var rows = dim[dim.length - 2];

var columns = dim[dim.length - 1];

if (rows != columns) {

throw new Error("To calculate Determinant the matrix should be square");

} else {

var res = [];

var counter = -1;

function innerDeterminant(innerData) {

if (isSquare(innerData)) {

res[++counter] = nCrossNDeterminant(innerData);

} else {

for (i in innerData) {

innerDeterminant(innerData[i]);

}

}

}

innerDeterminant(data);

if (res.length == 1) {

return res[0];

} else {

var te = dim.slice(0, dim.length - 2);

if (te.length >= 2) {

return miscellaneous.reshape(res, te);

} else {

return res;

}

}

}

}

function generateIdentityMatrix(dim) {

var rowsNumber = dim[dim.length - 2];

var colNumber = dim[dim.length - 1];

function create2DMatrix(rows, columns) {

var res = [];

for (var i = 0; i < rows; i++) {

var temp = [];

for (var j = 0; j < columns; j++) {

if (i == j) {

temp[j] = 1;

} else {

temp[j] = 0;

}

}

res[i] = temp;

}

return res;

}

if (rowsNumber != colNumber) {

throw new Error("The dimensions should be square dimension");

} else {

if (dim.length == 2) {

return create2DMatrix(rowsNumber, colNumber);

} else {

// it is a higher dimension matrix

var res = [];

var subset = dim.slice(0, dim.length - 2);

for (var i = 0; i < subset.length; i++) {

var temp = [];

for (var j = 0; j < subset[subset.length - 1 - i]; j++) {

if (i == 0) {

temp[j] = create2DMatrix(rowsNumber, colNumber)

}

}

if (i == 0) {

res[i] = temp;

} else {

res[i] = clone.deepCloneMatrix(res[i - 1]);

}

}

if (res.length == 1) {

return res[0];

} else {

return res;

}

}

}

}

function matrixInverse(data) {

function findAdjoint(data) {

var dim = miscellaneous.get_Dimensions(data);

var rows = dim[dim.length - 2];

var column = dim[dim.length - 1];

var ansToSave = [];

var yui = -1;

var timepass = [];

var ans = [];

var sign = 1;

for (k in data) {

var toOperate = clone.deepCloneMatrix(data[k]);

for (i in toOperate) {

sign = Math.pow(-1, Number(k) + Number(i));

var newData = [];

var tig = -1;

for (j = 0; j < data.length; j++) {

if (j != k) {

var temp = [];

var counter = -1;

var yak = data[j];

for (v in yak) {

if (v != i) {

temp[++counter] = yak[v];

}

}

newData[++tig] = temp;

}

}

toOperate[i] = sign * findDeterminant(newData);

}

ansToSave[k] = toOperate;

}

return transpose(ansToSave);

}

function isSquare(data) {

var dim = miscellaneous.get_Dimensions(data);

var rows = dim[dim.length - 2];

var columns = dim[dim.length - 1];

if (rows == columns && dim.length == 2) {

return true;

}

return false;

}

function findInverse(data) {

if (findDeterminant(data) == 0) {

throw new Error("Determinant has to be non Zero");

}

var res = elementMultiply.ElementMultiply(findAdjoint(data), 1 / findDeterminant(data));

var dim = miscellaneous.get_Dimensions(res);

var rows = dim[dim.length - 2];

var columns = dim[dim.length - 1];

for (var i = 0; i < rows; i++) {

for (var j = 0; j < columns; j++) {

if (res[i][j] == -0) {

res[i][j] = 0;

}

}

}

return res;

}

function innerMatrixInverse(data) {

if (isSquare(data)) {

res[++counter] = findInverse(data);

} else {

for (i in data) {

innerMatrixInverse(data[i]);

}

}

}

var res = [];

var counter = -1;

innerMatrixInverse(data);

if (res.length == 1) {

return res[0];

} else {

return res;

}

}

return {

transpose: transpose,

findDeterminant: findDeterminant,

generateIdentityMatrix: generateIdentityMatrix,

matrixInverse: matrixInverse

}

}

module.exports = complexMatrixOperations();