Eu gostaria de ser independente do software comercial para o meu trabalho científico. Acho que uma dependência de pacotes comerciais como o Matlab e suas caixas de ferramentas é insatisfatória, porque não sei se terei acesso ao Matlab no futuro e porque não gosto do idioma. Portanto, estou...
Considere com quase singular, o que significa que há um autovalor de que é muito pequeno. O critério de parada usual de um método iterativo é baseado no residual e diz respeito às iterações podem parar quando com n o número da iteração. Mas no caso que estamos considerando, pode haver um grande...
Suponhamos que tem uma malha que consiste em 2D triângulos não sobrepostos , e um conjunto de pontos { p i } M i = 1 ⊂ ∪ N k = 1 T K . Qual é a melhor maneira de determinar em qual triângulo cada um dos pontos se
Para fins de ensino, eu precisaria de uma função contínua de uma única variável "difícil" de aproximar com polinômios, ou seja, seria necessário um poder muito alto em uma série de potências para "encaixar" bem essa função. Pretendo mostrar aos meus alunos os "limites" do que pode ser alcançado com...
Sou um grande defensor do desenvolvimento orientado a testes em computação científica. Sua utilidade na prática é simplesmente impressionante e realmente alivia os problemas clássicos que os desenvolvedores de código conhecem. No entanto, existem dificuldades inerentes ao teste de códigos...
Gostaria de saber se existe uma maneira rápida de calcular a distância euclidiana de dois vetores na oitava. Parece que não há função especial para isso, então devo usar apenas a fórmula com
Estou executando uma tarefa de classificação de texto com R e obtenho uma matriz de termo de documento com tamanho 22490 por 120.000 (apenas 4 milhões de entradas diferentes de zero, menos de 1%). Agora, quero reduzir a dimensionalidade utilizando o PCA (Principal Component Analysis). Infelizmente,...
Comecei a usar python como a linguagem de programação para realizar todas as minhas tarefas no CFD. Eu tenho uma experiência muito pequena em programação. Sou formado em engenharia mecânica e estou cursando o ensino superior em engenharia aeroespacial. Às vezes, o aspecto computacional do CFD se...
O que mais detesto no MPI é lidar com tipos de dados (isto é, mapas de dados / máscaras) porque eles não se encaixam muito bem com o C ++ orientado a objetos. boost::mpino entanto, apenas suporta o MPI 1.1 no site: O boost :: mpi é uma interface amigável em C ++ para a Interface de transmissão...
Eu li sobre ciência da computação na Wikipedia, mas meu entendimento não é muito claro. A ciência computacional envolve programação? Quão diferente é a ciência computacional da computacional _ ____ , onde o espaço em branco pode ser qualquer disciplina (ciência dos materiais, engenharia, química,...
Ao calcular a fatoração QR na prática, utiliza-se as reflexões de Householder para zerar a parte inferior de uma matriz. Eu sei que, para calcular valores próprios de matrizes simétricas, o melhor que você pode fazer com as reflexões de Household é obtê-lo na forma tridiagonal. Existe uma maneira...
Na programação de cálculos de matriz densa, existe alguma razão para escolher um layout de linha principal do layout acima da coluna? Eu sei que, dependendo do layout da matriz escolhida, precisamos escrever o código apropriado para usar as memórias de cache efetivamente para fins de velocidade....
Provavelmente é uma questão de nível de estudante, mas não posso fazer isso exatamente para mim. Por que é mais preciso usar grades não uniformes nos métodos numéricos? Estou pensando no contexto de algum método de diferenças finitas para o PDE da forma . E suponha que eu esteja interessado em uma...
Como esta pergunta está relacionada à computação, decidi postar aqui. Espero que seja visto como apropriado. Acabei de começar a executar modelos atmosféricos e oceânicos e percebo que preciso de mais núcleos, memória e espaço em disco do que minha área de trabalho atual. Minha instituição possui...
Suponha minAvec(U)subject to Ui,j≤max{Ui,k,Uk,j},i,j,k=1,…,nminAvec(U)subject to Ui,j≤max{Ui,k,Uk,j},i,j,k=1,…,n\begin{align*} \min A &\mathrm{vec}(U) \\ &\text{subject to } U_{i,j} \leq \max\{U_{i,k}, U_{k,j}\}, \quad i,j,k = 1, \ldots, n \end{align*} onde UUU é uma matriz n × simétrica...
Existem razões para escolher um Runge Kutta implícito de alta ordem (IMRK) em vez do intervalo de tempo do BDF? BDF parece muito mais fácil para mim como estágio IMRK precisa linear resolve por passo de tempo. A estabilidade para BDF e IMRK parece ser um ponto discutível. Não consigo encontrar...
Estou usando o Valgrind + Callgrind para criar um perfil para o solucionador que escrevi. Conforme o manual do usuário do Valgrind, compilei meu código com as opções de depuração para o compilador: "Sem as informações de depuração, a melhor ferramenta que o Valgrind poderá fazer é adivinhar a...
Eu tenho um conjunto de dados e quero encontrar o parâmetro tal que minimize a soma isso é m k ∑ i = 1 | m - x i | .x1, x2, … , Xkx1,x2,...,xkx_{1}, x_{2}, \ldots, x_{k}mmm∑i = 1k∣∣m -
Tanto quanto sei, os solucionadores multigrid usam smoothers iterativos como Jacobi, Gauss-Seidel e SOR para atenuar o erro em várias frequências. Poderia ser utilizado um método de subespaço de Krylov (como gradiente conjugado, GMRES etc.)? Eu não acho que eles sejam classificados como...
Estou tentando resolver uma equação de Poisson 2D por diferenças finitas. No processo, obtenho uma matriz esparsa com apenas variáveis em cada equação. Por exemplo, se as variáveis fossem U , a discretização