コード例 #1
0
        public void InitializeCSR()
        {
            // list non-anchored nodes and give them sequential ids
            activeNodes = nodes.FindAll(nd => !nd.anchored);
            int id = 0;

            foreach (Node nd in activeNodes)
            {
                nd.altId = id++;
            }

            // in each node make a list of elements to which it belongs
            foreach (Element elem in elems)
            {
                foreach (Node nd in elem.vertices)
                {
                    nd.AddElem(elem);
                }
            }

            // find neighbour nodes
            Parallel.ForEach(nodes, nd => nd.InferConnectivityInformation());

            // count total number of neighbours in all nodes
            int count = 0;

            foreach (Node nd in activeNodes)
            {
                count += nd.neighbors.Count;
            }

            // allocate CSR
            // each neighbor contributes 3 rows and 3 columns to CSR matrix, so nnz = count * 9
            // the size of the matrix is (number of active nodes)*(3 coordinates)
            csr = new CSR_System(activeNodes.Count * 3, count * 9);

            // 3) create CSR indices
            count = 0;
            foreach (Node nd in activeNodes)
            {
                int row_nnz = nd.CreateCSRIndices(count, csr.cols);
                csr.rows[nd.altId * 3]     = count;
                csr.rows[nd.altId * 3 + 1] = count + row_nnz;
                csr.rows[nd.altId * 3 + 2] = count + row_nnz * 2;
                count += row_nnz * 3;
            }
        }
コード例 #2
0
ファイル: Form1.cs プロジェクト: rhazel/FiniteElements
        private void backgroundWorker1_DoWork(object sender, System.ComponentModel.DoWorkEventArgs e)
        {
            Stopwatch sw = new Stopwatch();

            // initialize Paralution
            CSR_System.initParalution();
            do
            {
                mre.WaitOne(Timeout.Infinite);
                sw.Restart();
                model.AssembleLinearSystem();               // assemble
                msAssemble = sw.ElapsedMilliseconds;
                sw.Restart();
                model.Solve(useGPU, useDoublePrecision);    // solve
                sw.Stop();
                msSolve = sw.ElapsedMilliseconds;
                backgroundWorker1.ReportProgress(0);
            } while (true);

            // stop paralution
        }